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Preface 


In recent years, power converters have become an integral part of 
modern electronics, providing efficient and reliable energy con- 
version for a wide range of applications. With the advancement 
of digital control technology, there is a growing interest in using 
microcontrollers such as Arduino and STM32 to implement con- 
trol algorithms for power converters. 

This book aims to provide a comprehensive guide to under- 
standing and implementing digital control of power converters 
using Arduino and STM32 microcontrollers. It covers the funda- 
mental principles of the design and implementation of control 
algorithms, and practical examples of real-world applications. 

Whether you are an electrical engineering student, researcher, 
or practicing engineer, this book will help you develop the neces- 
sary skills to design and implement digital control systems for 
power converters. Through a combination of theoretical concepts 
and hands-on examples, you will gain a deep understanding of 
the principles and techniques required to successfully implement 
digital control systems in your own projects. 

The author hopes this book will serve as a valuable resource for 
anyone interested in learning about digital control of power con- 
verters using Arduino and STM32 microcontrollers. Happy read- 
ing and happy experimenting! 
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CHAPTER 1 


Digital Control with 
Arduino 


1.1 INTRODUCTION 


In this chapter, we'll be showing you how you can design and 


implement a digital controller on the Arduino to control or regu- 
late a continuous-time system. First, we'll start off by discussing 
the advantages of using digital control as opposed to the tradi- 
tional analog control and will then briefly go over the Z-transform 
and discretization of continuous-time systems. Then we'll move 
on to the six steps for designing and implementing a digital con- 
troller on the Arduino and lastly, we'll finish with designing a 
digital controller for a direct current (DC) motor speed control 
system. So, this chapter outline is as follows: 


Advantages of digital control 


Z-Transform 


Discretization of continuous-time systems (Tustin’s method) 


Implementation of discrete systems on the Arduino 
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* Digital controller design (Six steps) 


e DC motor speed control example: obtaining the transfer 
function 


e DC motor speed control example: controller design 


* DC motor speed control example: Arduino implementation 


Now keep in mind like we said in this chapter's description this is 
not an introductory course on signal processing or controls. If 
this is your first-time hearing words such as Laplace and 
Z-transforms or sampling stability, then you should first seek out 
this information elsewhere. This chapter is a practical tutorial 
series made to bridge the gap between the theory and practice of 
discrete-time systems and control theory. Again, our goal is to 
give you the tools needed to design a controller, how you will inte- 
grate this into your design is 10096 your responsibility to figure 
out. Therefore, the disclaimer for this chapter is as follows. 


* This is NOT an introductory course on LINEAR SYSTEMS, 
CONTROL THEORY, DIGITAL SIGNAL PROCESSING, 
ELECTRONICS, ETC. 


* This is a practical tutorial course that will allow students, 
engineers, and hobbyists to apply their known theoretical 
knowledge of control theory and discrete-time systems. 


e We will not hold ourselves responsible for any explanation 
of the basic fundamentals of linear systems and signals, 
Arduino coding, and basic electronics. It is required you 
have some exposure to these topics coming into reading this 
chapter. So, let's get on with this chapter. 


1.2 ADVANTAGES OF DIGITAL CONTROL 


Before we discuss why digital control is useful, let's take a look at 
analog control or continuous-time controllers. Here we have a 
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Driver Plant 


Error 


Setpoint Output 


Controller 


Sensor 


FIGURE 1.1 Analog or continuous-time control system. 


system block diagram of a control system where C(s) is our analog 
controller driving our plant, P(s), to a desired set point as shown 
in Figure 1.1. 

Now decades ago, even sometimes today C(s) was designed 
using purely discrete analog components typically consisting of 
operational amplifiers. Therefore, we can say the traditional feed- 
back systems (decades ago) relied on analog controllers or com- 
pensators for closed-loop control. Whereas electronic-based 
systems, they were mainly done using op-amps. A good example 
of this is a switch mode power supply topology known as the buck 
converter as depicted in Figure 1.2. 


losuag 


PWM Driver 


Desired 
V, 
Ò 9 output voltage 


Controller/Compensator 


FIGURE 1.2 The buck converter (step-down converter) as an electrical 
example. 
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A buck converter is a power electronic device whose job is to 
step down the input voltage to a desired output voltage to supply 
power to all sorts of electronic loads. The output voltage is 
sensed through H(s) and compared to a desired output voltage 
to then generate an error signal. This error signal is then fed to 
our controller or C(s) which generates and sends a continuous 
control voltage signal to the pulse-width modulation (PWM) 
driver driving the plant. The buck converter is duty cycle con- 
troller with the relation V, = D x V;, where D varies between 0 
and 1. The feedback controller C(s) drives the PWM duty cycle 
value to keep the error at zero. However, let's not overlook how 
it operates. The point we're trying to make is that this system is 
entirely analog controlled as there is no digital or discrete cou- 
pling between the blocks. Now if we take a look at the controller 
C(s), we can see that it's a type II compensator made out of just 
one op-amp and a few resistors and capacitors as illustrated in 
Figure 1.3. 

The inverting and non-inverting inputs of the op-amp are effec- 
tively the summing junction as seen in the feedback control system 
blocks. Again, don't get too caught up with understanding what's 
happening. The takeaway is that the analog control systems rely 


C3 Vout (0 


Ve (9 


FIGURE 1.3 The type II compensator (adds two poles and one zero into 
control loop). 
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only on analog components for their controller implementation 
made entirely out of analog components, forming a 10096 continuous- 
time coupled control loop and thus there are no time delays, sam- 
pling delays, and quantization errors and they are also easier to 
analyze and design using classical control techniques such as root 
locus or Bode plot analysis. Now that's all great. However, again 
why use digital control if analog controls seem so sufficient? Well 
with the use of analog control come some disadvantages. First, 
they're mainly composed of discrete analog components so if you 
wanted to tweak your controller or play around with various 
designs, this would mean physically replacing the components which 
is inconvenient if the components are soldered onto a board or not 
immediately available. Second, the passive components such as 
resistors and capacitors degrade and vary in value over time and 
temperature and thus the transfer function of your controller is 
never fixed, and the worst-case scenario is that if they are somehow 
damaged due to environmental temperatures, their values could 
have drifted so far off that they cause the controller to go unstable. 
Some systems incorporate multiple control loops so the more con- 
trollers there are the more space is required on a circuit board, and 
this can lead to bulkier designs and also increase the costs of the 
implementation. Some systems are also time variant or non-linear 
such that they require an adaptive or a non-linear controller. So, 
the controller's parameters will change in accordance to the varia- 
tion of the system in order for the control loop to stay stable or 
robust. However, this is not possible with analog controllers 
because the parameters set by discrete components are hard wired 
and fixed. Therefore, the disadvantages of analog/continuous-time 
controllers are summarized as the following: 


e Since they're composed of analog components (resistors, 
capacitors, op-amps, etc.), the parameters are fixed. 


e Analog components degrade over time and temperature, thus 
causing the overall transfer function to vary unintentionally. 
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e Some feedback systems require multiple control loops at 
high orders; this leads to requiring more board space and 
bulky designs (multiple op-amps and passives), and overall 
higher cost. 


e For systems whose dynamics are time-variants or non- 
linear, an adaptive or a non-linear controller may be 
required, and this is very difficult to achieve using discrete 
analog components considering the components are hard- 
wired (soldered) and fixed before deployment. 


So, if we go from an analog control scheme to a digital control 
scheme, you'll notice that we introduce a couple of new blocks. 
Unlike in the analog control scheme, the feedback signal is now 
sampled discreetly and fed into a microcontroller and then after 
the error signals are generated, a discrete control algorithm con- 
verts its error signal to a control signal which then leaves a micro- 
controller to drive the plant as shown in Figure 1.4. 

The C(z) is our discrete equivalent of our analog controller 
C(s) and can be implemented through firmware or embedded 
code and therefore we can easily adjust the coefficients that 
define the discrete transfer function C(z) and thus have more 
flexible control over our feedback system's behavior. So here are 
the main advantages of digital control. The first advantage is that 
it allows a designer or a control engineer to readily adjust the 


i Setpoint 


” P(s) > Output 


Plant/Process 


| Digital | Sensor 


Sensed output 


FIGURE 1.4 The digital control scheme blocks. 
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control parameters through firmware which is inexpensive, saves 
time, and provides a flexible degree of tuning. Second, since digi- 
tal controllers are typically implemented on a controller or 
microprocessor, they take up less board space and allow you to 
implement multiple control loops and add various features such 
as saturation or hysteresis with just a few lines of code, and, last, 
they allow for adaptive or non-linear control laws for systems 
that are time varying or non-linear in their transient or steady- 
state behavior. Therefore, the advantages of digital control can be 
summarized as follows: 


e Allows the designer or control engineer to readily adjust the 
controllers parameters via firmware - inexpensive, saves 
time, and provides a flexible degree on tuning. 


e Since digital controllers are typically implemented on a 
microcontroller or microprocessor, they take up less board 
space, allowing you to implement multiple control loops and 
add various features such as saturation and hysteresis with 
just a few lines of code. 


e Allows for adaptive or non-linear control laws for systems 
that are time variant or non-linear in their transient or 
steady-state behavior. 


So, this was just a brief motivation for why you should learn to 
implement digital controllers and their advantages over analog 
controllers. The rest of this chapter will focus on bridging the gap 
between the theory and the practical implementation in much 
greater detail. 


1.3 OVERVIEW OF THE ZTRANSFORM 


In this section, we will be briefly looking over the Z-transform and 
its role in the digital controllers implementation. In continuous- 
time feedback systems, we represent our controllers transfer func- 
tion C(s) in the s-domain which is obtained by taking the Laplace 
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E(s) U(s) 


+ oo 
Controller C(s) = I C (t)e —st dt 
— oo 


«-- Control Output 
C(s) = U(s) p 
E(s) «-- Error Input 


Laplace Transform of c(t) 


FIGURE 1.5 The Laplace transform of the controller. 


transform of the impulse response c(t), where s is a complex num- 
ber, o + jw, and w = 2nf. The famous Laplace transform, L{}, takes 
a signal or system in the time domain and brings it to its complex 
frequency-domain representation. Doing so, it converts differen- 
tial equations to algebraic expressions and the convolution opera- 
tion does multiplication and helps make analyzing systems easier 
as equations below and as shown in Figure 1.5. 


+0 


C(s)= [dear (1.1) 


C(s)= (s) (1.2) 


where C(s) is the Laplace transform of c(t), U(s) is the control out- 
put, and E(s) is the error input. 

In other words, it makes analyzing linear systems and designing 
controllers a lot easier. In the continuous-time or analog-domain, 
the controller C(s) is derived from a differential equation that 
describes the dynamics of the physical system where in most cases, 
it is the dynamics of an electrical analog circuit. For example, the 
type II compensator used for power supply regulation has the fol- 
lowing transfer function as shown in Figure 1.3: 


V.(s) 14 GRs 
C(s)= = 
6) Vals) | (&C + RC;)s- RG Cs (1.3) 
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To get the transfer function, we can apply the Kirchhoff's 
Current Law (KCL) on the red node in Figure 1.3. The inverting 
and non-inverting inputs of the ideal op-amp should have equal 
voltage, and the current flowing to the inverting point of ideal op- 
amp is equal to zero; therefore, the red node voltage is also equal 
to Veas depicted in Figure 1.3: 


Vef = V: Veer — V. Veet + Vice = Vout 


1 1 +0+ =0 (1.4) 
[sac pe R, Ry R 
Cs Cis 
1 
5 [nete aget [ect (15) 
Vout DUM de R, R 
Cs 
1+C,R 

ue lio ane 2 (1.6) 
Vout (RoC, + RC;)s + RC C,R,s 


Now don't worry too much about how this compensation was 
derived. Just understand the concept as understanding the spe- 
cific controller is way beyond the scope of this chapter. So, we 
cross multiply the terms, so that they are written as such: 


-Va(s)(1- C55) = V. (s)( (RC -RC;)s- RGCRSs) (1.7) 


Then, after taking the inverse Laplace transform or in other 
words going back to the time domain, we get the following dif- 
ferential equation: 


dV, (t) 
dt 


dV, (t d*V, 
=(R,C, + RC;) a HT ES ) 
(1.8) 


Vou (t)-C,R, 


Now keep in mind that when we take the Laplace transform or 
the inverse Laplace transform, you dont actually have to solve in 
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the integrals. You just look at the Laplace transform table pro- 
vided. So, this differential equation (1.8) is what describes physi- 
cally the relationship between the input and the output of the 
compensator. Now this dynamic is implemented using discrete 
components such as op-amp resistors and capacitors. So, how do 
we implement this type of dynamic digitally on a microcontroller 
for example? In order to answer that, we must briefly go over the 
Z-transform. Just like how the Laplace transform puts a continu- 
ous signal in terms of its 's' or its complex frequency domain. The 
Z-transform puts a discrete or digitized signal in terms of a z' or 
complex frequency domain or ‘z’ is equal to ‘e’ the power of ‘sT’ 
(eT). Now, C(z) is the Z-transform of c[k], where c[k] is a discrete- 
time signal sampled at our sampling time, T, as illustrated in 
Figure 1.6. 

So, the Z-transform is essentially the discrete version of the 
Laplace transform and it follows similar properties of the Laplace 
transform and it is a mathematical operation that takes a discrete- 
time signal and converts it to its complex frequency-domain 
representation: 


c(2)= Yelk)" (1.9) 


Where z= e 
Most importantly it converts convolution in the discrete-time 
domain to multiplication in the complex frequency domain. 


T=S ling ti 
ampling time «ti 


Oo 
t — k 


A-D Sampler 


| 


FIGURE 1.6 The Z-transform. 
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C(z) = E(z)U(z) Multiplication C(s) = E(s)U(s) Multiplication 
c[n] = e[n] «u[n] Convolution c(t) = e(t) «u(t) Convolution 
z-Domain s-Domain 


FIGURE 1.7 The Z-transform as a Laplace transform for discrete signals 
sampled at T. 


The Z-transform is essentially the “Laplace transform’ for discrete 
signals sampled at a fixed sample rate T as given in the following 
equations and is shown in Figure 1.7: 


c[n]- e[n ]*u[n]-» C(z) - E(z)U(z) (1.10) 
c(t) - e(t)*u(t)  C(s) - E(s)U(s) (1.11) 


It is also essential to know the mapping between the s-plain 
and the z-plane which comes from the following expression: 


z-et (1.12) 


s-o-jo (1.13) 


Now c is a real component and o is the imaginary component. 
It is apparent from this formula that the left half plane of the s- 
plane maps to inside the unit circle as shown in Figure 1.8, where 
a point can be represented by its magnitude, e”, and its angle, wT, 
as represented in the following equation: 


z=e =e el =e" T (1.14) 


So, if there was a point on the origin of the s-plane, it would 
map to 1 on the real axis of the z-plane since the magnitude is 1 
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Zeer! 
Irs} s=0+wj 
Im{z} 
oT Twj 
Z=e en $ 
[e] o m d Angle 
Magnitude 


Refs} p> Re{z} 


S-PLANE Z-PLANE (Unit Circle) 


FIGURE 1.8 Mapping between s-plane and z-plane. 


and the angle is 0, and if we were to move further negative in the 
real s-plane, this would correlate to moving closer to the origin 
of the z-plane However, if we move to the right of the s-plane, 
this would correlate to a magnitude greater than 1 and, thus, we 
would be outside the unit circle, and since we also have an 
imaginary component, the angle is non-zero. Now, if we keep 
the imaginary component the same and move the real compo- 
nent to the left half-plane, we will find ourselves back inside the 
unit circle in the z-plane, and, however, the angle will remain 
the same as represented in equations below and as shown in 
Figure 1.8: 


z=e" |-0=1<0° (1.15) 

z=e" p (en « 1) —0ZoT (1.16) 
AST _ (,0T 

£e" [atc >) Zor (1.17) 


This is too crucial to the digital controller design process but 
it’s good to know how the routes match between the Laplace and 
the Z-transform. A stable system should have all its poles within 
the unit circle. A pole outside the unit circle would indicate an 
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r[n] y(t) 


yin] T 
x 
Sampler y(t) 


Sensor 


FIGURE 1.9 Closed-loop feedback diagram of digital control. 


unstable system. So, taking a look at our closed-loop feedback dia- 
gram, a controller C(z) is now digital and is represented in the 
z-domain preceding it as a sampler which is sampling at a sam- 
pling time of T, and after it comes our zero-order hold (ZOH) 
block. We're just theoretically bringing our discrete signal back to 
the continuous-time domain to drive the plant as illustrated in 
Figure 1.9. 

Now, we'll go more into both of these blocks in later sections of 
this chapter. So, let's say we've designed C(z), that is great. Now 
how do we implement this? To answer that let's say C(z) equals the 
following transfer function which you can clearly see relates are 
discrete error input E(z) and our controller's discrete output U(z) 
as below equation: 


c(z)- U(z) bh *z'b*z?b; tz" by 


= (1.18) 
E(z) dbe d $e d th oa 


Now, if we cross multiply the turns, we get the following 
expression: 


U(z)(ao tz laz?a +...+z ay) 
= E(z)(by + z b 2 +27 by +...+z by) (1.19) 


Now if we distribute U(z) and E(z) on the respective sides of the 
equation and then we take the inverse Z-transform of the 
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expression, or in other words going back to the discrete-time 
domain, we get the following difference equation: 


Z'4X(z)}= OE (1.20) 


ayu[ n|+au[n-1]+a,u[ n-2]+...+ayu[n-—M 
= boe| n|+be[ n-1]+by,e[ n-2]+...4+bye[n-M | (1.21) 


This is called a difference equation which is essentially the dis- 
crete version of the differential equation and this equation is what 
is coded and used to develop the code in the firmware to imple- 
ment a digital controller in a real-time system. Now, notice the 'a' 
and ‘b’ coefficients in front of each discrete-time term in the dif- 
ferential equation case. These had to be using resistors and capaci- 
tor values because this is a difference equation and they can be 
implemented digitally. So, these coefficients can simply be adjusted 
in the code. Therefore, as the title of this section indicates, this is 
a very brief overview of the Z-transform; in the later sections, we 
will dive much deeper into making intuitive sense out of it and 
youll see where it’s relevant in the digital controller design 
process. 


1.4 DISCRETIZATION OF CONTINUOUS-TIME 
SYSTEMS 


In this section, we'll be going over the discretization of analog or 
continuous-time systems using Tustin's method. Whenever we 
have a pure continuous-time feedback loop with the continuous- 
time plant, sensor, and controller, we design our controller, C(s), 
using classical control techniques as shown in Figure 1.1. The two 
most popular techniques being root locus and the frequency 
response or, Bode plot analysis. Using the root locus, we observe 
how close loop poles move in the s-plane as a function ofthe gain 
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of our control or C(s). Using the frequency response analysis 
technique, we observe the Bode plot of the entire feedback loop 
and check the gain and phase margins to assess the stability and 
the transient performance. So, if a continuous-time control or 
C(s) is designed using these classical control techniques, how do 
we go about designing a digital controller, C(z), to drive a 
continuous-time plant such as a power supply, a motor, or a 
heater system. The best way to derive a digital controller is to first 
derive a continuous controller C(s) that can appropriately drive 
our plant and then knowing that z equals e to the sT (z = e), we 
can convert controller transfer function, C(s), to a discrete con- 
troller or C(z) such that it mimics the continuous-time control- 
ler’s behavior in the discrete/digital domain in regard to its 
frequency response as you can see in Figure 1.4. So essentially 
our goal is to design C(z) such that its frequency response matches 
closely with C(s) frequency response. So, let’s say we have a sys- 
tem or plant that we want to control with feedback and we design 
the following continuous-time controller or controller C(s) that 
has a pole in —5 and zero at —1 as represented in the following 
equation: 


c(s)- = 


$45 (1.22) 

So, we know that z equals e to the sT as below: 
z-e (1.23) 
If we arrange in the equation the ‘s’ as a function of Z’, we get 


that s equals a natural logarithm of z over T, therefore, we have the 
following equation: 


,In(z) (1.24) 
T 
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So, substituting for ‘s’ in C(s), we get C(z) that equals natural 
logarithm of ‘z’ over T plus 1 over natural logarithm of ‘z’ over T 
plus 5 as follows: 


T +1 
T 


This equation is suddenly a non-linear equation due to the nat- 
ural logarithm term and therefore we cannot take the inverse 
Z-transform of this, hence, we cannot work with this equation. 
So, what we'll need to do instead is approximate z equals e to the 
sT (eT) as a first-order linear equation. There are many approxi- 
mation methods that have been proposed but one in particular 
does a better job than the rest, and that's Tustin’s method which is 
also known as a bilinear transformation as given in the following 
equations: 


Euler's method: z =e" &14sT (1.26) 
1 
Backward difference: z =e” = (1.27) 
1-sT 
l+s— 
Tustin's method (bilinear transformation) z-e6 x » 
1-s— 
2 
(1.28) 


Now, if we map the left half plane of the s-plane to the z-plane 
using each of these approximations, you can see that Euler's 
method and the backward difference methods map very poorly in 
comparison to Tustin's method as illustrated in Figure 1.10. 


Digital Control with Arduino m 17 


Euler's method Backward difference Tustin's method 


Im(z) Im(z) Im(z) 


> Reiz) Me Reiz) > Re(z) 


FIGURE 1.10 Mapping of the left-hand s-plane (Re{s}< 0) to the z-plane 
(unit circle). 


Therefore, using Tustin's method, the mapped z-plane for val- 
ues Re{s} < 0 matches the closest to the original z-plane of z = e. 
So, if we want our digital controller's frequency response to match 
very closely to the frequency response of our continuous-time 
controller, it is in our best interest to go with Tustin's method. So, 
if we rearrange this new approximation of ‘s’ as a function of Z’, 
we get that ‘s’ equals 2 over T times Z minus 1 over Z plus 1 as the 
following equation: 


z= m 2 »s— (1.29) 


If we now substitute for 's' and C(s) and, for example, set the 
sampling time T to be 0.02 seconds, we get the following z-domain 
transfer function: 


+1 $= 
C(s)=-— T z4l (1.30) 
RES (sampling time) = 0.02 seconds 
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c(z)- 0.9619z —0.9429 (1.31) 
z —0.9048 

This transfer function C(z) is the discrete equivalent to the 
continuous-time transfer function C(s). Now we can use MATLAB 
to verify that their frequency responses match close enough. So, 
lets verify with MATLAB that C(z) has similar frequency response 
as C(s). Well use the code below, where the first code line is 
declaring our continuous-time transfer function C(s) and then to 
get the discrete transfer function C(z), reuse the ‘c2d’ function 
with the method set to ‘tustim and the sampling time, Ts, set to 
0.02. 


Cos c Ef CE 31;[1 515; $$ C(s) is the transfer 
function 
Ts = 0.02; $$ Sampling time is 0.02 s 
C_z = c2d(C s,Ts,'tustin'); $$ Discretizing C(s) 
using Tustin's 


method 
bode(C s); $$ Plotting the frequency 
response C(s) 
hold on; 
bode(C z); $$ Plotting the frequency 


response C(z) 


Then, we're plotting the Bode plots for both C(s) and C(z) on 
the same plot to see how they overlap. Now after running the code 
on MATLAB, we will get the figure showing the frequency 
responses of both C(s) and C(z), and they are matched very closely 
with one another as shown in Figure 1.11. 

Although they do not perfectly align, in our case of designinga 
digital controller that behaves like the continuous one, it'll do just 
fine. Hopefully, this section gave you enough information to know 
how to discretize a continuous-time system. In the next section, 
we will look into how a digital system like the one derived here is 
implemented in code on the Arduino Integrated Development 
Environment (IDE). 
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FIGURE 1.11 The Bode diagram of C(s) and C(z). 


1.5 ARDUINO IMPLEMENTATION OF 
DISCRETE-TIME SYSTEMS 

In this section, we will go over to see how you can implement a 
discrete-time system on the Arduino. Now we're going to step 
back from any discussion about feedback control and just see how 
we can implement any arbitrary discrete-time system on a micro- 
controller so let's say we have a second-order system, H(s), that we 
wish to discretize as below: 


20 
H(s)2 ————— (1.32) 
(s) s? 2.5s 4-20 


Now, if we apply a step input for this system which we can sim- 
ulate with MATLAB using the following code we can observe how 
the output behaves and we get a very typical under-damped 
second-order response as shown in Figure 1.12: 


20 = Digital Control of Power Converters 


Step response 
1.5 T T T T T 


Amplitude 


0 Tu za TS a. — 


0 1 2 3 4 5 6 
Time (seconds) 


FIGURE 1.12 The step response of second-order continuous-time sys- 
tem, H(s). 


H s = t£([20],[1 2.5 20]); 
step(H_s); 


We can also observe the frequency response of the system using 
the ‘bode’ function on MATLAB by the following code: 


H s = t£([20],[1 2.5 20]); 
bode (H_s) ; 


In the Bode plot as depicted in Figure 1.13, we will have a natu- 
ral or resonance frequency of about 4.07 rad/s. 

So, it’s a pretty slow and low bandwidth system and because we 
want to discretize this to implement on Arduino, we use the ‘c2d’ 
function on MATLAB as shown in the previous section and we 
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Bode diagram 
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FIGURE 1.13 The frequency response of the system, H(s). 


set the method to 'tustin' and the sampling-time to T which is 
equal to 0.1 seconds as the following code: 

Hs = t£([20], I1 2.5 201) ; 
T = DUE. 

H z = c2d(H s,T,'tustin'); 


This is the discrete transfer function we get as the equation 


2z-l 
s=2— 0,042372! - 0.084752 + 0.04237 
H(z)-H(s) Tz+1= z d (1.33) 
T =0.1sec z' —1.61z + 0.7797 


We'll call this discrete-system, H(z), and if we plot the transient 
response of both the analog and digital systems, you can see that 
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FIGURE 1.14 The transient response of both the analog and digital 
systems. 


their responses line up with each other as shown in Figure 1.14, 
and the MATLAB code is as below: 


t£([20],[1 2.5 20]); 
rl 

H z = c2d(H s,T,'tustin'); 
subplot (2,2,1); 

step(H s); 

subplot (2,2,2); 

step(H z); 

subplot (2,2,[3,41); 
step(H z, H s); 


Hs 
T 


O I 


Also, if we plot the frequency responses of both systems, we can 
see that they line up relatively well with some error at the higher 
frequency as shown in Figure 1.15, and it should be said that when 
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FIGURE 1.15 The frequency response of both the analog and digital systems. 
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it comes to design the controllers and digitizing systems, it isn't 
necessary that you get a controller that behaves exactly 10096 like 
its analog equivalent. 

The following MATLAB code is used for plotting the frequency 
responses of both systems: 


z-EEC[20] 5; [1-25:5-20]); 
= xus 

H z = c2d(H s,T,'tustin'); 
subplot (2,2,1); 

bode (H s); 

subplot (2,2,2); 

bode (H z); 
subplot (2,2, [3,4]); 

bode (H_z, H $); 


The goal is to design a controller that will robustly control the 
open-loop system. Now, we want to write this transfer function in 
its causal form with the following equation: 


2 (z) _ 0.04237z? + 0.08475z + 0.04237 
z’ —1.61z 4 0.7797 


(1.34) 


So, because it's a second-order system, we divide the numerator 
and the denominator by Z to get the following causal expression: 


-1 -2 
H(z)- Y(z) _ 0.04237 + 0.084752 + 0.042372 ass 
X(z) 1-1.61z"+0.7797z" 


Then as before we cross multiply the terms to separate X(z) and 
Y(z) as below: 


¥(z)(1-1.61z"' +0.7797z") 
= X(z)( 0.04237 +0.084752"' 0.042372") (1.36) 
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Y(z)-Y (z)r.eiz ' +Y (z)0.7797z? 
= X(z)0.04237+ X(2)0.08475z ' + X(z)0.04237z ^ (1.37) 


Also, after distributing X(z) and Y(z) on both sides, respec- 
tively, and then taking the inverse Z-transform to go back to the 
discrete-time domain, we get the following expression: 


y[n ]-1.61y[n-1 «0.7797 y[n-2 | 
= 0.04237x| n ]-0.08475x| n—1]--0.04237x|n-2] (1.38) 


What you'll notice over and over again is that z-! raised to the 
-1 correlates to the previous sample from one discrete-time step 
ago and that z2 raised to the —2 correlates to the sample from two 
discrete-time steps ago. So, we then take our prior output terms 
and move them over to the right side as the following equation: 


y|[n]- 0.04237x[n ] -0.08475x| n—1]+0.04237x| n-2 | 
*161y[n-1]-07797 y[n-2 ] (1.39) 


Now, we have a difference equation that is causal, meaning that 
its output only relies on the current and the previous input and 
output values. So, now let's go and implement this difference 
equation on the Arduino and see it in action. We'll first start off or 
we do know the script by first declaring variables and we're setting 
our sampling time to T equals 100 which is one hundred millisec- 
onds or 0.1 seconds in this case. We will then create two variables 
called last time and current. time, which we'll use to keep a con- 
stant sample time interval, 0.1 seconds, in every iteration, and 
then we'll have to declare our input and output variables for stor- 
ing the current and previous values as given by the following code: 


//Declaring variables 
int T - 100; //sampling time in milliseconds 
(T = 0.19) 
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unsigned long last time; 
unsigned long current time; 


float x; / / Xx [In] 
float y; / / Y in] 
float x 1; / / x [n- 1] 
float y 1; / / y [n-1] 
float x 2; / /x [n-2] 
float y 2; / /y [n-2] 


Also, for our void setup function which runs only once at the 
beginning ofthe program, we set all our variables to zero since we 
want our system response to have a zero initial condition in order 
to mimic a step response, we'll have to generate a unit step func- 
tion which is done here in the code below and the waveform of 
x[n] is illustrated in Figure 1.16: 


//Test Input Discrete Step Signal 
vaoid x_input () { 


if (current time == 1000){ //After 1000ms or 1 
second, step x[n] from 0 to 1 
x = 1.05; 
} 
} 
x[n] 
1.0 
— n 
T 
t=1sec 


FIGURE 1.16 Testing the input step signal. 
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The if statement condition is that once the elapsed time has 
reached 1000 milliseconds or 1 second, the input x[n] should be 
equal to 1. So essentially, we're going from x equals 0 to x equals 1, 
the moment we reach 1 second. Here, as illustrated in the code 
below, is the main void loop where we use a millis() function to 
keep track of the current elapsed time in milliseconds. We then 
compute the delta time or time difference between the current 
time and the previous time, and once this delta time value reaches 
our sample time of 0.1 seconds, we execute the if statement and 
then we save our current time as the last time. 


//Main Loop 
void loop() { 
current time = millis(); //stores the program 
running time (ms) 
int delta time = current time - last time;  // 
computes AT 
if (delta time >= T) { //when AT = T, read 
input x[n] and compute y[n] 
elapsed time - elapsed time « T/1000.0; // 
compute elapsed time in seconds 
x input(); //obtain input x[n] 
y = 0.0424*x + 0.0848*x 1 + 0.0424*x 2 + 
1.61*y 1 - 0.7797*y 2; 


X.2 - x1; //store x[n-1] as x[n-2] for 
next iteration 

y2^-2y1; //store y[n-1] as y[n-2] for 
next iteration 

x 1 = x; //store x[n] as x[n-1] for 
next iteration 

Y loy; //store y[n] as y[n-1] for 
next iteration 

last time - current time;  //restarts AT or 


delta time back to 0 


28 m Digital Control of Power Converters 


This method allows us to have artificial control over our 
sample time interval despite the higher sampling frequency of 
the Arduino’s ADC. So, in our ‘if’ statement, we will include 
our step input function, x input, and then our difference equa- 
tion. Notice that this line of code is exactly the one we derived 
before in Equation (1.31), and then lastly, we store the current 
values as the previous or delayed values for the next iteration, 
and there you have it as you can see in the above code. That is 
all you need to implement a discrete-time system on the 
Arduino. Let's run this code on the Arduino IDE and monitor 
the input and the output. So here, we have the Arduino code we 
just discussed. 


float T - 100.0; 

unsigned long last time; 

unsigned long current time; 

double elapsed time; 

double x; 

double y; 

double x 1; 

double y 1; 

double x 2; 

double y 2; 

void setup() { 
Serial.begin(9600); 


x = 0.0; 
y = 0.0; 
x 1 = 0.0; 
y_1 = 0.0; 
x2 = 0.0; 
y 2 = 0.0; 


void loop() { 
current time = millis(); 
int delta time - current time - last time; 
if (delta time >= T) { 
elapsed time = elapsed time + T/1000.0; 
x input(); 
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y = 0.0424*x + 0.0848*x 1 + 0.0424*x 2 + 
1.61*y 1 = 0.7797*y 2; 

x 2 = x ly 

y_2 = yl; 

x 1 x 

y_l = Y; 


Serial.print (elapsed_time) ; 
Serial.print (" DE 
Serial.println(y); 


last time - current time; 
} 
} 
void x_input() { 
if (current_time == 1000) { 
x = 1.0; 
} 
} 


One neat thing about the Arduino IDE is that you can moni- 
tor any value on the Serial Plotter, so we’ll want to monitor the 
step input x and the output y and see how they behave; hence, 
from Tools menu, we select the Serial Plotter option. So, let’s run 
the program as shown in Figure 1.17, we're getting the slow 


3.0 


0.0 


FIGURE 1.17 The slow under-damped response with Arduino IDES 
Serial Plotter. 
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under-damped response that we expected but in order to verify 
that this result is theoretically correct, let's confirm it with 
MATLAB Simulink. 

We'll first create a discrete transfer function block and type in 
the numerator and denominator coefficients in Equation (1.34), 
and then we'll set the sample time to 0.1 seconds as depicted in 
Figure 1.18. 

Next, we'll add the step input block which will step up from 0 
to 1 at 1 second with sample time equal to 0.1 seconds as illus- 
trated in Figure 1.19. 

Then lastly, we'll add a scope to plot and observe the output and 
the input and run a 10-second simulation as shown in Figure 1.20. 
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FIGURE 1.19 Adding the step input block. 
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FIGURE 1.20 Adding a scope and running the simulation. 
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FIGURE 1.21 The discrete transfer function step response in scope. 


From the scope, we can see that this is what we expect from the 
theoretical model of our discrete transfer function as shown in 
Figure 1.21. 

Now we want to compare this with the implementation result. 
So, we'll need to monitor the time and output data, so that we can 
create a two-dimensional variable on MATLAB. Now going back 
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to the Arduino, compile and then upload the program from the 
Tools menu, we can use the Serial Monitor to display the values 
numerically. By running the program, you can see the elapsed 
time and output columns on the Serial Monitor. 

Currently, we have completed the initial steps. Next, we will 
copy the values up to 10 seconds into Notepad++ and then import 
them into MATLAB. To do this, we will create a data variable that 
includes both the time and output data. Additionally, we will insert 
a row above the time value of 0.1 and set it to 0, as we want to 
include a row for time zero (time = 0 s), as illustrated in Figure 1.22. 

After that, we do rename the data variable as Output Data as 
shown in Figure 1.23. 

Then, we use the From Workspace block from Sources in 
Simulink to plot the experimental output data and then use a 
scope to view this signal alongside the theoretical step input. 
Therefore, we do double click on the From Workspace block and 
set Data as Output Data and Sample time to 0.1 seconds as 
depicted in Figure 1.24. 
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FIGURE 1.22 Creating a data variable that contains time and output 
data and adding a row above time 0.1 s. 
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FIGURE 1.24 Setting the From Workspace block properties. 


Then we connect it to a scope to compare the theoretical and 
experimental results, and if we were to plot the theoretical simu- 
lated output against the implementation output, we can see that 
they lie almost exactly on top of one another as shown in Figures 
1.25 and 1.26, respectively. 
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FIGURE 1.25 Theoretical simulated output against the implementation 
output. 
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FIGURE 1.26 Exact overlap of theoretical simulated output and imple- 
mentation output. 
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This goes to show that you can use the Arduino or any micro- 
controller to implement a difference equation or discrete-time 
system. So now that we're confident that we can implement a 
discrete-time system in the real time. Let’s move on to discussing 
how we can implement the feedback control algorithm on the 
Arduino. 


1.6 DIGITAL FEEDBACK CONTROL 
ALGORITHM 


In this section, we'll build upon what we learned in the previous 
section and now discuss how to implement feedback control 
structure on the Arduino. So, in the last section, we learned how 
to implement a discrete-time system such as C(z) in this case as 
shown in Figure 1.4. However, how do we code up the feedback 
control algorithm so that we can use an Arduino in a control loop 
to digitally regulate a continuous-time plant? Or, how do we 
implement the digital control scheme inside the dashed line as 
shown in Figure 1.4? In Figure 1.27, you'll see the portion of the 
feedback loop we wish to implement digitally and, in the code 
below, you'll see the main loop block that we discussed in the 


e[n] u[n] 


setpoint 


Controller 


X T Output y(t) 


output y[n] Sampler 


FIGURE 1.27 Digital implementation for the portion of the feedback 
loop. 


36 m Digital Control of Power Converters 


previous section for implementing the controller transfer func- 
tion, C(z), as the following equation: 


U(z) b,-bzi-bz?cr.cbz" 
C(z)= = -= 5 — (1.40) 
E(z) lta -a5z"-...ta,z 
void loop() { 
current time = millis(); //stores the program 
running time (ms) 
int delta time = current time - last time; // 


computes AT 
if (delta time >= T) { //when AT = T, read 
input x[n] and compute y[n] 
sensed output(); //obtain a current sample, 
output yin] 
e = setpoint - output y; //e[n] = setpoint 
- output yin] 


u = -a l*u 1- . -a n*u n + b 0*e + b m*e m; 
/ /C (z) 

un = u (n-1); //store u[n-(n-1)] as 
u[n-n] for next iteration 

em =e (m-1); //store e[n-(m-1)] as 


e[n-m] for next iteration 


u_l = u; //store u[n] as u[n-1] 


for next iteration 

e1l=e; //store e[n] as e[n-1] 
for next iteration 

last time = current time; //restarts AT or 
delta time back to 0 


So, let's work our way from the feedback path, so, if you take a 
look at the ‘if? statement in the code, you'll notice that it in itself is 
the sampler as it only executes every sample time of T. So, for 
example, if T equals one second, then the code inside the ‘if’ state- 
ment will iterate or update every one second. That being said if we 
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calla function which may be declared somewhere else in the script 
and reads in some discrete or analog sensed output data and if we 
place that in the ‘if? statement, then we're artificially sampling the 
output signal every T second. So, this code so far takes care of the 
sampling block. Next, we take our sensed output and then sub- 
tract it from a set point or a desired output to generate an error 
signal and this is done by doing a simple subtraction. Although 
we don't show it here, the variables we're introducing are declared 
at the beginning of the script. So, after we've obtained our error, 
we want to generate a control single to drive our plant. So next we 
implement the controller as a difference equation where the error 
eis our input and u is our control output. Here in the above code, 
we have it written as a pseudo code for the fallen transfer function 
with ‘@ and ‘b’ coefficients in Equation (1.30). If you notice this is 
exactly the same implementation we covered in the previous sec- 
tion. Nothing is different about the way it is implemented; all 
that's different is what we're feeding in as our input and getting 
back as our output. The difference equation structure is exactly 
the same. So, that is all, and this is all you need to implement the 
feedback control structure as shown in Figure 1.4. Now, you may 
wonder what is a zero-order hold (ZOH) block for as you can see 
in Figure 1.4, and why did we include it? Well let's step back and 
assess its theoretical function if we have a continuous signal, x(f), 
and pass it through a sampler which samples every T second, then 
we theoretically end up with the discrete signal, x[n]. Now, what 
the zero-order hold (ZOH) block does is hold the sample value for 
T seconds so that we get a staircase looking version ofthe continu- 
ous signal as depicted in Figure 1.28. 


x[n] x,(0) 


x(t) 
lis = n 
T 


Sampler ZOH 


a] 


FIGURE 1.28 The zero-order hold (ZOH) block operation. 
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ó[n] 


1 ó(t) 
< > 
7 n « 0 >t 
Discrete time Continuous time 
impulse function impulse function 


FIGURE 1.29 The impulse function is used only for theoretical deriva- 
tion and analysis. 


So, how do we implement a ZOH block? The answer is we don’t 
have to. It’s already intrinsically part of the firmware or code. 
Keep in mind that this signal x[n] does not exist in our physical 
reality and it is simply a mathematical representation of a sampled 
signal. Nowhere, inside the microcontroller, you will find a signal 
where the value only shows up in for testimony at the beginning 
of the sample period and then suddenly disappearing in between. 
This kind of representation is only useful for theoretical or math- 
ematical analysis. The same goes for the impulse function, as 
illustrated in Figure 1.29. 

They only exist for theoretical derivation and analysis purposes. 
Therefore, all the variables in the ‘if (delta time >= T) {F} state- 
ment are held constant for the entire T sampling period and only 
change or update in the next iteration. So, the ZOH block is intrin- 
sically present in the void loop() function, it is not a function we 
implement and it’s not something we have to implement directly. 
Again, the reason we include the sampler and the ZOH block in 
our diagram is because they mathematically and theoretically rep- 
resent what is happening to our signals and they are a model rep- 
resentation of what is happening to our signal, theoretically. So, 
there you have it and that is all you need to implement the follow- 
ing control structure as represented in the above code, and you can 
also see in Figure 1.30. 


setpoint 


Controller 


oS output_y(t) 
« 


output y[n] Sampler 


void loop() ( 


current time - millis(); 


int delta time - 


current time - last time; 


if (delta time »- T) ( 


sensed output(); 


e 
u 


u, 


e 


UL. 
e: 


1 
1 


setpoint - output y; 


-a 1*u 1 +. 


u (n-2); 
e (m-1); 


u; 
e; 


+ -a n*u n + b_O*e + .. 


last_time = current_time; 


} 


FIGURE 1.30 Arduino code for implementing the control structure block diagram. 
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So, now that we are in most of the theory out of the way. In the 
next few sections, we'll dive into the design step process and even- 
tually demonstrate this process using a DC motor speed control 
example. 


1.7 DIGITAL CONTROLLER DESIGN (SIX STEPS) 

In this section, we'll be going over the six design steps you can use 
to design a digital controller to control your continuous-time sys- 
tem, so, here are the six design steps. 


1. Obtain the plant’s transfer function, P(s) 
2. Choose an approximate sampling time, T 


3. Design controller, C(s), while taking into account the effects 
of ZOH/sampling 


4. Use bilinear transform/Tustin’s method to discretize C(s) to 
C(z) 


5. Take the inverse Z-transform of C(z) to obtain the difference 
equation 


6. Implement the difference equation on the Arduino and run 
the system! 


We want to let you know that this step-by-step procedure isn't the 
only way to go about designing a digital controller. This is just a 
method that makes sense to us and hopefully will make intuitive 
sense to you. So, the first step is to obtain the continuous-time 
transfer function of the system you wish to control and we will 
refer to this transfer function as P(s). Now this step is very open 
ended because there are many ways to obtain a transfer function 
of an open-loop system. For example, the method used to obtain 
the transfer function ofa power supply will differ from the method 
used to obtain the transfer function of a DC motor or a heater 
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system. So, for now, we'll assume you know what the open-loop 
transfer function is for the system you want to control. There's 
plenty of literature and papers out there. We'll show you how to 
get the transfer function of a quad-copter DC motor, heater sys- 
tem, power supply, etc. In the coming section, we'll be demon- 
strating this step for a DC motor but it is your job to figure out the 
transfer function of the system you wish to control. So, let's say 
you have your open-loop transfer function. So, the next step is to 
choose the sampling frequency or sampling time for your digital 
controller. So, there is actually no clear-cut method or equation to 
choose a sampling frequency or sampling time. If you Google 
this, you'll see that there are several ways to go about choosing a 
sampling time. However, Tustin's method works better when the 
sampling frequency (1/T) is much higher than the Nyquist fre- 
quency. So, what is the Nyquist frequency? Nyquist frequency (f,) 
is the minimum sampling frequency needed to prevent introduc- 
ing aliasing or introducing distortion/error in your closed-loop 
system. The Nyquist frequency is defined as being strictly twice 
the bandwidth of your open-loop system as the following 
equation: 


fa 22X fow or T, sw (1.41) 


But like we said Tustin’s method works better when the sam- 
pling frequency is much higher than the Nyquist frequency. So, 
as a rule of thumb is preferable to use a sampling frequency 
that's 5-10 times faster than the bandwidth (f,,,) of your open- 
loop system. The bandwidth of your system is typically defined 
as being the frequency at which the gain of your system is about 
3 dB is less than the DC gain. So, the bandwidth frequency or 3 
dB frequency is essentially the point where the DC gain begins 
to roll off, especially at 3 dB below the DC gain. So, in order to 
find the bandwidth of your system, it is required to know the 
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FIGURE 1.31 The bandwidth frequency or 3 dB frequency. 


transfer function so that you may plot the Bode plot and locate 
this minus 3 dB point on the magnitude graph as shown in 
Figure 1.31. 

Next, we move on to design the controller, C(s). So, as we 
mentioned in the earlier section when needing a digital control- 
ler, it's much easier to design it as a continuous mode controller 
and then discretize it to implement it digitally. So, we got the 
open-loop transfer function P(s) and we've chosen our sampling 
frequency so can we start the controller, C(s), design process 
yet? No, we cannot because well be designing a continuous- 
time controller using continuous-time design methods such as 
root locus or the frequency response analysis, we need to take 
into account the continuous-time effects the sampling and the 
ZOH function will have in the control loop as depicted in 
Figure 1.32. 

If we take a continuous-time signal x(t) and pass it through the 
sampler and the ZOH block, we get the expected continuous 
staircase signal which we'll call x(t), and if we were to overlap 
these signals in time you can see the hold delay, and it's 
discontinuous-time hold delay that we wish to model and include 
in our control loop before designing our controller as illustrated 
in Figure 1.33. 


setpoint r(t) y(t) 


To} 


So far we know P(s)... 


But can we design C(s) yet ? 


r[n] 


We MUST take into account ^ T 
Wa | JS 
the continuous-time effects 
of ZOH/sampling in the Sampler y(t) 


transmission loop! 


y(t) 


FIGURE 1.32 The continuous-time effects of ZOH/sampling in the transmission loop. 
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x(t) x(t) 


3l 


We want an approximate transfer function G, G that 
x(t) 


relates an original continuous-time signal X(t) to 
x(t) its ZOH reconstructed signal x(t) such that: 
' xe Xs) 


FIGURE 1.33 Approximate transfer function to model the hold delay. 


So, we want an approximate transfer function, model this hold 
delay which we'll call Gzo4(s, that will relate an original 
continuous-time signal x(f) to its ZOH reconstructed signal x,(s) 
as shown in Figure 1.33. So, let's briefly derive Gz,,,(s). If we pass a 
continuous-time signal x(t) through a sampler, we get the follow- 
ing sample signal x"(t), which we can mathematically represent 
with the following equation: 


v (t)=x(t) (tk) (1.42) 


If we take a Laplace transform of x*(t), we get the following 
expression: 


cx" (t) -x9- Ea) (1.43) 


Now driving this specific equation is beyond the scope of this 
chapter. So just take our word for it that this is the Laplace trans- 
form of x" (ft) as shown in Figure 1.34. 

So, if we were to pass x" (t) through the ZOH block, we get x;(t). 
However, x,(t) can be mathematically represented as a convolution 
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x(t) x*(t) 


VA 


t Sampler t 


FIGURE 1.34 Brief derivation of Gz5,,(s) - Part 1. 


between x*(f) and a function we'll call g(t) where g(t) is the follow- 
ing rectangular function: 


g(t)-u(t)-u(t-T) (1.44) 


If we took the Laplace transform of g(t), we get that G(s) equa- 
tion as below: 


£lg(t)j -G(s)- E (1.45) 


xD) xD 


ttt qe > TT 
LI t 7 = t 


In the time-domain, X , (1) is obtained by performing convolution between X * (t) and 8 ()., where 


g(t) =u(t) -u(t— T) 


l e 
L{g(t)}=G(s) = —- 
M 
&(t) Convolution 


x*(t) @g(th=x,(t) 


T 


FIGURE 1.35 Brief derivation of Gzo,(s) - Part 2. 
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Now we know from our linear signals and systems funda- 
mentals that the convolution in the time domain translates to 
multiplication in the frequency domain so we can express 
x (t) convolves with g(t) as X'(s) multiplied by G(s) which gets 
us x,(s), which is mathematically expressed as the following 
equations: 


x()es()ex()—9x()e()-X() ^ a9 


1v .2zk 
Dy PE G(s) =Xx,(s) (1.47) 


Now of this expression, if we take a closer look at the periodic 
train function, we only really care about the band-limited compo- 
nent which in itself is simply X(s) or the Laplace transform of our 
original continuous-time signal x(t). Therefore, realistically we 
only care about band-limited component (k = 0), which is simply 
X(s). So, we can equate the periodic train function to just X(s) as 
illustrated in Figure 1.36. 

Thus, our expression has been simplified to the following 
equation: 


1 
PX (5) G(s)=%.(s) (1.48) 
However, realistically we ONLY care about 
N " bandlimited component (k = 0), which is 
ON '" Ah simply X(s)... 
s= jw* à *s-jo 
SB jgk 


FIGURE 1.36 Brief derivation of Gzoy(s) - Part 3. 
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We know that G(s) equals to Equation (1.45). So, substituting 
this we get the following transfer function that relates the stair- 
case function x,(f) to our original continuous-time signal x(t): 


y X.(s) 3 ise 


(1.39),(1.36) x()^ sr 


(1.49) 


So, Gzoy(s) is approximately equal to the following equation: 


T e 


ST (1.50) 


Grou (s) = 


So, this is the transfer function block that models the hold delay 
in the continuous-time domain as shown in Figure 1.37. 

Now, we'll take this block and we'll place it in the closed-loop, 
and now we can finally go ahead and design C(s) while taking into 
account the continuous-time effects of the discrete parts. Therefore, 
we can now design C(s) using classical design techniques because 
we accounted for the discrete effects in our control loop as depicted 
in Figure 1.38. The closed-loop transfer function now looks like 
the following equation: 


MUN c(s)P(s) (1.51) 
R(s) 1«C(s)P(s)Gzou(s)H(s) 


x(t) 
^ ~ 
t UA 
/ qe > — NR 


Sampler On 


X(s) X.(s) 


FIGURE 1.37 Modeling the hold delay in the continuous-time domain. 
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FIGURE 1.38 Adding the hold delay block in the closed-loop. 


Sensor 


Now, after designing the continuous-time controller C(s), we 
want to discretize it using Tustin's method (bilinear transforma- 
tion). So, knowing that the sampling time is T and knowing that 
we're going to use Tustin's method we get the following z-domain 
transfer function: 


s= 
C(s) T z+1 =C(z) (1.52) 
T =chosen sampling time 


(1.53) 


C(z) U(z) by zb +z°b, +... +z “by 
 E(z) ay +z 'a+z7ayt...+2 “ay 


From our transfer function equation, (1.53), we cross multiply 
the terms as below: 


U(z ( Ay tz a, +z a, +...42 May) 
= E(z)(by +2 "b, +z°b, e. ez" by) (1.54) 


Next, we take the inverse Z-transform aka going to discrete- 
time domain, doing so; we end up with a different equation of our 
controller that relates our error input to the controller output as 
the following expression: 
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ayu| n ]- aiu| n1] au|n-2 ]...- auu[n- M 
= be|n *be[n-1]*be[n-2]«...« bue[ n- M | (1.55) 


Once we've got that we can go straight to implementing it on 
the Arduino, and this is the code for the feedback control struc- 
ture that we discussed in the previous section, and that is it. Those 
are all the preliminary steps needed to design and implement a 
simple digital controller for your system. In the last few sections, 
will apply these six steps to design a digital controller for a DC 
motor speed control system. 


1.8 DC MOTOR SPEED CONTROL 


1.8.1 System Modeling 


In this section, were going to dive into designing a digital feed- 
back control system to regulate the speed of a real DC gearmotor. 
The block diagram of the control loop is shown in Figure 1.39. 

In this case, we will have the PWM motor driver block after the 
ZOH block and the plant P(s) is the DC gearmotor, H(s) is the 
encoder, and output is the velocity (6(t)) in rad/s. We want to 
implement the rotary encoder that is used to acquire the motor 
speed which then gets sampled by the Arduino as a discrete signal 
which we'll call output speed. We then compare it with a desired 
speed variable to generate an error signal, e. This error signal then 
goes through our controller C(z) which again is implemented as a 


Arduino 
Vit) 
qe | : P(s) > O(t) 


Desired | 
speed (rad/s) 


Controller ZOH PWM motor 
driver 


KE 
f H(s) k 
Output_speed[n] Sampler 


Encoder 


DC Gearmotor 


FIGURE 1.39 The block diagram of the control loop. 
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FIGURE 1.40 The schematic of the setup. 


difference equation to generate a control signal. This control signal 
is realized as a PWM signal that comes out of the Arduino and 
goes into the PWM motor driver. The motor driver imposes an 
input voltage on the gearmotor to bring the actual speed, e(t), 
toward desired speed and thus we have our feedback loop. Here in 
Figure 1.40 is the apparatus for the DC motor control system. We 
have our Arduino Mega 2560 which is acting as our digital con- 
troller followed by a PWM motor driver which drives our gearmo- 
tor. Our motor conveniently has a built-in encoder that sends 
pulses data to where do we know which contains an algorithm to 
convert this data to an output speed value. The motor driver 
receives this power from the DC power supply and the real output 
voltage is set to 12 volts. The schematic of the setup is depicted in 
Figure 1.40. 

So, we want to design a digital controller to effectively regulate 
the speed of this DC gearmotor. So, the following are the six-step 
procedures as discussed in the previous section. 
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1. Obtain the plant’s transfer function P(s) 
2. Choose an approximate sampling time T 


3. Design controller C(s) while taking into account the effects 
of ZOH/sampling 


4. Use bilinear transform/Tustin’s method to discretize C(s) to 
C(z) 


5. Take the inverse Z-transform of C(z) to obtain the difference 
equation 


6. Implement the difference equation on the Arduino and run 
the system! 


Our first task is to obtain the approximate transfer function of our 
DC gearmotor. So, in order to do that we must drive a transfer 
function based on the electro-mechanical system model of a 
generic DC motor as illustrated in Figure 1.41. 

Therefore, in order to design a speed control system for this 
gearmotor using classical control techniques, we must first derive 
its transfer function, P(s), using an equivalent system model that 
approximately models its dynamics. Now running out all the 
physics for the mathematical model derivation is beyond the scope 
of this chapter. So, we'll go straight into the transfer function of 


X Fixed 
^ field 
"N R L 
R "o 4 a {i 


Armature 
circuit 


Electro-mechanical System Model 
34:1 Metal DC Gearmotor 
Nominal: 12VDC (220 RPM) 


FIGURE 1.41 The DC motor transfer function - Part 1. 
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this model that relates the input voltage, V(s), to the output speed, 
0 (s) , as the following equation [4]: 


(s) se(s)- K rad/s 
P(s) v(s) V(s) (je B | (1.56) 


where J (kg.m?) is the moment of inertia of the rotor, b (N.m.s) is 
the motor viscous friction constant, K is equal to the electromo- 
tive force constant (K,(V/rad/s)) or motor torque constant (K,(N.m/ 
Amp)), R(ohm) is the electric resistance, and L(H) is the electric 
inductance. From Figure 1.41, we can derive the following govern- 
ing equations based on Newton's second law and Kirchhoff’s volt- 
age law: 


J6+b0 - Ki - Ki (1.57) 
di ; 
L' *Ri-V-K-V-KÓ (1.58) 


By applying the Laplace transform to Equations (1.57) and 
(1.58), we will have the following equations: 


(Js b)s0(s) = Kr (s) (1.59) 
(Ls+R)I(s)=V(s)—Ks0(s) (1.60) 


By eliminating I(s) from Equations (1.59) and (1.60), we get the 
open-loop transfer function in Equation (1.56). Now, if you wish 
to see where we got this transfer function from, you can find the 
detailed derivation in the link provided in [4]. So, the five param- 
eters of the DC motor transfer function are as below as shown in 
Figure 1.42. 
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Ü b = motor's viscous friction constant (N.m.s) 
K = electromotive force constant (V/rad/s) 

R = electric resistance (Q) 


L = electric inductance (H) 


O(s) K 


V(s) (Js+b) (Ls R) - K2 


FIGURE 1.42 The DC motor transfer function - Part 2. 


Jis the rotors moment of inertia (kg.m?) 

b is the motor's viscous friction constant (N.m.s) 

Kis the electromotive force or motor torque constant (V/rad/s) 
Ris the electric resistance (Q) 


L is the electric inductance (H) 


So, now the big question is how do we find these values for our 
gearmotor? Well, if you're lucky all this information can be avail- 
able on the motor's data sheet. However, almost 9996 of the time 
this is never the case, so, we'll have to somehow figure out these 
values using our tools and our intuition. So, let's try to figure out 
K, now K is the electric motor force or torque constant. This is 
essentially the DC gain that relates our input voltage to our output 
speed in radians per second (rad/s). Now, according to the data 
sheet, the nominal speed at the rated voltage of 12 volts is 220 
RPM now, 220 RPM can also be expressed as 23 rad/s since we 
have: 


2m (rad) 
220 RPM = 220 x ——— = 23 rad/s (1.61) 


60(s) 
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So, dividing 12 V by 23 rad/s, we get an approximate constant of 
about 0.52 V/rad/s. Now what about the resistance? Now luckily, we 
have an LCR meter which can measure electrical resistance of 
complex impedance. So, if we probe across the leads of the motor, 
we can get an approximate value of the internal resistance. Now 
the resistance and the inductance actually vary with the angle of 
the rotor. But again, were just looking for an approximate value. 
So, for now, we'll just say that the resistance is around 4.4 ohms. 
Now using the same meter which is capable of measuring induc- 
tance values at several frequencies we measure the inductance to be 
around 6.3 mH, at an applied low frequency of 100 Hertz. So far, 
we've gotten three values down, great, but now how do we go about 
finding / and b. Now, unless we have some very expensive precision 
tests' apparatus, finding even the approximate values for these 
parameters is quite difficult. Luckily, MATLAB Simulink has a 
parameter estimation feature in the optimization toolbox. The 
design optimization toolbox comes with a parameter estimation 
tool, where the user can estimate the parameters of a simulated 
model of their system using real measured input and output data. 

So, the first step to use a tool is to create a Simulink model 
based on the differential equations that describe the dynamics of 
the system or to create a linear model that captures the dynamics 
of a DC motor (Equation 1.56 or Equations 1.59 and 1.60) as 
shown in Figure 1.43. 


i L 
JZ +b0=Ki EE >I 


L9. +Ri=V-Kò 
dt 


16 


FIGURE 1.43 The Simulink model based on dynamic equations. 
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Now, if you've never played with Simulink or modeled on 
Simulink, we highly suggest that you get exposed to it because it’s 
very useful for control design and simulating real-time systems. 
There are plenty of tutorials all over YouTube and Udemy on mod- 
eling using MATLAB Simulink, so we won't be discussing that 
here as it’s beyond what this chapter focuses on. The second step is 
to obtain the input and output data we want to use. So, in this 
Arduino script, we've created a test input signal that will vary the 
input voltage applied to the motor in a step-like fashion and will 
want to monitor the elapsed time, the applied input and the out- 
put speed. 


int sampling time - 50; 

unsigned long previous time, current time; 
float encoder count - 0; 

boolean encoder A, encoder B; 

byte state, statep; 

float previous angle, current angle - 0; 
float motor speed = 0; 

float supply voltage 
float elapsed time = 0; 
const int pwm output - 7; 
float control; 

void Encoder State() { 
encoder A - digitalRead(3); 
encoder B - digitalRead(2); 


Il 
H 
N 


if encoder A == HIGH) && (encoder B == HIGH) ) 
state - 1; 

if encoder A -- HIGH) && (encoder B -- LOW)) 
state - 2; 

if encoder A -- LOW) && (encoder B -- LOW)) 
State - 3; 

VE encoder A -- LOW) && (encoder B -- HIGH)) 
state - 4; 

switch (state) 

{ 


case 1: 
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{ 


if (statep == 2) encoder count--; 
if (statep == 4) encoder count-«-*; 
break; 
case 2: 
if (statep == 1) encoder count-«-; 
if (statep -- 3) encoder count--; 
break; 
case 3: 
if (statep == 2) encoder count-«-*; 
if (statep -- 4) encoder count--; 
break; 
} 
default: 
if (statep == 1) encoder count--; 
if (statep == 3) encoder _count++; 
break; 


} 
} 
statep = state; 
} 
void get_speed() { 
current angle = (encoder_count*360.0) / (1632.67) ; 
motor speed = 0.01745*((current angle 
- previous angle)/0.05); 
previous angle - current angle; 
} 
void step input () { 
control = 0; 
if (current time >= 2000 && current time <= 
6000) { 
control = 12; 


} 
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if (current time > 6000 && current time <= 
10000) { 

control = 6; 
if (current time > 10000 && current time <= 
14000) { 

control = 8; 
if (current time > 14000) { 

control = 0; 


void setup () { 

Serial.begin(9600) ; 

previous time = millis(); 
pinMode (2, INPUT) ; 

pinMode (3, INPUT) ; 

attachInterrupt (digitalPinToInterrupt (2), 
Encoder State, CHANGE); 

attachInterrupt (digitalPinToInterrupt (3), 
Encoder State, CHANGE); 

encoder A - digitalRead(3); 

encoder B - digitalRead(2); 


LE encoder A == HIGH) && (encoder B == HIGH) ) 
statep = 1; 

LE encoder A == HIGH) && (encoder B == LOW) ) 
statep = 2; 

if encoder A == LOW) && (encoder B == LOW) ) 
statep = 3; 

LE encoder A == LOW) && (encoder B == HIGH) ) 
statep = 4; 

pinMode (pwm output, OUTPUT); 


} 


void loop() { 
current time = millis(); 
int delta time - current time - previous time; 
if (delta time > sampling time)[( 
elapsed time - current time/1000.0; 
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step input(); 

get speed(); 
analogWrite(pwm output, (control/ 
supply voltage)*255.0); 


Serial.print (control); 
Serial.print(" Wy ss 
Serial.print(elapsed time); 
Serial.print(" "); 
Serial.println(motor speed); 
previous time - current time; 
) 

) 


Let's run the script and then from Tools menu, we select the 
Serial Monitor option. Now we'll take 15 seconds rows of data and 
then on MATLAB, we'll create 7 data variables that contain R, L, 
K, J, b, and also, the measured Input Data and Output Data 
which the first column is the elapsed. time and the second col- 
umn is control and motor speed, respectively, as depicted in 


Figure 1.44. 


Now, before the estimation tool can run, we'll need to use some 
initial values for all the parameters. So earlier we approximated 
what K, R, and L could be. So, that's what we'll be using for their 


initial values as illustrated in Figure 1.45. 


| R 


Za Variables - Output Data 
; NERIK 
294x2 double 


FIGURE 1.44 Creating variables that contain the measured data. 
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FIGURE 1.45 Obtaining parameter values. 


However, for b and J, well just pick some random values to 
begin with (J = 0.001, b = 0.1). Now, we are opening up the 
Simulink, and the file containing our motor model. 

We'll go to the apps tab and then open up parameter estimator. 
From this window, we want to create a New Experiment. Next, we 
want to pull the measured Output_Data and Input_Data from the 
workspace and then we want to select the parameters (Select 
Parameters), since we wish to estimate all those five parameters 
(J, K, L, R, and D). Therefore, we select all of them and click on the 
OK button. 

Now, in order to get more realistic parameter estimation, we 
want to set the Min and Max boundaries so that the tool doesn't 
vary any parameters too far from a realistic value. So, we'll set the 
Min and Max for K, R, and L to what we think they could be, and 
because we don't know b and J, we'll leave its max alone but we'll 
want to set the minimum value to zero because realistically none 
of these parameters should ever go negative as below. 


J= 0.001, Minimum = 0, Maximum = Inf 

K = 0.52, Minimum = 0.47, Maximum = 0.54 

L = 0.0063, Minimum = 0.006, Maximum = 0.007 
R= 4.4, Minimum = 3.8, Maximum = 5 


b = 0.1, Minimum = 0, Maximum = Inf 
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Now, if you click on Plot button in that window, you should see the 
measured input voltage and output speed data we got from the 
Arduino. 

Now let's click on the Estimate button, you can see that the tool 
is running some sort of optimization algorithm to curve fit the 
simulation output to the measured output and it's doing this by 
iterating and varying the parameters. Now, once we've come to a 
close enough fit, the tool stops estimating and our parameter val- 
ues should have automatically been updated by clicking on the 
Exp icon in Experiments sub-window. 

Now, using the following script as given below, we can get the 
approximate transfer function for our DC gearmotor as the fol- 
lowing equation: 


(1.62) 


2 (s) _ 0.47 
3.425e — 5s? +0.01859s +0.2211 


J 0.004893; 

b 5.8147e-5; 

K = 0.47; 

R 3.8; 

L 0.006999; 

num = K; 

den = [J*L, (J*R)+(b*L), (b*R) +(K*2)]; 
sys = tf(num,den); 


So, step 1, which is obtaining the plant's transfer function P(s), is 
complete. So, let's stop here for now, until now, we've gotten an 
approximate transfer function. In the next section, we'll move on to 
selecting an appropriate sampling time and design your controller. 


1.8.2 Controller Design 


In this section, we will continue the controller design process for 
DC motor speed control example. So, we dedicated the previous 
section to obtaining an approximate transfer function for DC 
gearmotor which is shown in Figure 1.46. 
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FIGURE 1.46 Obtaining an approximate transfer function for DC 
gearmotor. 


So, moving on to step 2, we must find a sampling time, T, for 
our control loop. So, as we discussed a couple of sections ago, the 
sampling frequency should be at least 5-10 times higher than the 
bandwidth of the open-loop system, P(s). Using the following 
MATLAB code, we can obtain the Bode plot for our system. 


s = tf('s'); 

P s= 0.47/(3.425e-5*s^2 + 0.01859*s + 0.2211); 
bode (P_s) 

grid on; 


1 


From observing the magnitude plot, we can see that the DC 
gain is around 6.55 dB and the bandwidth of the system is at the 
frequency 3 dB below the DC gain, so, moving down 3 dB below, 
we get to 3.55 dB and the frequency at this gain 12.1 rad/s or 1.93 Hz. 
So, the bandwidth of our system is 1.93 Hz or around 2 Hz. So, 
we'll multiply this by 10 and arrive at a sampling frequency of 20 
Hz which correlates to a sampling time of 50 milliseconds as 
shown in Figure 1.47 and equation below: 


Tsis sübscands (1.63) 
f 20 


20 


Bode Diagram 


"m 
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FIGURE 1.47 Choosing a sampling time, T. 
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FIGURE 1.48 The continuous-time block diagram for our control loop. 


So now that, we have our sampling time and step 2 is complete. 
Now let’s move on to designing a continuous-time controller C(s), 
as shown in Figure 1.48, which is the continuous-time block dia- 
gram for our control loop and here after C(s) we have a PWM driver 
block with gain of 1, the P(s) is the transfer function of DC gearmo- 
tor obtained in the previous section, and H(s) is the transfer func- 
tion of encoder with the gain of 1. Because the encoder is simply 
providing pulse count data that is to be interpreted by the Arduino 
as a speed value in rad/s, we can just set the encoders transfer func- 
tion, H(s), as a gain of 1, so we have the following equation: 


H(s)=1 (1.64) 


So, this is our closed-loop transfer function that relates a set 
point to our output speed as the following equation: 


P (s) n C(s)P(s) 
Garl ) R(s) 1+C(s)P(s)Gzon (s)H(s) 


(1.65) 


Also, knowing that each H(s) is equal to 1, Equation (1.54), we 
get the following expression: 


(s C(s)P(s C(s)P(s 
Cul) FOG a I2 


T(s)- C(s)P(s)Gzou(s)H(s) (1.67) 
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FIGURE 1.49 Adding explicit transfer functions for P(s) and Gzoy(s). 


With T(s), is being our loop gain. Here as depicted in Figure 
1.49 is a closed-loop block diagram except with the explicit trans- 
fer functions shown for P(s) and Gyo,(s). 

Therefore, the explicit transfer functions for Gzoy(s) are given 
in the following equation: 

1-e^77T 


Grou (s) x. , T=0.05sec (1.68) 


Now, its quite apparent that designing a controller from this 
point can be quite difficult because G;54,(s) has an exponential 
term and therefore it's irrational. We should approximate Gzo,(s) 
further to a rational function to make the control design easier for 
MATLAB. The control design tools in MATLAB tend to work 
more nicely with the rational transfer functions as opposed to 
transfer functions that have exponential terms in them. So, using 
an approximation technique known as the Pade approximation, 
the exponential term can be approximated as a rational function 
with ever-increasing orders of precision, so the higher the order, 
the better the approximation as the equation below: 


l x4 : x? J x^4 1 x... 
9 72 1008 30240 
l x? l x 1 x* ! x. 
9 72 1008 30240 


(1.69) 


Pade Appr.:e* = 
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Now, of course, doing this by hand is impractical but luckily 
MATLAB has a function (available with Control System Toolbox) 
called pade(sys,n) with sys = transfer function, n = approximation 
order, which if you give it the irrational transfer function and a desired 
order n, it will output the approximate rational transfer function. So, 
using the following MATLAB script, we can obtain the following 
rational transfer function for G,o,,(s) and we'll call that Gzopp(9): 


40s? +960000s 
s^ +240s* +24000s* +960000s (1.70) 


Gzoup (s) = 


T = 0.05; $sampling time T = 0.05 seconds 

s = t£('s'); 

G zoh = (1-exp(-s*T))/(s*T); %original ZOH/ 
sampling delay transfer function 

G zohp = pade(G zoh,3); % 3'?» order approximation 


So now that, we know that both P(s) and Gzopp(s) are rational 
and we can go ahead and design C(s) as shown in Figure 1.50. 
Now, for controlling the speed and position of a DC motor, 
typically a proportional integral (PI) controller is used in the 
industry and most conventional DC motors have their output 
speed regulated using a PI controller as the following equation: 
K; 


C(s)=K,+— (1.71) 
S 


+ 0.47 : 
R(s) = C(s) 3.425e-05s? + 0.018595 + 0.2211 O(s) 


Controller DC Gearmotor (P(s)) 


40s? + 960000s 
s! + 2405? + 240005? + 960000s 
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FIGURE 1.50 The rational transfer functions for P(s) and Gzopp($) in 
the control loop. 
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So now, we can design C(s) using classical control techniques such 
as root locus, frequency response analysis, or some other complicated 
control design method but look at the transfer function we're dealing 
with in total they make up a six-order system. Now if you want to do 
root locus or frequency response analysis on a six-order system, you 
can. However, we think it's a lot more convenient to use the control 
design tools that are available with MATLAB. Now luckily because 
we're designing a controller that's a part of the proportional integral 
derivative (PID) family, we can use the pid Tuner(sys) app that's avail- 
able with the control system toolbox and the code is as below. 


Ty 2:005; $sampling time T - 0.05 seconds 

Ss: erwi 

P S= 0.47/(3.425e-5*s^2 + 0.01859*s + 0.2211); 
$DC gear motor’s transfer $function P(s) 

G zoh s = (1-exp(-s*T))/(s*T); %ZOH/Sampling’s 
continuous-time effect transfer function G zoh(s) 
G zohp s = pade(G zoh s,3); %3% order Pade 
approximation of G zoh(s) 

sys = P s*G zohp s; %total open-loop gain 
pidTuner(sys); %launch PID Tuner App on MATLAB 


We declare our transfer functions P(s) and Gyo p(s) and then 
we multiply them to get the total effect of open-loop transfer func- 
tion. Then we store it in a variable which we'll call ‘sys’ which is 
short for system. We then call the pidTuner(sys) function with 
'sys being the input argument. Now, before we run the script, we 
want to set some design requirements. We would like the settling 
time for the output transient to be less than 2 seconds. The over- 
shoot and undershoot to be under 596 of the steady-state value and 
the steady-state error to be less than 396 of the steady-state value. 
So, the design requirements are: 


e Settling time less than 2 seconds 
e Overshoots and undershoots should be under 596 


e Steady-state error should be less than 3% 
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So now let’s run the script, so the PID tuner raptured open up 
automatically and by default it seems like it chose a PI controller 
on its own, and what you see now is a step response of the closed- 
loop system’s output and upon the top that are sliders that we can 
use to adjust the response time and transition behavior such that 
it meets our criteria, and on the bottom right should be the K, and 
the K; gain values that it correlates to. So, we'll play around with 
this plot until we think it looks good. So, this response seems to 
work and if we click on the Show Parameters tab, we can see the 
gain values, the settling time, overshoot percentage, and the sta- 
bility margins as depicted in Figure 1.51. 

So far, they all meet our criteria and we think we'll go with this 
design, so, this is our continuous-time PI controller with K, that is 
equal to 0.427 and K; is equal to 5.114 and we have the following 
equation: 


5.114 
S 


(1.72) 
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FIGURE 1.51 Clicking on the Show Parameters tab. 
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FIGURE 1.52 Complete design of continuous-time control loop. 


So now that, we have got our continuous-time controller, and 
step 3 is complete as illustrated in Figure 1.52. 

Now, we have to discretize it using Tustin's method by running 
the following MATLAB script that contains a c2d function with 
the method set to Tustin. 


T 0.05; $sampling time T - 0.05 seconds 
S t£('s'); 

C s = 0.427 + (5.114/s); %continuous-time 
controller C(s) 

C z = c2d(C_s,T,’tustin’); tdiscretized 
controller C(z) 


We get the following discrete transfer function C(z), and we 
divide the numerator and the denominator by ‘z’ to get it into his 
causal expression as the equation below: 


c(z)- U(z) _ 0.5548z -0.2991 _ 0.5548-0,2991z" 


E(z) z-1 1-z1 om 


So now, step 4 is complete and now we have to take the inverse 
Z-transform to get the difference equation. Now we take C(z) 
from Equation (1.73) which relates our error E(z) to our control 
signal U(z) and cross multiply the terms and, after distributing, 
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we then take the inverse Z-transform to get a difference equation 
as the following equations: 


U(z)( 1-2) =E(z)(0.5548—0.29912") (1.74) 


U(z)-U(z)z ' = E(z)0.5548— E(z)0.2991z * (1.75) 


(1.64) — 5 u[n]-u[n-1] = 0.5548e[n ]--0.2991e[n -1] 
(1.76) 


Then we bring all the delayed input and output terms to the 
right side as the equation below: 


u(n |- u[n-1]«0.5548e[ n | -0.2991e| n-1] (1.77) 


Now, we've obtained the controller's difference equation that 
relates our discrete error signal e[n] to our discrete control output 
u[n]. So now, step 5 is complete. Now, all we have to do is to imple- 
ment the controller on the Arduino and see the closed-loop speed 
control in action and we'll do this in the next section. 


1.8.3 Arduino Demonstration 


In this final section, we'll finally implement the digital controller 
designed in the previous section to regulate the speed of our DC 
gearmotor. So, we dedicated the last two sections to completing 
steps 1-5. Now all that remains is implementing the controller on 
the Arduino to demonstrate the closed-loop control. However, 
before we go ahead and code up the script, we are going to first 
simulate this control structure to see what type of behavior we 
should expect. Shown here in Figure 1.53 is the Simulink block 
diagram for our closed-loop feedback system. Notice that there's a 
saturation block after the controller. Now what does this block 
will do is to limit the control signal to be within the bounds of 
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FIGURE 1.53 The Simulink block diagram for our closed-loop feedback 
system. 


0-12 volts and that’s because our real power supply for our PWM 
driver is set to 12 volts. 

Also, because the control signal of 12 V correlates to a PWM 
duty cycle of 100% so it’s inherently the maximum voltage we can 
apply. So, a good way to test the feedback system is by stepping the 
desired set point to a few speed values and observing how well the 
output speed tracks the set point. So, using the staircase signaling 
block, we create the set points signal where every four seconds, the 
set point changes abruptly to a few different values and we're 
observing the plots for the desired set point and the measured 
speed as shown in Figure 1.53. Let's run the file, while observing 
the plot, you can see that it tracks the set point speed very well and 
it also meets our design right here in regard to the settling time, 
overshoot, and steady-state error as shown in Figure 1.54. 

So, in simulation, our results look good. Now, let’s go imple- 
ment this exact same test on the Arduino setup. So, here is the 
Arduino script. 
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FIGURE 1.54 Plot of tracking the set point speed. 


int T = 50; // sampling time T - 50ms 
unsigned long previous time, current time; 

// variables to hold time values 

float encoder count - 0; // initialize encoder 
count to 0 

boolean encoder A, encoder B; 

byte state, statep; 

float previous angle, current angle - 0; 

// set current and previous angle //values to O0 
float output speed - 0; // initialize 
output speed to 0 

float supply voltage - 12; // the supply 
voltage for our PWM motor driver //is 12V, so we 
set this to 12 

float elapsed time = 0; 

float setpoint; // declare setpoint variable 
float e, e 1, u, u_l = 0; // delcare e[n], 
e[n-1], u[n], and u[n-1] 
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const int pwm output - 7; // label digital 
pin 7 as the PWM output 
void Encoder State()[ // this function is 
used for encoder pulse counting 
encoder A - digitalRead(3); // read pulse 
data from encoder's pin A 
encoder B - digitalRead(2); // read pulse 
data from encoder's pin B 
if encoder A -- HIGH) && (encoder B -- HIGH)) 
state - 1; 
if encoder A -- HIGH) && (encoder B -- LOW)) 
state = 2; 
if encoder A == LOW) && (encoder B == LOW) ) 
state = 3; 
if encoder A -- LOW) && (encoder B -- HIGH)) 
state - 4; 
switch (state) 
{ 
case 1: 
{ 
if (statep == 2) encoder count--; 
if (statep == 4) encoder count-«-*; 
break; 
} 
case 2: 
{ 
if (statep == 1) encoder _count++; 
if (statep == 3) encoder count--; 
break; 
} 
case 3: 
{ 
if (statep == 2) encoder _count++; 
if (statep == 4) encoder count--; 
break; 


} 


default: 


{ 
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if (statep == 1) encoder count--; 
if (statep == 3) encoder count--*; 
break; 


) 
) 
statep - state; 
) 
void get_speed() { 
current angle = (encoder _count*360.0) /(1632.67) ; 
// 1632.67 is the total //number of counts per 
revolution. Because it’s a 48 CPR encoder and a 
34.1 //gearmotor, you multiply 48 x 34.1 = 
1632.67 
output speed = 0.0175*((current angle - previous _ 
angle)/0.05); // 0.0175 rad/s = 1 degree. 
The formula is essentially calculating theta dot 
// or dteta/dT, where dT - 0.05, is the sampling 
time in seconds 
previous angle - current angle; 
) 
void step setpoint () { // a stair-case setpoint 
input signal to test the //transient response of 
the control loop 
setpoint = 0; 
if (current time >= 4000 && current time <= 


8000) { 
setpoint = 5; // after 4 seconds, 
setpoint = 5 rad/s for 4 seconds 


} 


if (current time > 8000 && current time <= 
12000) { 
setpoint = 20; // setpoint = 20 rad/s 
for 4 seconds 
} 
if (current time > 12000 && current time <= 
16000) { 
setpoint = 12; // setpoint = 12 rad/s 
for 4 seconds 
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if (current time > 16000) { 
setpoint = 0; 
} 


void setup () { 
Serial .begin(9600) ; 
previous time = millis(); 


pinMode (2,INPUT); // set digital pin 2 
as an input 
pinMode (3,INPUT); // set digital pin 3 


as an input 
attachInterrupt (digitalPinToInterrupt (2), 
Encoder State, CHANGE); 

attachInterrupt (digitalPinToInterrupt (3), 
Encoder State, CHANGE); 

encoder A - digitalRead(3); 

encoder B - digitalRead(2); 


if encoder A -- HIGH) && (encoder B -- HIGH)) 
statep - 1; 

if encoder A -- HIGH) && (encoder B -- LOW)) 
statep - 2; 

if encoder_A == LOW) && (encoder_B == LOW)) 
statep = 3; 

LE encoder A == LOW) && (encoder B == HIGH) ) 
statep = 4; 

pinMode (pwm_output, OUTPUT); // set pwm_ 


output or pin 7 as the PWM signal’s output pin 
} 
void loop() { 

current time = millis(); // store the elapsed 
time in milliseconds 

int delta time = current time - previous time; 
// calculate the delta //time or time difference 
in milliseconds 

if (delta time > T) { // if the time 
difference is equal to 50 ms, execute the if 
Statement. Doing so, keeps a consistent sampling 
interval of //50ms 
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elapsed time = current time/1000.0; // stores 
the elapsed time in seconds 


step setpoint(); // call out desired 
setpoint speed signal 

get speed(); // obtain our output 
Speed signal in rad/s 

e - setpoint - output speed; // calculate 
the error e[n] 

u =u 1 + 0.555*e - 0.299*e 1; // calculate 


the control signal u[n], this is our difference 
equation derived from C(z) 
if (u >= 12.0) { // if statement is used to 
saturate the control signal value between 0(V) 
and 12 (V) 

u = 123 


else if(u <= 0) { 

Ui = 05 

} 
analogWrite (pwm output, (u/supply_voltage) *255.0); 
// maps the control signal values linearly to 
PWM duty cycle %’s. 0 - 12 correlates to 0% 
- 100%. 


u_l = u; / store u[n] --» u{n-1] for next 
iteration 
el=e; // store e[n] --> e{n-1] for next 
iteration 
Serial.print (elapsed_time) ; // Serial print 


the elapsed time(s) and the output_speed or 
measured speed (rad/s) 

Serial.print(" ") ; 
Serial.println(output speed); 


previous time - current time; // store 
current time --» previous time for next loop 
iteration 


} 
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You can see that the controller is implemented as a difference 
equation as we expected inside the void loop function. You can 
also see that we've declared an input set point function (step. set- 
point) that mimics the input set point signal from the simulation. 
We've also included some code to saturate or limit the control 
output between 0 and 12. So, let's monitor the set point signal and 
the measured output speed on the Serial Plotter and observe the 
DC gearmotor's behavior as depicted in Figure 1.55. 

So, it seems to be working as expected. You'll notice that there's 
a lot of noise in the measured output. This is mainly due to the 
noise in the encoder speed measurement, but again it's very mini- 
mal and the overall transient and steady-state performance is 
pretty good. With built-in encoders, we dont really expect them 
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FIGURE 1.55 Observe the DC gearmotors behavior on the Serial 
Plotter. 
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to have very precise or noiseless performance, but overall, this 
control structure works and does exactly what we expected it to. 
So, we did the MATLAB Simulink simulation, and then we did 
the Arduino implementation on the real DC gearmotor. So now, 
let's compare the two results with one another and see how closely 
they match or mismatch. We'll modify our Arduino script to out- 
put the elapsed time and the output speed and collect 17 seconds 
with the data using the Serial Monitor. 

Then on MATLAB, we'll create a variable called Measured . 
Speed and enter in two columns of that, this elapsed time (column 
1) and measured output speed (column 2) data as shown in 
Figure 1.56. 


Wá Variables - Measured S 
Measured Speed >% 
HH 333x2 double 


OO YAU RU MN 
o 
: 


o 

UJ | 

z 
lolololololololololololo 


FIGURE 1.56 Creating a variable string with this measured output 
speed data. 
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From Workspace 


Read data values specified in timeseries, matrix, or structure format from the 
MATLAB workspace, model workspace, or mask workspace. 


MATLAB timeseries format may be used for any data type, complexity, or fixe: 
Waveform dimensions. To load data for a bus signal, use a MATLAB structure that 
Siondcalor matches the bus hierarchy and specify timeseries for each leaf signal. 


Speed(rad/s) 


For matrix formats, each row of the matrix has a time stamp in the first 
column and a vector containing the corresponding data sample in the 
subsequent column(s). 


For structure format, use the following kind of structure: 


var.time=[TimeValues] 
var.signals. values=[DataValues] 


—€— var.signals.dimensions- [DimValues] 
Workspace 


Parameters 
Data: 
Measured Speed 


Output data type: [ Inherit: auto 


Sample time (-1 for inherited): 
[0.05 
< 


9 


FIGURE 1.57 Feeding data from the workspace and plot against the 
simulated data - Part 1. 


Then on Simulink, we'll feed this data to From the Workspace 
block and fill in its Data as Measured Speed and set the Sample 
time as 0.05 and plot on the scope against the simulated output 
data and the simulated set point as shown in Figures 1.57 and 
1.58, respectively. 

As shown in Figure 1.59, the theoretical result and measured 
result match very well with the exception of the noise on the mea- 
sure data and you can see that even the transient parts of the mea- 
sure data match pretty close with the simulated data. 

So, this comparison proves three things. First, the transfer 
function we approximated for our actual DC gearmotor was a 
good enough approximation to use for our controller design step 


* 0.55482-0.2991 
z-1 


Waveform 
generator 


Discrete 
transfer Fcn 


Measured Speed 


From 
Workspace 


H- 


Saturation 


Zero-Order 
hold 


0.47 


3.425e-5s2+0.01859s+0.2211 


Transfer Fen 


FIGURE 1.58 Feeding data from the workspace and plot against the simulated data - Part 2. 
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Ready Sampie based [T=17.000 


FIGURE 1.59 Comparison of the theoretical result and measured result. 


process. Second, the closed-loop control structure on Simulink is 
a very accurate simulation of the real physical closed-loop system, 
and, third, the six-step controller design process is theoretically 
and practically a sound method to use to design a simple digital 
controller for your continuous-time system, so that is all. 


CHAPTER 2 


Theoretical Overview 
of the Buck Converter 


2.1 BUCK CONVERTER DESIGN 


In this chapter, we will cover the basic design of a buck converter 
and implement closed-loop control for it. This system has both 
voltage and current sensors. We do not dive deeply into theoretical 
analysis of buck converters, because this book is practically ori- 
ented. So, we will not consider the derivation of the buck converter 
equations and we will just present them. The schematic of the basic 
buck converter is shown in Figure 2.1. 

We have an input voltage, Vin, a MOSFET, Q, a diode, D, an 
inductor, L, a capacitor, C, and a resistive load. The load can be 
active, for example, a battery or other type of loads, but in this 
case, a resistor has been used for analysis. We have also the output 
current, L,, inductor current, I, and output voltage, V,. The 
parameters of the buck converter which are constant and designed 
for a particular case are listed in Table 2.1. 
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FIGURE 2.1 The schematic of basic buck converter. 


TABLE 2.1 Parameters of the Buck Converter 


Parameters 
Vin 


V, 


out 


La = ij) 


Values 


24V 
16V 
2.133 A 
7.5 Q 
0.02 V... 
0.1 Ion 

3 mQ 
300 mQ 
1.15 V 
0.075 V 
200 uH 
40 uF 
100 kHz 


For the sake of analysis, we choose the inductor current ripple to be 
1096 of inductor current, using the following equations for duty cycle, 


D, and inductor, L, we get the duty cycle and inductance values: 


C Vwa + m Ty Vout 


Vin T Viwa P. Vies 


= 0.7094 


(2.1) 
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> (Vin = V cessat A irr, a Vout )x D 


= 243.3 uH (2.2) 
AT pk-pk fow 


Because of the size of the inductor, we choose the inductor L to 
be 200 uH and let the inductor current ripple to be 15% of induc- 
tor current. Additionally, the output voltage ripple, which is 2% of 
the output voltage, and the actual capacitor value (40 uF) signifi- 
cantly exceed the calculated value for capacitor C (0.83 uF), as 
shown in the following equation: 


> AT, pk-pk 


> PE 0.83 uF (2.3) 
SAV: accu 


Then we have the following small signal transfer function 
which is crucial in designing the buck converter [5]: 


to (s) a (1+ sCrg )(Vin + Visa —Vas sat ) 


d(s) (LCR + LCrggn )s? + (CRr, + Cresrty + CRrigsp + L) SF TL + Rout 
(2.4) 


We should note that these equations cannot be used to other 
converters. Also, the output voltage to duty cycle small signal 
transfer function is as below: 


Vo (s) Rout (1 T SCTesr \(Vin + Visa T Vee ) 


d(s) (LCR+LCresr )s* +(CRr, + Cresan + CRresn +L) s ri + Rout 
(2.5) 


Notice that the parasitic parameters, that is, equivalent series 
resistance of capacitor and inductor resistance, are denoted by ry, 


84 m Digital Control of Power Converters 


and r, respectively. If we substitute the parameters in Table 2.1 in 
the above equations, we will have the following relations: 


io(s) _ 3.009 10 55-- 25.07 (2.6) 
d(s) 6.003 x10? s? +0.0002912s +7.801 


vo(s) - 2.257 x10°s +188.1 (2.7) 
d(s) 6.003 x10? s? +0.000291s + 7.801 


For designing a PI controller for current control loop or voltage 
control loop, we can use MATLAB. In the previous chapter, we 
used MATLAB functions (pidTuner and c2d tools) for designing a 
digital PID controller for our control loop. To design a PI control- 
ler for the current controlloop using MATLAB, the following code 
can be applied: 


S z-tfí'st); 
Gid = (3e-6*s425.07)/(6e-8*8^240.000291*$47.801); 
pidTuner (Gid); 


By executing the above code, the PID tuner window is popped up. 
We select the type of controller as PI and its form as Standard and 
then regulate the Response Time and Transient Behavior sliders to 
get desired Step Plot. Then, if we click on Show Parameters button, we 
can see our controller parameters. Finally, we click on the Export 
button and give it a name for our PI controller. In this case, we chose 
‘Cont’ as its name and click on the OK button as shown in Figure 2.2. 

Now, in the command window, if we type the name (‘Cont’) 
and then press the Enter key, you can get the controller transfer 
function as the following equation: 


E 
T, i 
G(s) [1e e saoe +202 (2.8) 


is 
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22 dB @ 1.14e-04 rad/s 
(89.2 deg © 383 rad/s 
[stable 


J 
0.02 0025 


Controller Parameters: Kp = 3.108e-05, Ti = 2.578e-07 


FIGURE 2.2 The PID tuner window. 


K, =3.11x10”,T, =2.58x107 (2.9) 


If you want to implement analog controller using op-amp, the 
above equation is enough. However, for use in a digital controller, 
the following substitutions must be applied (Tustin’s method): 


_2(z-1)_2(1-2") 
m (z+1) T. (127) us 


where T, is the sampling time and is chosen as 100 us, so we have: 


T,(z+1 
G.(z) = orn] 
K,(z-1)* K,T, (z 1)/ (27;) 
(z-1) (2.11) 
(K, - K,T. /(2T;))z+K,pT, /(21;)-K, 


z-l 


0.0061z 4- 0.006 
z—1 
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0.006z ! -- 0.0061 
G.(z)= cem (2.12) 


We can also use the c2d function in MATLAB to get the 
discrete-time transfer function from G,(s) using the following 
code in the command window as depicted in Figure 2.3. 


c2d(Cont,100e-6,'tustin') 


Correlating the equation above with the standard form, we will 
have the following equation: 


G, (z) = Bz" + Bz? -Bz! t B, _ Yn 


(2.13) 
Age? -Az ° -Az 41 e, 


Command Window 


>> Cont 
Cont = 
2 1 
Ep 4o ee 
TE s 


with Kp = 3.11e-05, Ti = 2.58e-07 
Continuous-time PI controller in standard form 
>> c2d(Cont,100e-6,'tustin') 
ans = 
h Ts* (z+1) 
Rp p up a ecce ) 
T3 2* (z-1) 
with Kp = 3.11e-05, Ti = 2.58e-07, Ts = 0.0001 


Sample time: 0.0001 seconds 
Discrete-time PI controller in standard form 


fs >> 


FIGURE 2.3 The PI controller transfer function in continuous- and 
discrete-time domain. 
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A, =1, B, = 0.0061, B, = 0.006 (2.14) 


The equivalent equation is given as follows: 


Yn = A3 Yn-3 + A2 Yn-2 + Agi + Bz€n-3 + Bye,» + Bye, 4 Boe, 
(2.15) 


To design a PI controller for voltage control loop, the same as 
the current control loop, we can use pidTuner and c2d tools in 
MATLAB. To design a PI controller for the voltage control loop 
using MATLAB, the following code can be applied. 


s = t£('s'); 

Gvd = (2.257e-5*$4188.1)/ (6e-8*s^240.000291*s 
+7.801); 

pidTuner (Gvd) ; 


By executing the above code, the PID tuner window is popped 
up. We select the type of controller as PI and its form as Standard 
and then regulate the Response Time and Transient Behavior slid- 
ers to get desired Step Plot. Then, if we click on Show Parameters 
button, we can see our controller parameters. Finally, we click on 
the Export button and give it a name for our PI controller, in this 
case, we chose ‘Cv’ as its name and click on the OK button as 
shown in Figure 2.4. 

Now, in the command window, if we type the name (‘Cv’) and 
then press the Enter key, you can get the controller transfer func- 
tion as the following equation: 


aN 


| 

T, : 

6.()=K (1+ ES = 414x106 +22? (2,16) 
s S 
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Amplitude 


LOK Reni Qep |p i 

E 228 G Mes 04 rad; 
82 deg @ 388 rad/s 
Stable 


J 


i L- n 
0 0005 0.01 0015 0.02 0.025 


{2 Controller was re-tuned with the new PID form. Controller Parameters: Kp = 4.143e-06, Ti = 2578e-07 


FIGURE 2.4 The PID tuner window. 
K, 2414x10 5,T; 22.58x107 (2.17) 


If you want to implement analog controller using op-amp, the 
above equation is enough. However, for use in a digital controller, 
the following substitutions must be applied (Tustin's method): 


(2.18) 


where T; is the sampling time and is chosen as 100 us, so we have: 


7 T,(z+1) 
G,(z) = sp BE 
K,(z-1)+K,T, (z+1)/(2T7;) 

(z-1) (2.19) 
(K, - K,T, /(27,))z +K, /(2T;)-K, 
z-1 

0.00081z 4-0.0008 

z-1 
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Command Window 


>> Cv 
Cv = 
1 X 
Epow (£k * eS) 
*X s 
with Kp = 4.14e-06, Ti = 2.58e-07 
Continuous-time PI controller in standard form 
>> c2d(Cv,100e-6,'tustin') 
ans = 
1 Ts* (z+1) 
Kp * (1 + ---- * -------- ) 
Ti 2* (z-1) 


with Kp — 4.14e-06, Ti — 2.58e-07, Ts — 0.0001 


Sample time: 0.0001 seconds 
Discrete-time PI controller in standard form 


fs >> 


FIGURE 2.5 The PI controller transfer function in continuous- and 
discrete-time domain. 


0.0008z ^! +0.00081 
G.(z)= EE (2.20) 


We can also use the c2d function in MATLAB to get the 
discrete-time transfer function from G,(s) in the command win- 
dow, with the code below as depicted in Figure 2.5. 


c2d(Cv,100e-6,'tustin') 


Correlating the equation above with the standard form, we will 
have the following equation: 


G, (z) = B^ + Byz? -Bz! + Bo - Yn 


(2.21) 
Ase -Az °? -Az I €, 
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A, =1, By = 0.00081, B, = 0.0008 (2.22) 


The equivalent equation is given as follows: 


Yn = Az3Yn-3 + A2 Yn-2 + A1 Yn-1 + Bz€n-3 + Bye,» + Bye, 4 + Boe, 
(2.23) 


In continuous conduction mode (CCM), we keep the inductor 
current ripple very low (usually lower than 20%). We select the 
inductor, L, lower than the calculated value in Equation (2.2), for 
practical reasons. Now, we are going to introduce another method 
for controller design. Here, we are going to set the controller as 
simple as possible. We use the phase margin to design the control- 
ler. So, we set up the phase margin that we wanted and designed the 
controller to meet that phase margin at a specific crossover fre- 
quency. The crossover frequency in control systems refers to the 
frequency at which the magnitude of the open-loop transfer func- 
tion is equal to 1 or 0 dB. The switching frequency, on the other 
hand, refers to the frequency at which a switching device (such as a 
transistor in a power converter) turns on and off. There is no direct 
relationship between the crossover frequency and the switching 
frequency in control systems. These are two separate concepts that 
relate to different aspects of a system. However, the choice of 
switching frequency in a power converter can indirectly affect the 
crossover frequency of the control system that is regulating the 
converter. For example, a higher switching frequency in a power 
converter may allow for faster response times in the control system, 
which could potentially impact the crossover frequency. Ultimately, 
the crossover frequency is determined by the design of the control 
system itself and is not directly tied to the switching frequency. In 
this case, the crossover frequency is chosen to be f, = 2 kHz. So, for 
the current controller, we have the following equations: 


s = jo, = j2z f. = j12566.4 (2.24) 
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i, (j12566.4 
(2.6) > PU - 6.2367 114.6 (2.25) 
d(j12566.4) 
Gain = 6.236 (2.26) 
Phase = -114.6 (2.27) 


The following equations are used to determine the compensa- 
tion required to do that we specify a phase margin in crossover 
frequency and calculate the required compensation phase: 


Prem = 59° (2.28) 
Perossover + Peomp +180 = Ppry (2.29) 
Peomp = Pem — Perossover — 180° (2.30) 
Peomp =59 —(-114.6')-180" =-6.4 (2.31) 


We usually choose phase margin between 40 and 60 degrees. 
However, we do not want the system to be unstable, so we have 
selected higher phase margin in this case which is 59 degrees. If 
the compensator phase is positive, you will need a lead-lag com- 
pensator. If the compensator phase is negative, you will need a lag 
compensator, which in this case is a PI controller. Using the above 
equations, a lag compensator is adequate, so the PI compensator is 
of the form: 


G.(s)=K, (241) (2.32) 
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G,(jo.) - K; (2 | (2.33) 

ZG, (ja. ) --—tan'! (2) (2.34) 

6.4 - tan"! | 2 | (2.35) 
0, 

(0; = 0, X tan(6.4°) (2.36) 

@; =12566.4 x tan(6.4 ) = 1409.55rad/s (2.37) 


The value for K, is found using the following equations: 


2 
Qj 1 1 
G, (jo, )|=K,,{| — | +1= e (2.38) 
Ua) =K, (2) Gain 6.236 
Ke i E. (2.39) 


(tan (5.7) +1 


Since, in the crossover frequency, the gain of compensated sys- 
tem should be equal to 0 dB or 1. Therefore, the controller transfer 
function is of the form: 


G.(s) =K, [241] =o.159{ M28 41) —0.1594- 224.12 
s 


s 
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We will have the identical crossover frequency for voltage con- 
trol loop as well. The following equations are applied for the volt- 
age controller. At 2 kHz, the following info is obtained: 


s = jo, = j27 f. = j12566.4 (2.41) 
v, (j12566.4) : 

(2.7) > —— — = 46.797 -114.6 (2.42) 
d(j12566.4) 

Gain — 46.79 (2.43) 

Phase = -114.6 (2.44) 


The following equations are used to determine the compensa- 
tion required to do that we specify a phase margin in crossover 
frequency and calculate the required compensation phase: 


Prem = 59° (2.45) 
Perossover + Peomp +180 = Ppry (2.46) 
Peomp = PEM 7 Perossover ~180° (2.47) 
Peomp =59 -(-114.6°)-180° =-6.4 (2.48) 


We usually choose phase margin between 40 and 60 degrees. 
However, we do not want the system to be unstable, so we have 
selected higher phase margin in this case which is 59 degrees. If 
the compensator phase is positive, you will need a lead-lag com- 
pensator. If the compensator phase is negative, you will need a lag 
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compensator, which in this case is a PI controller. Using the above 
equations, a lag compensator is adequate, so the PI compensator is 
of the form: 


G.(s)=K, (241) (2.49) 
s 
©; 
G.(j@)=K,| —+1 2.50 
Go)» (2s) (2.50) 
ZG. ( jo.) = -an2 (2.51) 
O- 
6.4 =—tan! e ) (2.52) 
©, 
0; =, X tan(6.4°) (2.53) 
@; =12566.4 x tan(6.4 ) = 1409.55rad/s (2.54) 


The value for K, is found using the following equations: 


2 
Q; 1 1 
G,(jo,)-K,]|] — | +1= 2.55 
(o. - K, (2) Gain 46.79 pem 
i L9. - 50931 (2.56) 


(tan(6.4°)) +1 
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Since, in the crossover frequency, the gain of compensated sys- 
tem should be 0 dB or 1. Therefore, the controller transfer func- 
tion is of the form: 


G.(s)=K,| +1]=0.021 PE? d NO ceo 
P| : à 


(2.57) 


We can use Google Colab for Bode plots of closed-loop current 
and voltage control with PI controllers designed by MATLAB pid- 
Tuner tool and phase margin-based frequency response calcula- 
tions. Google Colab is a free cloud-based Jupyter notebook 
environment that allows users to write and execute Python code 
in a browser. It provides access to free GPU and TPU resources, 
making it a popular choice for data scientists and machine- 
learning researchers to work on computation-intensive tasks. 
Google Colab also allows for easy sharing and collaboration on 
projects, the Google Colab code is as below. 


!pip install control 


$matplotlib inline 

import numpy as np 

import matplotlib.pyplot as plt 

from control.matlab import * 

Gid = tf([3e-6, 25.07],[7.33e-8, 0.0002912, 
7.8011) 

Gvd = t£f([2.257e-5, 188.1],[7.33e-8, 0.000291, 
7.801] ) 


print ("Gid = ", Gid) 
print ("Gvd = ", Gvd) 
Gid = <TransferFunction>: sys[0] 


Inputs (1): ['ul0]'] 
Outputs (1): ['yl0]1'] 
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3e-06 s + 25.07 


7.33e-08 s^2 + 0.0002912 s + 7.801 


Gvd = <TransferFunction>: sys[1] 
Inputs (1): ['u[0]'] 
Outputs (1): ['yl0]']l 


2.257e-05 s + 188.1 


7.33e-08 s^2 + 0.000291 s + 7.801 


The above codes are used for printing the transfer functions of 
Gid and Gvd as given in Equations (2.6) and (2.7). The code below 
is used for the Bode plot of Gid as illustrated in Figure 2.6. 


w - np.logspace(1,5) 


mag,phase,omega = bode (Gid,w) 
plt.tight layout () 


Bode plot for sys[0] 
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3 
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FIGURE 2.6 The Bode plot of Gid. 
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The Gci_mat is the transfer function of current controller 
obtained by the pidTuner tool of MATLAB as given in the code 
below. 


Gci mat = tf([3.11e-5, 120.5],[1, 01) 
print("Gci mat - ", Gci mat) 

Gci mat = <TransferFunction>: sys[3] 
Inputs (1): ['ul[0]'] 

Outputs (1): ['yl01'] 


3.11e-05 s + 120.5 


The feedback loop transfer function is the product of the plant 
and controller transfer functions and the following code is used 
for the Bode plot of it as shown in Figure 2.7. 


Bode plot for sys[4] 


Magnitude [dB] 
I 
N 
V 


Phase [deg] 


10! 10? 10° 10° 10° 
Frequency [rad/s] 


FIGURE 2.7 The Bode plot of Gpci mat. 
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Bode plot for sys[15] 


Magnitude [d8] 


Phase [deg] 


10! 10? 10 10! 10° 
Frequency [rad/s] 


FIGURE 2.8 The Bode plot of Tpci_mat. 


Gpci mat = Gid*Gci mat 

w - np.logspace(1,5) 
mag,phase,omega = bode(Gpci mat,w) 
plt.tight layout () 


The closed-loop transfer function and the code for the Bode 
plot of it are depicted in Figure 2.8. 


Tpci mat = Gpci mat/(1«Gpci mat) 

w - np.logspace(1,5) 
mag,phase,omega - bode(Tpci mat,w) 
plt.tight layout () 


TheGci pmisthetransfer function ofcurrent controller obtained 
by the phase margin calculations as given in the code below. 


Gci pm = tf([0.159, 224.12],[1, 01) 
print ("Gci_pm = ", Gci pm) 
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Gci pm = <TransferFunction>: sys[10] 
Inputs (1): ['u[0]'] 
Outputs (1): ['yl0]1'] 


0.159 s + 224.1 


The feedback loop transfer function is the product of the plant 
and controller transfer functions and the following code is used 
for the Bode plot of it as shown in Figure 2.9. 


Gpci pm - Gid*Gci pm 

w - np.logspace(1,5) 
mag,phase,omega - bode(Gpci pm,w) 
plt.tight layout () 


Bode plot for sys[11] 
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FIGURE 2.9 The Bode plot of Gpci_pm. 
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Bode plot for sys[18] 
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FIGURE 2.10 The Bode plot of Tpci pm. 


The closed-loop transfer function and the code for the Bode 
plot of it are depicted in Figure 2.10. 


Tpci pm = Gpci pm/(1«Gpci pm) 

w - np.logspace(1,5) 
mag,phase,omega - bode(Tpci pm,w) 
plt.tight layout () 


The code below is used for the Bode plot of Gvd as illustrated in 
Figure 2.11. 


w = np.logspace(1,5) 
mag, phase,omega = bode (Gvd,w) 
plt.tight layout () 
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Bode plot for sys[1] 
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FIGURE 2.11 The Bode plot of Gvd. 


The Gcv_matis the transfer function ofcurrent controller obtained 
by the pidTuner tool of MATLAB as given in the code below. 


Gcv mat = tf([4.14e-6, 16.05],[1, 01) 
print("Gcv mat - ", Gcv mat) 

Gcv mat = <TransferFunction>: sys[42] 
Inputs (1): ['ul[0]'] 

Outputs (1): ['yl0]'] 


4.14e-06 s + 16.05 


The feedback loop transfer function is the product of the plant 
and controller transfer functions and the following code is used 
for the Bode plot of it as shown in Figure 2.12. 
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Bode plot for sys[44] 
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FIGURE 2.12 The Bode plot of Gpcv_mat. 


Gpcv_mat = Gvd*Gcv_mat 

w = np.logspace (1,5) 

mag, phase,omega = bode (Gpcv mat,w) 
plt.tight layout () 


The closed-loop transfer function and the code for the Bode 
plot of it are depicted in Figure 2.13. 


Tpcv mat = Gpcv mat/(1«Gpcv mat) 

w - np.logspace(1,5) 
mag,phase,omega - bode(Tpcv mat,w) 
plt.tight layout () 


The Gcv pm is the transfer function of current controller 
obtained by the phase margin calculations as given in the code 
below. 


Gcv pm = tf£([0.021, 29.6],[1, 01) 
print("Gcv pm - ", Gcv pm) 
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Bode plot for sys[48] 
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FIGURE 2.13 The Bode plot of Tpcv. mat. 


Gcv pm = <TransferFunction>: sys[50] 
Inputs (1): ['u[0]'] 
Outputs (1): ['yl0]1'] 


0.021 S + 29.6 


The feedback loop transfer function is the product of the plant 
and controller transfer functions and the following code is used 
for the Bode plot of it as shown in Figure 2.14. 


Gpcv pm - Gvd*Gcv pm 

w - np.logspace(1,5) 
mag,phase,omega = bode(Gpcv pm,w) 
plt.tight layout () 


The closed-loop transfer function and the code for the Bode 
plot of it are depicted in Figure 2.15. 
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Bode plot for sys(51] 
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FIGURE 2.14 The Bode plot of Gpcv_pm. 


Bode plot for sys[55] 
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FIGURE 2.15 The Bode plot of Tpcv. pm. 
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Tpcv pm = Gpcv pm/(1«Gpcv pm) 

w = np.logspace(1,5) 
mag,phase,omega = bode (Tpcv_pm, w) 
plt.tight layout () 


As you can see in Bode plots in Google Colab for closed-loop 
current and voltage control using MATLAB pidTuner tool, and 
phase margin calculations, the designed PI controller for current 
and voltage control loops is stable. For further analysis, we can 
compare the operation of designed PI controllers in MATLAB 
Simulink. As illustrated in Figure 2.16, we are going to compare 
the step or transient response of closed-loop current and voltage 
control with PI controllers designed by MATLAB pidTuner tool 
and phase margin calculations. 

The results of transient response of current and voltage control 
loops with PI controllers designed by MATLAB pidTuner tool and 
phase margin calculations are plotted on the scopes as shown in 
Figures 2.17 and 2.18, respectively. 


€) 3.116-55*120.5 3e-65*25.07 
s 66-8s2«0.00029125*7.8 


Step | Transfer Fen Transfer Fen1 


0.159s*224.12 
s 


Transfer Fcn2 


Current Control Loop 
using MATLAB/ 
Phase Margin 


3e-6s+25.07 
6e-8s2+0.0002912s+7.8 


Transfer Fcn3 


4.14e-6s+16.05 p| —2257e-5s+188.1 
Q s 6e-852+0.0002912s+7.8 
Transfer Fcn4 Transfer FenS 


0.0215*29.6 p| — 2257e55*188.1 
s 6e-852+0.0002912s+7.8 | 


>i 


Voltage Control Loop 
using MATLAB/ 
Phase Margin 


> 


Transfer FcnG Transfer Fcn7 


FIGURE 2.16 Comparing the transient response of current and voltage 
control loops with PI controllers designed by MATLAB pidTuner tool 
and phase margin calculations. 
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FIGURE 2.17 The step response of current control loop (up: MATLAB, 
down: Phase Margin). 


File Tools View Simulation Help * 
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FIGURE 2.18 The step response of voltage control loop (up: MATLAB, 
down: Phase Margin). 
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Although the PI controllers with both methods are stable, how- 
ever, the PI controllers designed by pidTuner tool of MATLAB 
have better transient response than phase margin calculations. 
Therefore, we will use the PI controllers designed by pidTuner tool 
of MATLAB. 


2.2 THE CONTROLLER DESIGN WITH OCTAVE 


In this section, we are going to discuss about how to obtain 
the digital controller coefficients for designing closed-loop 
buck converter controllers using STM32 microcontroller. To 
do that, we can also apply the Octave software for designing 
the PI controller, using phase margin calculations. Therefore, 
you should first initialize the buck converter parameters as the 
code below. 


sclear all 

pkg load control 

Sz tfttt's) 
i 


i 
Initialization of converter parameters 


Ro = Vo/Io; 

fsw = 100000; 

rL = 300e-3; 

rESR = 3e-3; 

Vds = 0.075; 

Vfwd = 1.15; 

ib - Io; 

D = (Vfwd + iL*rL + Vo)/(Vin - Vds + Vfwd); 
L = 200e-6; 


C = 40e-6; 
delta_iL = 0.1 * iL; 
delta_Vo = 0.02 * Vo; 


E! 


>= (Vin - Vds - iL*iL - Vo)*D/(fsw*delta iL); 
S 


o? o9 oe 
Q 
V 
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The parasitic values are inductor resistance, r,, equivalent series 
resistance of capacitor, fgsp saturation voltage of MOSFET, V,,, and 
forward voltage drop of diode, Vwa which have very little effect on 
output of the system, but it is useful to include them in practice. All 
of the commands in Octave can be executed in MATLAB as well. 
Then we initialize the crossover frequency which is very important 
in controller design. In power supply design, we mostly focus on 
crossover frequency of the converter frequency response. So, in 
power supply design, we are focusing on frequency-domain. If you 
choose the crossover frequency high, your converter will be very 
sensitive and it will respond basically to anything. So, we here do 
not focus on disturbances, however, we will mostly focus on robust- 
ness, so, if we change the set point, it will not affect the stability of 
the converter. We define crossover frequencies for current and volt- 
age separately as given in the following code. 


$Initialization of the crossover frequency 
9- 9- 9- 9- 9- 9- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- O- 9- 9- O- 9- O- O- 9- 9- O- 9- 9- O- 9- 9- O- 9- O- 
O "O 7O 'O 7O 7O 'O 7O 'O 7O 7O 'O OOOO 7O 7O 7O 7O OOOO OOOO OOOO OOOO OOOO OO 
wc i = 2*pi*2000; 

fc i = wc i/(2*pi); 


wc complex i - wc i*li; 
wc v = 2*pi*2000; 

fc v = wc v/(2*pi); 

wc complex v = wc v*li; 
targ PM deg i - 5 
targ PM deg v = 5 


The codes for small signal output current to duty cycle and out- 
put voltage to duty cycle transfer functions are as below. 


n 
3 
D 


11 Signal Transfer Functions 


ae oe 
o? oe 
o? oe 
oe? oe 
oe? oe 
oo 
oo 
oo 
oo 
oo 
oo 
oe 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oe 
oo 
oo 
oo 
oo 


o? oe 
o? oe 


o? oe 
oe? oe 


o? oe 
o? oe 
o? oe 


o? oe 
oe oe 


o? oe 
o? oe 
oe? oe 
oe? oe 


o? oe 
oo? oe 


o? oe 
oe? oe 


o? o9? o9? oe 
o? oe 
o? oe 
o? oe 
oP oe 
o? oe 
o? oe 


Output Current to Duty Cycle beginning 
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GIoD num = (1«s*C*rESR)*(Vin - Vds + Vfwd); 
GIOD den (L*C*Ro + L*C*rESR)*s^2 + (C*Ro*rL 

+ C*rESR*iL + C*Ro*rESR + L)*s + rL + Ro; 

GIoD = (GIoD num/GIoD den); 

$Output Current to Duty Cycle ending 

$Output Voltage to Duty Cycle beginning 

GVoD num = (1-«s*C*rESR)*(Vin - Vds + Vfwd)*Ro; 
GVoD den (L*C*Ro + L*C*rESR)*s^2 + (C*Ro*rL + 
C*rESR*rL + C*Ro*rESR + L)*s + rL + Ro; 

GVoD = (GVoD_num/GVoD_den) ; 


O 
& 


o? o9 oe 


The realization of complex small signal output current to duty 
cycle and output voltage to duty cycle transfer functions at the 
crossover frequency is represented in the following code. 


$Transfer functions to realize the gain and 
phase at crossover frequency 


oe oe 


oe oe 
oe oe 


oe oe 
oe oe 


o? o9? oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 


Output Current to Duty Cycle gain phase 
beginning 

g - wc complex i; 

GIoD num g = (1«g*C*rESR)*(Vin - Vds + Vfwd); 
GIoD den g- (L*C*Ro + L*C*rESR)*g^2 + (C*Ro*rL 
+ C*rESR*rL + C*Ro*rESR + L)*g + rL + Ro; 
GIoD_g= (GIoD num g/GIoD den g); 

sOutput Current to Duty Cycle gain phase ending 


6 


o 


sOutput Voltage to Duty Cycle gain phase beginning 
g = wc complex v; 

GVoD num g = (1«g*C*rESR)*(Vin - Vds + Vfwd)*Ro; 
GVoD den g = (L*C*Ro+ L*C*rESR)*g^2 + (C*Ro*rL + 
C*rESR*iL + C*Ro*rESR + L)*g + rL + Ro; 

GVoD g = (GVoD num g/GVoD den g); 


o 


sOutput Voltage to Duty Cycle gain phase ending 


110 m Digital Control of Power Converters 


The code required for current and voltage controller design is 
given as follows. 


$Current Controller design begin 

i gain = abs(GIOD g); 

i phase = angle (GIoD g)*(180/pi); 

targ PM rad i - targ PM deg i*(pi/180); 

phase comp i = targ PM rad i - (i phase*(pi/180)) 
- pi; 

fi i = fc i * tan(-1*phase comp i); 

Kpi = (1/i gain)/sqrt((fi i/fc i)^2 + 1); 

wi i = fi i*2*pi; 

$Current Controller design ends 

$Voltage Controller design begin abs (GVOD g); 
v gain - abs(GVoD g); 

v phase = angle (GVoD g)*(180/pi); 

targ PM rad v - targ PM deg v*(pi/180); 

phase comp v - targ PM rad v - (v phase* (pi/180)) 
- pi; 

fi v = fc v * tan(-1*phase comp v); 

Kpv = (1/v gain)/sqrt((fi v/fc v)^2 + 1); 

wi v - fi v*2*pi; 

$Voltage Controller design ends 


The codes for current and voltage controller transfer functions 
and also current and voltage closed-loop transfer functions are 
illustrated below. 


$EZSELEESESLAELLELELLLELLEELAELLAELESEELESEELEA ELLA EE 
EESEEEEEEEEEEEEEEEEEEEEEEEEES 

Gci = Kpi*(wi i/s + 1); 

Gcv = Kpv* (wi v/s + 1); 


tGcl = 1s 
Ti = Gci*GIoD; 
Tv = Gcv*GVOD; 
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$Create a figure to hold the plots 


sfigure; 

$Plot Bode plot for the system 

sbode (Ti); 

sbode (Tv) ; 
UIILIELIEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEREEE 
SESECEEEEESEEESEECESEEEEEEESEES 


The following codes are used for calculating gain and phase 
margins using closed-loop transfer functions. 


$Calculate gain and phase margins 

[gain margin, phase margin, crossover _ 
frequencies, phase crossover] - margin(Ti); 
$Display Current Bode Plot the results 
fprintf('Current, Gain margin: $.2f dB\n’, 
20*1og10 (gain margin)); 

fprintf('Current Phase margin: $.2f degrees \n’, 
phase margin); 

fprintf('Current Crossover frequencies: $.2f 
rad/s\n’, crossover frequencies); 
fprintf('Current Phase at crossover: $.2f 
degrees\n’, phase crossover); 

[gain margin, phase margin, crossover _ 
frequencies, phase crossover] - margin(Tv); 
$Display Voltage Bode Plot the results 
fprintf('Voltage Gain margin: $.2f dB\n’, 
20*1og10 (gain margin)); 

fprintf('Voltage Phase margin: $.2f degrees\n’, 
phase margin); 
fprintf('Voltage Crossover frequencies: $.2f 
rad/s\n', crossover frequencies); 
fprintf('Voltage Phase at crossover: $.2f 


, . 
degrees\n’, phase crossover); 
00.00.00000000000000000000000000000000000000000000 
$$1515t55555590505000000500000050000050000005000095 050000950 
000.00000000000000000000000000 
$59515t5655559900000000500000500000 
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For converting controller transfer functions from continuous 
to digital mode or from s-domain to z-domain, we use the sam- 
pling time of 100 us as follows. 


Q 
H- 
Q 
ll 


= c2d(Gci, 100e-6); 
Gv d = c2d(Gcv, 100e-6); 


So, in this case, the duty cycle is updated every 100 us. We run 
the program and go to the command window and get current and 
voltage closed-loop gain and phase margins with current and 
voltage controllers with and without current and voltage control- 
lers and then get the controllers transfer functions in z-domain as 
below. 


>> loop design octave 

Current, Gain margin: Inf dB 

Current Phase margin: 59.06 degrees 

Current Crossover frequencies: NaN rad/s 
Current Phase at crossover: 12563.57 degrees 
Voltage Gain margin: Inf dB 

Voltage Phase margin: 58.94 degrees 

Voltage Crossover frequencies: NaN rad/s 
Voltage Phase at crossover: 12569.17 degrees 
>> Gel 

Transfer function 'Gci' from input ‘ul’ to output 


0.1594 s + 225.5 


YL :o9g————Óeg-oR- 

S 
Continuous-time model. 
>> Gcv 


Transfer function 'Gcv' from input ‘ul’ to output 


0.02127 S + 30.16 
Vie seme) Seas ss 


Continuous-time model. 
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>> Gid 
Transfer function 'Gi d' from input 'ul' to 
output .. 
0.1594 z - 0.1369 
Vis. S PSCCRSEUSASCEE 


Sampling time: 0.0001 s 
Discrete-time model. 


>> Gv_d 
Transfer function ‘Gv_d’ from input ‘ul’ to 
output .. 
0.02127 z - 0.01825 
Vis... BemeeeEeREEeR&CIÓM 


Sampling time: 0.0001 s 
Discrete-time model. 


These z-domain transfer functions are actually what we want 
to use in our codes. The ‘J, value in the initialization of converter 
parameters at the beginning does not affect so much on the 
designed robust system. So, the output current set point could be 
changed in the firmware as you like in the reasonable range and it 
will not make the system unstable. As mentioned before, we will 
use the PI controllers obtained by pidTuner tool of MATLAB; 
however, we can also use the PI controllers gained from phase 
margin calculations by Octave programming. The whole code in 
Octave software is as below. 


$clear all 

pkg load control 

tes) 4 

Initialization of converter parameters 


n 
I 


wm 
[9] 
Il 
< 
O . 
™~ 
H 
[9] 


0.1 * iL; 
0.02 * Vo; 


(Vin - Vds - iL*iL - Vo)*D/(fsw*delta iL); 


100000; 
3e-3; 


delta Vo 


(V£wd + iL*rL + Vo)/(Vin - Vds + Vfwd); 


200e-6; 
40e-6; 


delta iL 


fsw 
$C >= delta iL/(8*delta Vo*fsw); 
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rESR 
D 


oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oe 
ole 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
ole 
oe 
ole 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oe 


$Initialization of the crossover frequency 


oo 
oo 
ole 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
ole 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
ole 
oo 
oo 
oo 
oo 
ole 
oo 
oo 
oo 
oo 
ole 
oo 
oo 
oo 
oo 
oo 
oec 


wc i*li; 
wc v*1li; 
59; 
targ_PM_deg_v = 59; 


2*pi*2000; 


we_v/(2*pi) ; 
wc complex v 


2*pi*2000; 
wc i/(2*pi); 


wc i 
fc x 
wc complex i 
WC V 
fc v 
targ PM deg i 


oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 


$Small Signal Transfer Functions 


oo 
oo 
oo 
oo 
oo 
oo 
oe 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oe 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 


oe 
ole 
ole 
oo 
oo 
oo 
ole 
oo 
oo 
ole 
ole 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
oo 
ole 
oe 
oo 
oo 
oo 
oo 
oo 
oo 


(C*Ro*rL + 
(C*Ro*rL + 


o Duty Cycle ending 


o Duty Cycle beginning 
(1+s*C*rESR)* (Vin - Vds + Vfwd); 


(L*C*Ro « L*C*rESR)*s^2 « 


C*rESR*iL + C*Ro*rESR + L)*s + rL + Ro; 
o Duty Cycle beginning 


(1+s*C*rESR)* (Vin - Vds + Vfwd) *Ro; 


(L*C*Ro + L*C*rESR)*s^2 + 


C*rESR*rL + C*Ro*rESR + L)*s + rL + Ro; 


(GIoD num/GIoD den) ; 


GIoD_den 


GVoD den 


$Output Voltage 


$Output Current 
GVoD num 


$Output Current 
GIOoD num 
GIoD 
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GVoD = (GVoD_num/GVoD_den) ; 
$Output Voltage to Duty cycle ending 


oe oe 
oe oe 


oe oe 
oe oe 


o? o9 oe c 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 


Transfer functions to realize the gain and 
phase at crossover frequency 


9- 9- 9- 9- 9- O- 9- 9- 9- 9- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- 9- O- O- 9- O- O- 9- - O- 9- 9- O- 9- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- O- 
O 'O 'O 'O 'O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O OOOO OOOO ^O 7O 7O 7O 7O 7O 7O 7O 7O 7'O 7'O 'O OOOO 'O 
9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- 9- O- 9- 9- O- 9- 9- O- 9- 9- O- 9- O- O- O- 

O 'O 'O 'O 7O 'O 'O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O 7O OOOO OOO DO 

sOutput Current to Duty Cycle gain phase 
beginning 

g - wc complex i; 

GIoD num g = (1«g*C*rESR)*(Vin - Vds + Vfwd); 
GIoD den g = (L*C*Ro + L*C*rESR)*g^2 + (C*Ro*rL 
+ C*rESR*rL + C*Ro*rESR + L)*g + rL + Ro; 

GIOD g = (GIOD num g/GIoD den g); 


sOutput Current to Duty Cycle gain phase ending 
$Output Voltage to Duty Cycle gain phase 
beginning 

g = wc complex v; 

GVoD num g = (1«g*C*rESR)*(Vin - Vds + Vfwd)*Ro; 
GVoD den g = (L*C*Ro+ L*C*rESR)*g^2 + (C*Ro*rL + 
C*rESR*iL + C*Ro*rESR + L)*g + rL + Ro; 

GVoD g = (GVoD num g/GVoD den g); 

$Output Voltage to Duty Cycle gain phase ending 


o? oe 


oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
o? oe 


oe oe 
o? oe 


oe oe 
oe oe 


oe oe 
oe oe 


oe oe 
o? oe 


o? o9 oe c 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 
oe oe 


Current Controller design begin 

i gain = abs(GIoD_g); 

i phase = angle (GIoD g)*(180/pi); 
targ PM rad i - targ PM deg i*(pi/180); 


phase comp i - targ PM rad i - (i phase*(pi/180)) 
- pi; 

fi i - fc i * tan(-1*phase comp i); 

Kpi = (1/i gain)/sqrt((fi i/fc i)^2 + 1); 

wi i - fi i*2*pi; 


$Current Controller design ends 
$Voltage Controller design begin abs (GVOD g); 
v gain - abs(GVoD g); 
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v phase = angle (GVoD g)*(180/pi); 
targ PM rad v - targ PM deg v*(pi/180); 


phase comp v - targ PM rad v - (v phase* 
(pi/180)) - pi; 

fi v = fc v * tan(-1*phase comp v); 

Kpv = (1/v gain)/sqrt((fi v/fc v)^2 + 1); 
wi v - fi v*2*pi; 


Gci = Kpi*(wi i/s + 1); 
Gcv = Kpv* (wi v/s + 1); 
*$Gcl e l1; 

Ti = Gci*GIOD; 

Tv = Gcv*GVOD; 


$Create a figure to hold the plots 
$figure; 
$Plot Bode plot for the system 


sbode (Ti); 

sbode (Tv) ; 
SESESSEEEESEEEESESEEEESEEESSESEEEESEEESESEEESESSEES 
ESSESEESEESSESEES EES EESSEESSEES 

$Calculate gain and phase margins 

[gain margin, phase margin, crossover _ 


frequencies, phase crossover] - margin(Ti); 
$Display Current Bode Plot the results 
fprintf('Current, Gain margin: $.2f dB\n’, 
20*1og10 (gain margin)); 

fprintf('Current Phase margin: $.2f degrees \n’, 
phase margin); 

fprintf('Current Crossover frequencies: $.2f 
rad/sNn', crossover frequencies); 
fprintf('Current Phase at crossover: %.2f 
degrees\n’, phase crossover); 

[gain margin, phase margin, crossover 
frequencies, phase crossover] - margin(Tv); 
$Display Voltage Bode Plot the results 
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degrees\n’, phase crossover); 


100e-6); 


c2d(Gcei, 
Gv d = c2d(Gcv, 


Gid 


100e-6); 


CHAPTER 3 


Digital Control 
Implementation of 
the Buck Converter 


3.1 HARDWARE OUTLINE 


In this chapter, we are going to discuss about the hardware we will 
use to implement the buck converter controller. We will utilize 
the STM32 Nucleo-G474RE board and it is connected to the main 
converter board as shown in Figure 3.1. 

We have actually an AC-DC converter that converts single 
phase 230 V,. to 24 Va. and that is an input to the system. In the 
main converter, we have an auxiliary power supply and, of course, 
we have an inductor coil, and on the right-hand side of it, we have 
a diode and on the left-hand side, there is a Metal-Oxide- 
Semiconductor Field-Effect Transistor (MOSFET) and the 
MOSFET is significantly over sized and it does not require a heat 
sink, but the diode requires a heat sink. Then we have the output 
of the main converter and the resistive load is connected to the 
output. There are four 10 uF parallel capacitors and also we have 
an inductor that serves outgoing current to the output and it 
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FIGURE 3.1 The buck converter controller implementation. 


comes after the capacitors. Also, we have a shunt current-sensing 
resistor next to the inductor. There is another shunt current- 
sensing resistor of the main inductor. On the top right-hand side, 
we have a voltage sensor for the output. There is also a MOSFET 
driver in this case and we usually use bootstrap MOSFET driver 
for this particular application but we have not used it in this main 
converter board because they have some drawbacks, so, we use an 
isolated MOSFET driver. The converter is not considered isolated 
since the ground of input and output are the same. However, it 
helps go to any duty cycle we want. Next to the MOSFET driver, 
there is a power supply for it. It is a very common converter for 
this application. The STM32 Nucleo-G474RE board has high- 
resolution PWM and timer modules which are very useful in 
power electronics applications. The codes provided can be used 
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FIGURE 3.2 The schematic of buck converter and its controller. 


on other STM32 microcontrollers as well. The buck converter and 
its controller schematic using the STM32 Nucleo-G474RE devel- 
opment board are depicted in Figure 3.2. 

We have two current sensors and two voltage sensors in the 
schematic. The input voltage and current sensors are not used in 
this system. However, the output voltage and current sensors are 
utilized. We have also two input and output inductor current sen- 
sors and we can switch between them. We have used HCPL-3120 
as the MOSFET driver. There is an auxiliary power supply, 
LM2596S, which generates 5 V at the output and its maximum 
input voltage is 28 V. There is also a power supply, B0515S, exclu- 
sively for the MOSFET driver which converts 5 V-15 V. The only 
deviation from traditional buck converter is that there is an out- 
put inductor, 10 uH, which is an outgoing for the next section of 
buck converter, for the load basically. We have also connections 
from buck converter board to the STM32 Nucleo-G474RE 
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development board. The current sensor amplifier IC is INA180 
from Texas Instruments Incorporated (TT). 


3.2 FIRMWARE PERIPHERAL INITIALIZATION 
WITH CUBEMX 


In this chapter, we are going to initialize the firmware periph- 
erals to program the buck converter controller which is STM32 
Nucleo-G474RE development board. To get started, you should 
note that the initialization is done using CubeMX software and it 
is done easily and already available with that. Therefore, we open 
the STM32CubeMX and under File tab, we create a New Project. 
In the Board Selector tab and in the Commercial Part Number, 
we do type NUCLEO-G474RE and select it in the bottom window 
and click on the Start Project button. In the popped-up window, 
we click on the Yes button to initialize all peripherals with their 
default mode to load the IOC and rendering UI. We select the 
Clock Configuration tab, the maximum clock frequency is 170 
MHz. You can use pre-scaler to reduce the clock speed. However, 
in this case, we want to get maximum PWM resolution to real- 
ize the closed-loop control, so the maximum clock frequency 
is perfect. Then we go to the Project Manager tab, we do copy 
and paste the project directory in the Project Location also select 
a name, that iss BUCK. CONTROL in the Project Name. Then 
select STM32CubeIDE as Toolchain/IDE. We do not need to 
change anything in the Tools tab. Then we go back to the Pinout 
& Configuration tab. The STM32 microcontroller is capable of 
doing many tasks and we are going to initialize what is actually 
needed for the system. We start with the System Core and go to 
the DMA; however, before using DMA, the DMA is actually can 
be used in conjunction with the ADC and in the most conve- 
nient way in closed-loop controller is done by DMA. Then, in the 
ADCI Mode and Configuration window, we set INI and IN2 as 
single-ended. 

Then, we go back to the DMA and click on the Add button and 
select ADCI item and choose Circular as Mode and Word as 
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Data Width. We do not discuss deeply any available peripherals, 
so we only consider peripherals as initialization perspective. 
Then we go back to the ADCI Configuration window and in the 
Parameter Settings tab, we select independent mode for Mode, 
synchronous clock mode divided by 4 for Clock Pre-scaler, ADC 
12-bit resolution, Data Alignment is not important, scan conver- 
sion mode is disabled, end of sequence conversion since we have 
two ADC channels, continuous conversion mode is enabled, and 
DMA continuous requests is enabled. In the ADC Regular 
ConversionMode part, we do enable regular conversions and 
oversampling is disabled. In the number of conversions, we 
choose 2 since we have 2 channels and for the sampling time of 
Channel 1, we select the maximum cycle, because it actually 
makes sure that the system is much stable so we sacrifice the sam- 
pling time to get much stability, so the Rank 1 is completed. We 
do the same settings for Rank 2 and Channel 2. Then we come 
back to the scan conversion mode in ADC Settings part and you 
can seethat is enabled now. We leave the other settings unchanged 
for ADCI. The next one is actually the timers. So, we go to the 
TIMI in Timers section and this is specifically for PWM signal 
generation. We choose Internal Clock for Clock Source and select 
PWM Generation CHI for Channel 1. You have not inversed out- 
put, so you select CH1 and not CHIN where its polarity is 
inversed relative to CH1 and you can see the TIMI CHI is actu- 
ally in the pin PCO in this case. In the Configuration window and 
Parameter Settings tab, we choose Counter Mode as Up and we 
select No Division in Internal Clock Division (CKD) for maxi- 
mum resolution, also we do enable the auto-reload preload and 
you do not need to worry about the Break and Dead Time 
Management, it will play the role in the synchronous buck, boost 
or half-bridge converters and you need to consider the dead time 
management in those aspects. However, in this case, we will 
focus on Pre-scaler and Counter Period. The Counter Period and 
Pre-scaler are from 0 to 65535 (2!° — 1). The following equations 
can be used in this case: 
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E Fox 
is (ARR +1)x(PSC +1) oe 
Fax 
RR = 2 
^ (PSC +1)x(Fowm ) i A 


The Auto Reload Register (ARR) or Counter Period value is 
equivalent to the compare value and PSC is the Pre-scaler value. 
In this case, we assume the PWM frequency is 100 kHz and the 
clock frequency is 170 MHz. If, for example, the PSC is equal to 1, 
from Equation (3.2), we can calculate ARR which is 849. The next 
one is the timer interrupts, remember that in Chapter 2, to gener- 
ate the digital controller coefficients, we selected the timer inter- 
rupts to be 100 us which means its frequency is actually 10 kHz. 
We use TIM4 in this case to implement the timer interrupt. The 
first thing is to select Internal Clock as Clock Source and we will 
not have any outputs in this case, so, we will not select any chan- 
nels. To do the parameter settings, we should note that the equa- 
tion in this case is slightly different as below: 


= Fax 
ARR = (PSC sn) 1 (3.3) 


The Fi; in this case is 10 kHz, and the same procedure done 
before is followed here. We select the PSC which is 1 in this case, 
also the Ferg is 170 MHz, so the calculated ARR is to be 8499 and 
we keep the auto-reload preload to be Disable. So, everything is 
perfect and we can click on the Generate Code button. After gen- 
erating the code, the STM32CubeIDE is automatically opened. So, 
the successful code generation message window is popped up 
finally, and then we click on the Open Project button. So, we 
launch the program in the workspace directory. Finally, you can 
see successfully importing the program into the workspace 
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message. Then, we click on the project name to open it in 
STM32CubeIDE. We go to the main.c file in Src folder and see 
that the selected modules have been initialized. If you have forgot- 
ten to add some initialization modules, you can click on the 
BUCK CONTROL .ioc file and in the popped up Open Associated 
Perspective window, click on the Yes button. Then, a new window 
called the Question window is popped up and we click on the No 
button to load the IOC. You may have mistaken and it is possible 
and that is basically what you need for initialization. In the next 
section, we will dive into coding the firmware and several other 
aspects that need to be taken into account. 


3.3 MANUAL FIRMWARE INITIALIZATION 


In this section, we are going to implement the second part of ini- 
tialization of the firmware. This part is not done by CubeMX but 
it is very necessary for the functionality of the closed-loop control 
system. The first thing is the initialization of PWM and the timer 
actually we have done PWM with timer. In the previous section, 
we did timer interrupts on TIM4 and did the parameter settings, 
and then we go to the NVIC settings and enable the TIM4 global 
interrupt. We just save the changes and generate the code gain. 
Before doing anything, you should choose if the system giving 
you the PWM with the frequency that is actually required and, in 
this case, we have set it to be 100 kHz. The library we are using is 
the HAL library and it means Hardware Abstraction Layer and it 
is so convenient to use. Then, we go to start the PWM by writing 
the following code before the while (1) code: 


HAL TIM PWM Start (&htiml, TIM CHANNEL 1); 


Then, we are going to initialize the PWM duty cycle. We calcu- 
lated PWM duty cycle to be between 0 and 849. We start with a 
duty cycle of 350 which is about 41.2% duty cycle by adding the 
following code inside the while (1) loop and run the program as 
depicted in Figure 3.3. 


TIM1->CCR1 = 350; 
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FIGURE 3.3 Adding code to start the PWM with a duty cycle of 41.296. 


If you connect pin PCO to the oscilloscope, you can see the 
PWM signal with a frequency of 100 kHz and a duty cycle of 
41.296. We will use a maximum of 7096 of duty cycle in this pro- 
gram. Next, we are going to initialize the ADC and DMA which is 
actually done together in this case. As a first thing, we will utilize 
the ADC with DMA and add the following code before while (1): 


HAL ADC Start DMA(&hadcl, (uint32 t *) 

adcValues, 2); 

Also, in the Private variables section, you need 
to set up the code below. 

float Vactual 0.0; 

float Iactual 0.0; 

uint32 t adcValues[2]; 


ll 


It is better to initialize the timer interrupt as well so we add the 
following code before while (1). 


HAL TIM Base Start IT(&htim4) ; 


That code actually enables the timer interrupts to get started, 
the TIM4 counter is counting up and when it is overflowed it can 
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trigger the interrupt duty. We need a section to do the interrupt 
duty, so we just go to the bottom of the program in USER CODE 
4 section and add the following code for interrupt duty. 


void HAL TIM PeriodElapsedCallback(TIM 
HandleTypeDef* htim)( 
if(htim-»Instance == TIM4) { 
Vactual = adcValues [0] *0.008137207031; 
Iactual = adcValues [1] *0.001611328125; 


Once, initializing the timer interrupts, it starts to counting 
up at the beginning of initialization and when it overflows it 
executes the above code and then checks the Instance, if it is 
TIM4, then it is triggered and then it executes the codes inside 
of it. As you know, we have used resistor dividers to sense input 
or output voltages, there are applications such as solar or bat- 
tery charge controller, we are interested to sense the input volt- 
age; however, in this case, we have used output voltage sensing 
for our system. The following equations are used to calculate the 
V 


actual: 


RawAdcVal RR 

Vctual = (ee) x Vat ) x tR) (3.4) 
RawAdcVal 

tae eee eae eee (3.5) 

4096 100k 
Vactual = RawAdcVal x M x R + R, (3.6) 
2 icCKes R, 
Vias = RawAdcVal x Ze x ark (3.7) 


4096 100k 
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Vactual = RawAdcVal x 0.008137207031 (3.8) 


For the I 


actual? 


we have something similar: 


2 AdcRes 


iru Va 


AmpGain 


I actual — 


Rshunt ( 3. 9) 


RawAdcVal 
———————— Ix3.3 
4096 
50 


Tni = 0.01 (3.10) 


You should note that the voltage read in the STM32 microcon- 
troller ADC pin is as below: 


yos SONOS o (3.11) 
4096 


For calculating the actual voltage across the shunt resistor, we 
should divide the voltage read in STM32 microcontroller ADC pin 
by the amplifier gain, AmpGain, which is for INA180 differential 
current sense amplifier is 50. Then for obtaining the L etua we should 
divide the actual voltage across the shunt resistor by the shunt resis- 
tor value, which is 10 mO as you can see in the following equations: 


(4096) 


50 
Las = RawAdcVal x ————- (3.12) 
0.01 
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Las = RawAdcValx0.001611328125 (3.13) 


Then we can run the program and debug the code. So, we run 
the code in the debug mode to verify that the V etua and I,4,,; are 
really updated with the output voltage and output current. That is 
very important before going ahead to control the converter and 
first thing is to check that the duty cycle can be changed and the 
second part is to check the timer interrupt is being triggered and 
it is actually very important to execute the control code and the 
third part is ensuring that you can read the actual output voltage 
and current with the equations have been proposed because they 
are very important to obtaining closed-loop control. If one of 
them is not set, nothing will be worked basically. We have checked 
the PWM signal before but the others have not been checked yet 
and they will be demonstrated in this section. Therefore, for run- 
ning the program in debug mode, we do right click on the project 
name, select Debug As and then choose Debug Configurations... 
item. In the Debug Configurations window and in the Debugger 
tab, we set the Enable live expressions and then click on the Close 
button. 

Then, we click on the Debug icon (green insect). In the opened 
Confirm Perspective Switch window, we do click on the Switch 
button. From Windows tab in menu, we select Show View and the 
Expressions item. Then we click on the Add new expression, to see 
the expressions we want, V,.,, is one of them and I etua is being 
the other one and we want to show up them, so, we click on the 
Resume (F8) icon in menu. The duty cycle of the PWM signal has 
been set to as below equation: 


PWM duty cycle = m, 100 « 47% (3.14) 


The system is turned on and the values of V, and 1,4, are 
displayed as Table 3.1. 
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TABLE 3.1 Displayed V, and Lctuaı Values 
Expression Type Value 

actual float 10.8957205 
I float 1.44052732 


actual 


We see that values are not stable and changing because it is an 
open-loop control system and we have not added the controller 
codes to make a closed-loop controller yet, so, that is because it is 
uncompensated or open-loop control system. In this case, the 
output resistor is about 7.5 Q. Now, we connect the output voltage 
to the oscilloscope to display the output voltage on it, we can see 
that the value on the oscilloscope screen is different from the 
value displayed on live expressions in debug mode running. 
Because the oscilloscope used in this case has not been calibrated 
that well, so, there is a slight difference with what it senses so it 
reads 12.4 volts but in actuality when you use a multimeter, you 
will see more accurate result of 11.48. We did manual initializa- 
tion of our code and used the live expressions in debug mode run- 
ning the program. In the next section, we will implement the 
voltage and current closed-loop controller and will add them to 
our code. 


3.4 FIRMWARE IMPLEMENTATION OF 
VOLTAGE CONTROL 


In this section, we will discuss about the next part of firmware 
which is implementing the closed-loop control system, especially 
the function that will implement the closed-loop control system, 
that is, changing the duty cycle based on the center inputs. We 
will think about the broader scheme which in terms of codes 
in it and code the firmware in such a way that be easily used in 
other families or use a function or separate file, that you actually 
call to implement the controller, so that it can be used in other 
microcontrollers as well and you basically copy and paste certain 
file and change the duty cycle update and also do ADC update 
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from the DMA. The first step is that basically, you should create a 
header file, so we do right click on the Inc folder, select New, and 
then choose the Header File option. We do type a name, custom.h, 
in Header file row and click on the Finish button. We define a 
void function inside custom.h file as below, since it will not return 
anything and save it. 


#ifndef INC CUSTOM H_ 
#define INC CUSTOM H_ 
void CNTRL ROUTINE(void); 
#endif 


After that, we do right click on the Src folder, select New, and 
then choose the Source File option. We do type a name, custom.c, 
in Source file row and click on the Finish button. First, we are 
going to make the function we defined before, so we add the fol- 
lowing code in custom.c file. 


void CNTRL ROUTINE (void) { 


} 


A lot of extern variables need to be defined inside that function, 
so, we are going to add them from main.c file, in the Private vari- 
ables (PV) section, we add two variables, Vand Iep then you will 
need to set up arrays to store the variables for the pole zero com- 
pensator or digital controller so we define voltage and current 
arrays for errors, e_v and e_i, also we define voltage and current 
arrays for outputs, y_v and y_i, and initialize them to zero as 
below. 


/* USER CODE BEGIN PV */ 
float Vref = 16.0; 

float Iref = 2.133; 
float Vactual = 0.0; 
float Iactual = 0.0; 
float e v[5] = {0}; 
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float y v[5] (0); 
float e i[5] (0); 
float y i[5] = {0}; 
uint32 t adcValues[2]; 
/* USER CODE END PV */ 


Remembering from Equation (2.30), we have substituted e v 
instead of e, ande m instead of y, so e, = e v[0], e, , = e v[1], etc., 
and also y, = y. v[0], y, ; = y. v[1], etc., and y, is actually the duty 
cycle update. The closed-loop block diagram of the system is illus- 
trated in Figure 3.4. 

Then, we go to the custom.c file and add them as extern vari- 
ables as given in the following code. 


include <stdio.h> 
#include <stdint.h> 
#include <stdlib.h> 
#include main.h 

#include custom.h 

extern float Vref; 

extern float Iref; 

extern float Vactual; 
extern float Iactual; 
extern float e v[5]; 
extern float y v[5]; 
extern float e i[5]; 
extern float y i[5]; 
extern uint32 t adcValues[2]; 
void CNTRL ROUTINE (void) { 


} 


Vr,/Ir, + 


——. 


Buck Vo,,/Io,, 


converter 


Controller 


FIGURE 3.4 The closed-loop block diagram of the system. 
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Then, we should add necessary libraries to the custom.c file as 
you can see in the above code; however, including math.h library 
is not necessary in this case. So, we are going inside the CNTRL_ 
ROUTINE function. From Equation (2.30), we should note that 
the y,(y[0]) is used to update the duty cycle and y, ,(y[1]) is the 
history 1, y, ;(y[2]) is the history 2, and y, 4(y[3]) is the history 3. 
After implementing them, you need to update each one, so y, , 
will be equated to what y, is, y, ; will be equated to what y, , is 
and y, , will be equated to what y, , is as below: 


Yn-1 > Yn (3.15) 
Yn-2 > Yna (3.16) 
Yn-3 — Yn-2 (3.17) 


The same updates are applied for the errors, so, first, we do 
copy and paste the V,ctua equation code from interrupt duty func- 
tion inside the CNTRL ROUTINE function and apply the 
updates for outputs and errors as depicted in the following code. 


void CNTRL ROUTINE (void) { 
Vactual = adcValues[0]*0.008137207031; 


y v[3] = y vI21]; 

y v[2] = y vl11]; 

y_v[1] = y vl0]; 

e v[3] = e v[2]; 

e v[2] = e v[11; 

e v[1] = e v[01; 

e v[0] = Vref - Vactual; 

y v[0] = A1 v*y v[1] + A2 v*y v[2] + 


A3 v*y v[3] + BO v*e v[0] + B1 v*e v[1] + 
B2 v*e v[2] + B3 v*e v[3]; 
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if(y v[0] > y v max)( 


y v[0] = y v max; 
} 
if(y_v[0] < y_v_min) { 
y v[0] = y v min; 
} 
TIM1->CCR1 = (int) (MAX DUTY*y v[0]); 


Then we add the codes for e_v[0] and PI controller and define 
limits for maximum and minimum duty cycles, y_v[0]. Since the 
y. v[0] is of float type and is between 0 and 1, we should convert it 
to integer and we know that the MAX DUTY is 849 in this case as 
you can see in the above code. It is the great advantage of STM32 
microcontrollers that you can implement a controller mathemati- 
cally and allow you to define float-type numbers and it usually 
does the mathematical operations in less than 10 us. Now we are 
going to define coefficients and variables that have been utilized in 
the code inside the custom.c file as the following code. 


#define Al v 1 
#define A2 v 0 
#define A3 v 0 
#define BO v 0.001107 
#define Bl v 0.002056 
#define B2 v 0 
define B3 v 0 
#define y v max 0.75 
define y v min 0.0 
#define MAX DUTY 849 


In terms of protecting your microcontroller, one thing should be 
considered and it is the constant duty cycle in Equation (2.1), so we 
should note that when setting the maximum value of duty cycle. 
Therefore, we can choose y. v max to be a little bit more than D value, 
that is, 0.75 as you can see in the above code. So, it will not go above 
75% duty cycle. Next, we should call the function CNTRL_ROUTINE, 
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so we go back to the main.c file and call it inside the interrupt duty 
function and do comment the Voua and L etua expressions, as depicted 


in the following code, since they will be calculated inside the 
CNTRL ROUTINE function in custom.c file. 


/* USER CODE BEGIN 4 */ 
void HAL TIM PeriodElapsedCallback(TIM 
HandleTypeDef* htim)( 
if(htim-»Instance == TIM4) { 
CNTRL ROUTINE () ; 
// Vactual = adcValues [0] 
*0.008137207031; 
// Tactual = adcValues [1] 
*0.001611328125; 
} 
} 
/* USER CODE END 4 */ 


Before starting to build the program, we should remove the fol- 
lowing code inside the loop while (1). 


TIM1->CCR1=400; 


Then, we do build the program by clicking on the Build All 
(Ctrl+B) icon in the menu. The program is successfully built and we 
can now run it by clicking on the run icon in the menu. So, the code 
is successfully uploaded to the STM32 microcontroller, now we can 
see the output voltage using the oscilloscope and multimeter. As 
you can see on the oscilloscope, the voltage displayed is 17.9 V; 
however, on the multimeter screen, the voltage value is 16.36 V. You 
should note that the V,,, has been set to 16 V. In the next section, we 
will evaluate the system operation for different set points. 


3.5 FIRMWARE IMPLEMENTATION OF 

CURRENT CONTROL 
In this section, we will discuss about implementation of current 
closed-loop control specifically output current to duty cycle which 
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will be used to obtain the coefficients of the closed-loop control- 
ler and used for controlling the output current of buck converter. 
Before we actually start, we are going to change the setpoint to 
get the goal of robustness, so we change the set point V, to 10 V 
and see the results to ensure the operation of implemented voltage 
closed-loop control. Therefore, in the main.c file, we modify the 
code as below, and then run the program. 


/* USER CODE BEGIN PV */ 
float Vref - 10.0; 
float Iref - 2.133; 


As you can see, the voltage value on the oscilloscope is 11.5 V; 
however, the output voltage on the multimeter screen is 10.41 V 
which is much less than 11.5 V. As you can see, there is a slight off- 
set from what we actually set up. In essence of how to deal with this 
is basically to increase the average error, the average voltage error is 
to use a more precise voltage sensor that will also affect the loop 
and you will incorporate it with the loop design. Another thing you 
can do is that you simply adjust the reference voltage, V,.,, and make 
it slightly lower and cause the offset value to change. In many appli- 
cations that are actually utilized, you do want to be in a certain 
range, for example, if you set the output voltage to be 3.3 V, and you 
want to get about 3.3 V-3.5 V and you want it to be in certain range. 
You should know there are many factors to cause the output voltage 
to be exactly 10 V. Now, we are going to implement the output cur- 
rent closed-loop control, so, first of all in the custom.h file, we add 
the new function, CNTRL ROUTINE CURRENT, and change 
the name of the voltage control function, CNTRL ROUTINE, to 
CNTRL ROUTINE VOLTAGE as shown in the following code. 


#ifndef INC CUSTOM H 

&define INC CUSTOM H 

void CNTRL ROUTINE VOLTAGE (void); 
void CNTRL ROUTINE CURRENT (void); 
#endif 
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Then we are going to the file custom.c and we create the func- 
tion CNTRL ROUTINE CURRENT, exactly like the function 
CNTRL, ROUTINE VOLTAGE, and we use the I etua code line 
instead of V etua in this case, and also change any v, in arrays or 
variables to i, as you can see in the code below. 


void CNTRL ROUTINE CURRENT (void) { 
Iactual = adcValues[1]*0.001611328125; 


yeaa} x 

y il2] = y il1]; 

y il1] - i[01]; 

e i[3] i 

e i[2] =e i[1]; 

e i[1] = e i[01]; 

e i[0] = Iref - Iactual; 

y i[0] = A1 i*y i[1] + A2 i*y i[2] + 


A3 i*y i[3] + BO ite i[0] + B1 i*e i[1] + 
B2 i*e i[2] + B3 i*e i[3]; 
if(y i[0] > y i max)( 


y i[0] - y i max; 
} 
if(y i[0] < y i min)( 
y i[0] - y i min; 
} 


TIM1->CCR1 = (int) (MAX DUTY*y i[0]); 


We know that the output resistor value is 7.5 Q and we set the 
Lau t0 2.133 A, so, the output voltage will be 16 V. We have used 
596 resistors for resistor divider sensors to get more accurate out- 
put voltage it is better to use 196 resistors to achieve more accurate 
values. Then we are going to define the PI current closed-loop 
controller coefficients as well as maximum and minimum values 
for the duty cycles as depicted in the following code. 


define Al v 
#define A2 v 0 
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#define A3 v 0 
#define BO v 0.001107 
#define Bl v 0.002056 
#define B2 v 0 
#define B3 v 0 
define Al i 1 
#define A2 i 0 
define A3 i 0 
#define BO i 0.008334 
#define Bl i 0.01548 
#define B2 i 0 
#define B3 i 0 
#define y i max 0.75 
#define y i min 0.0 
#define y_v_max 0.75 
#define y_v_min 0.0 
#define MAX DUTY 849 


Then we go to the interrupt duty function and call the current 
controller function inside of it as given in the code below. 


/* USER CODE BEGIN 4 */ 
void HAL TIM PeriodElapsedCallback (TIM_ 
HandleTypeDef* htim) { 
if(htim->Instance == TIM4) { 
// CNTRL ROUTINE VOLTAGE () ; 
CNTRL ROUTINE CURRENT () ; 

} 
} 
/* USER CODE END 4 */ 


You should note that we do not call the voltage controller func- 
tion here, so we do comment it. Then we can build the program by 
clicking on the Build All (Ctrl+B) icon in the menu and it is built 
successfully. Since, the Z etua is 2.133 A, and the output resistor is 
7.5 Q, we should expect to have the output voltage in the range 
from 16 V to 16.5 V. So, we can run the program to upload it to the 
STM32 microcontroller and the download is verified successfully. 
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Now, we can see the output voltage value on the oscilloscope 
which is 18.1 V; however, it is 16.51 V on the multimeter screen. 
That is all about the current closed-loop controller, and basically 
all of the works have been done. We want to have precise current 
closed-loop controller, so we used a shunt in this case. However, 
when you are going to sense high currents above 15 A, it is not 
advisable to use a shunt because the shunt color is changing and 
also its resistance is changing at some point and affects the effi- 
ciency of the system. In certain cases, you want to utilize Hall 
effect current sensors and they are much better and they are 
highly efficient as well and also its resistance is almost negligible. 
So, in that case, Equation (3.13) for Leta will not suffice. How do 
we modify the equation in case when we use the Hall effect sen- 
sor? The Hall effect sensor essentially works with offset. It is usu- 
ally an offset value and sensitivity. For example, the ACS712 
current sensor is a Hall effect current sensor and it has an offset 
value because it has been designed to work with DC and AC sys- 
tems and we should scale the raw ADC value with the offset value 
and we should subtract offset from it however we prefer to sub- 
tract it from offset value since we will not pass 2.5 V in this case 
since we assume the current is flowing in reverse direction, that is, 
from IP- to IP+ in ACS712 (-20 A to 20 A). So, we have the follow- 
ing equations: 


ortserval -[ Renier Vaa ) 
Lactual = E (3.18) 
Sensitivity 
NE OffsetVal RawAdcVal x Vef (3.19) 


Sensitivity Sensitivity x 24¢** 


OffsetVal = 2.5 V (3.20) 
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Sensitivity =0.1V/A (3.21) 
Vie =3.3 V (3.22) 
2^4mes — 212 — 4096 (3.23) 
Lau = 25 - RawAdcVal x 0.008056640625 (3.24) 


So, we can use it in the function CNTRL_ROUTINE_ 
CURRENT in custom.c file as given in the following code. 


void CNTRL ROUTINE CURRENT (void) { 
// Tactual = adcValues [1] *0.001611328125; 
Iactual = 25-adcValues [1] *0.008056640625; 


You should note that the Hall effect current sensors are very 
prone to be affected by magnetic fields. In the next section, we will 
discuss about average current mode control. 


3.6 AVERAGE CURRENT MODE CONTROL OF 
THE BUCK CONVERTER 


In this section, we will deal with average current mode control 
and essentially what it means you actually control the current 
through the inductor, so instead of controlling the output volt- 
age or output current, you basically control the current flowing 
through the inductor and what that tells us, you have double 
loop also it is called double loop control, and what you get in 
double loop control is that the outer voltage loop makes a ref- 
erence for the inner current loop and it really presents many 
advantages in the converter control because you can essentially 
ensure that your converter constantly operates in continuous 
conduction mode, CCM, and additionally you can actually save 
your MOSFET from a lot of stress because inherently you are 
limited in the amount of the current goes through your inductor 
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by proxy MOSFET or diode as well, essentially the amount of 
current goes through the converter so you basically limit it in 
every cycle as opposed you have not checked it in other meth- 
ods, that is, the output voltage to duty cycle or output current 
to duty cycle where you control them like a single loop control. 
If you want to limit current, you essentially have to use other 
than the duty cycle method, we keep our fingers crossed, or, lit- 
erally, you have to use another if statement, that is, if the cur- 
rent above the limit reduces the duty cycle as opposed to double 
control loop inherently and you just limit the current by default 
as you are using the ADC of course. In this section, we will dive 
into practical implementation of this method and we will not 
go through the theoretical analysis because the average current 
mode control in theoretical is rather complicated and we will 
go to other approach so we will implement the average current 
mode in the firmware and we are going to start coding it. You 
should note that the outer control loop is quite slow, so, every 
time the outer control loop is executed, the current control loop 
is executed 10 times, so the voltage control loop only executes 1 
iteration for 10 iterations of current loop control. Therefore, in 
main.c file and Private variables section, we define a loop coun- 
ter variable (loop cnt) and initialize it to zero as given in the 
code below. 


/* USER CODE BEGIN PV */ 


float Vref - 10.0; 
float Iref - 1.2; 
float Vactual = 0.0; 
float Iactual = 0.0; 
float e v[5] = {0}; 
float y v[5] = {0}; 
float e i[5] = {0}; 
float y i[5] = {0}; 
int loop cnt - 0; 


uint32 t adcValues [2]; 
/* USER CODE END PV */ 
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Then, we go to the custom.h file and define the functions we 
require in this case as you can see in the following code. 


#ifndef INC CUSTOM H 

#define INC CUSTOM H_ 

void CNTRL ROUTINE VOLTAGE (void); 
void CNTRL ROUTINE CURRENT (void); 
void DOUB LOOP (void); 

void DOUB LOOP VOLTAGE (void); 
void DOUB LOOP CURRENT (void); 
#endif 


You should note that the first function is for the entire double 
loop controller, the second one is for the voltage control loop, and 
the third one is for the current control loop. 

In custom.c file, we define the loop counter variable (loop_cnt) 
as external integer and add the codes inside the entire double loop 
control function and call the double loop current control in every 
iteration and double loop voltage control in every 10 iterations as 
given in the following code. 


#include <stdio.h> 
#include «stdint.h» 
#include <stdlib.h> 
#include main.h 
#include custom.h 


#define Al_v 1 
#define A2 v 0 
#define A3 v 0 
#define BO v 0.001107 
#define Bl v 0.002056 
#define B2 v 0 
#define B3 v 0 
#define Al i 1 
#define A2 i 0 
define A3 i 0 
#define BO i 0.008334 
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#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 
extern 


extern 


igital Control of Power Converters 
Bl i 0.01548 
B2 i 0 
B3 i 0 
y i max 0.75 
y i min 0.0 
y v max 0.75 
y v min 0.0 
MAX DUTY 849 

float Vref; 

float Iref; 

float Vactual; 

float Iactual; 

float e v[5]; 

float y v[5]; 

float e i[5]; 

float y i[5]; 

uint32 t adcValues[2]; 


int loop cnt; 


void CNTRL ROUTINE VOLTAGE (void) { 
Vactual = adcValues[0]*0.008137207031; 


= y vl2]; 
-y vill; 

- y vl0]; 

= e v[2]; 

- e v[1]; 

- e v[0]; 

- Vref - Vactual; 

= Al v*y v[1] + A2 v*y v[2] + 

y v[3] + BO v*e v[0] + B1 v*e v[1] 


B2 v*e v[2] + B3 v*e v[3]; 
if(y v[0] > y v max)( 
y v[0] = y v max; 
} 
else if(y_v[0] < y_v_min) { 
y v[0] = y v min; 
} 
TIM1->CCR1 = (int) (MAX DUTY*y v[0])j; 
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void CNTRL ROUTINE CURRENT (void) { 
Iactual = adcValues[1]*0.001611328125; 
y il[3] = y il21]; 
y_il[2] = yvy. iul; 
y i[1] = i[0]; 
e i[3 Ego. 
e i[2] =e il[1]; 
e i[1] = e il[0]; 
e i[0] = Iref - Iactual; 
y i[0] = A1 i*y i[1] + A2 i*y i[2] + 


A3 i*y i[3] + BO ite i[0] + B1 i*e i[1] + 
B2 i*e i[2] « B3 i*e i[3]; 
if(y i[0] > y i max)( 


y i[0] = y i max; 
} 
else if(y i[0] < y_i_min) { 
y i[0] = y i min; 
} 
TIM1->CCR1 = (int) (MAX DUTY*y i[0]); 


} 


void DOUB LOOP (void) { 
loop_cnt++; 
if (loop ent == 10) { 
loop_cnt = 0; 

DOUB LOOP VOLTAGE () ; 


} 


DOUB_LOOP CURRENT () ; 


} 


void DOUB LOOP VOLTAGE (void); 
void DOUB LOOP CURRENT (void); 


You should note that the double loop current control function, 
DOUB LOOP CURRENT, is identical to the current control 
routine function, CNTRL. ROUTINE CURRENT, however, the 
double loop voltage control function, DOUB. LOOP VOLTAGE, 
is slightly different from the voltage control routine function, 
CNTRL ROUTINE VOLTAGE. We do copy and paste the codes 
inside the CNTRL. ROUTINE VOLTAGE function to inside the 
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DOUB LOOP VOLTAGE function and remove the last code line 
of it since the output of this function is the reference current, Isp 
and it will be used asa reference in the DOUB. LOOP CURRENT 
function, so we add the reference current code as given in the code 
below. 


void DOUB_LOOP_VOLTAGE (void) { 
Vactual = adcValues [0] *0.008137207031; 


y v[3] = y v[2]; 

y v[2] = y v[1]; 

y ví[1] = y vl01]; 

e v[3] = e vI21; 

e v[2] = e v[11; 

e v[1] = e v[01; 

e v[0] = Vref - Vactual; 

y v[0] = A1 v*y v[1] + A2 v*y v[2] + 


A3 v*y v[3] + BO v*e v[0] + B1 v*e v[1] + 
B2 v*e v[2] + B3 v*e v[3]; 
if(y v[0] > y v max)( 


y v[0] = y v max; 

} 

else if(y v[0] « y_v_min) { 
y v[0] = y v min; 

} 

Iref = Iref max*y v[0]; 


} 


void DOUB LOOP CURRENT (void); 


The Iref max is the maximum current which is essentially 
required and that is the maximum inductor current you want, 
and we never get the inductor current above that value. Then we 
go back to the custom.c file and define Iref max as 2.133 and 
change y. v max from 0.75 to 1.0 as you can see in the following 
code. 


//#define y v max 0.75 
#define y_v_max 1.0 
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define y v min 0.0 
#define Iref max 2.133 


Then we do copy and paste the codes inside the CNTRL_ 
ROUTINE CURRENT function to inside the DOUB LOOP 
CURRENT function and you have used Iç which has been 
obtained directly from the double loop voltage control function, 
DOUB LOOP VOLTAGE, as given in the following code. 


void DOUB LOOP CURRENT (void) { 
Iactual = adcValues[1]*0.001611328125; 


yil3]-2y. as l; 

y_il2] = Lg 

y_i[l] = y_ mo 

e i[3] =e i[2]; 

e i[2] =e il[1]; 

e i[1] =e i[0]; 

e i[0] = Iref - Iactual; 

y i[0] = A1 i*y i[1] + A2 i*y i[2] + 


A3 i*y i[3] + BO ite i[0] + B1 i*e i[1] + 
B2 i*e i[2] « B3 i*e i[3]; 
if(y i[0] > y i max)( 


y i[0] = y i max; 
} 
else if(y i[0] < y_i_min) { 
y i[0] = y i min; 
} 
TIM1->CCR1 = (int) (MAX DUTY*y i[0]); 


So, the reference current, I „p basically keeps the inductor cur- 
rent in CCM, which is literally constant, remember that I, is run- 
ning for 10 cycles and the function DOUB_LOOP_CURRENT 
runs for 10 times at first before start up the function DOUB_ 
LOOP VOLTAGE, so you should note that adequate Ie should be 
there and not be too high that is why we initiate it to 0.2 A first, 
until it gets updated inside the function DOUB_LOOP_VOLTAGE, 
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so we do the initialization of I „p inside the main.c file as given in 
the code below. 


/* USER CODE BEGIN PV */ 
float Vref 10.0; 
float Iref 02:5 


Also, as we mentioned earlier, the function we want to call is 
the DOUB LOOP function in every 100 us inside the interrupt 
duty function. Therefore, we go to the main.c file and find the 
interrupt duty function and call the DOUB LOOP function 
inside of it as you can see in the following code. 


/* USER CODE BEGIN 4 */ 
void HAL TIM PeriodElapsedCallback(TIM 
HandleTypeDef* htim)( 
if(htim-»Instance == TIM4) { 

DOUB LOOP(); 

// CNTRL ROUTINE VOLTAGE () ; 

// CNTRL ROUTINE CURRENT () ; 

} 
} 
/* USER CODE END 4 */ 


Then, we can build the program and run the program to upload 
the code to the STM32 microcontroller. Just keep in mind that the 
system is not modifying the duty cycle directly as a result of the 
sensed output voltage, it takes and generates a reference current 
value from the output voltage and the reference current value is 
what used to actually modify the duty cycle directly, so it is almost 
indirect control instead of direct control. As you can see on the 
oscilloscope, the output voltage is not constant and it is going up 
and down from 11.3 V to 11.5 V, so you need to modify the voltage 
and current controller coefficients. So, we go back to the code and 
inside the main.c file we change the set point, V, value, from 10 V 
to 15 V. You should note that we can only change the set point Ve 
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and we cannot change the Ip manually and it is changed auto- 
matically and dynamically based on what Vep is, so we run the 
program gain. There is some offset in the oscilloscope and it is 
changing from 16.7 V to 16.9 V. So, that is basically what the aver- 
age current mode control is and it is very useful in a vast majority 
of applications such as battery chargers and, of course, you cannot 
use average current mode control if you are not operating in con- 
tinuous conduction mode, CCM. If the inductors are not large 
enough, they will not be operating in CCM most time. We should 
note that we are getting close to the rated, that is, maximum refer- 
ence current, Iref max, which is 2.133 A in this case. The last set- 
point, V. that we want to use is 16 V, so we go to main.c file and 
set it as given in the code below. 


/* USER CODE BEGIN PV */ 
float Vref - 16.0; 
float Iref 0. 2 


Also, we modify slightly Iref max from 2.133 to 2.2 in custom.c 
file as you can see in the following code. 


//#define y_v_max 0.75 
#define y_v_max 1.0 
#define y_v_min 0.0 
#define Iref max 2.2 


Then, we run the program again, and we see on the oscillo- 
scope that the output voltage goes up and down from 179 V to 
18.1 V. That is the implementation of average current mode con- 
trol. The noise in the output voltage is mostly because of the AC/ 
DC converter, using the output of linear voltage regulators as an 
input to the buck converter will give better results. It is very 
advantageous to control the current through the inductor and 
average current mode control can do that very well. If we modify 
the coefficients, we can reduce the output voltage oscillations, for 
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example, we change the BO v and B1 v and also BO i and Bl i 
values to one-tenth of their previous values in custom.c file as 
given in the code below. 


#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 


Al v 
A2 v 
A3 v 
BO v 
B1 v 
B2 v 
B3 v 
Alii 
A2 i 
A3 i 
BO i 
B1 i 
B2 i 
B3 i 


OOo 000n000000nHidn 


By doing that, you actually increase 
coefficients are more than what actually they are, and when you 
reduce them, you actually reduce the settling time and you will 
have more stable output. Then we run the program, we can see 
that the output voltage is now more stable and it is not going up 
and down so much. We again change more the BO. i and Bl. i val- 
ues to one-tenth of their previous values as you can see in the fol- 


lowing code. 


#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 
#define 


Al v 
A2 v 
A3 v 
BO v 
B1 v 
B2 v 
B3 v 
Alii 
A2 i 


O n|ooo0o0o0o0o0nHn 


.0001107 
.0002056 


.0008334 
.001548 


the settling time when 


.0001107 
.0002056 
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#define A3 i 0 
#define BO i 0.00008334 
#define B1 i 0.0001548 
#define B2 i 0 
#define B3 i 0 


Then, we run the program and can see that the output voltage 
is now more stable and it is fixed at 18.1 V. 


CHAPTER 4 


Digital Control 
Implementation with 
PLECS 


4.1 INTRODUCTION 


PLECS is a simulation software tool used for the design and anal- 
ysis of power electronic systems. It stands for piecewise linear 
electrical circuit simulation and is developed by Plexim GmbH. 
PLECS provides a user-friendly interface for modeling complex 
power electronic systems and simulating their behavior in real- 
time. Key features of PLECS include: 


e Intuitive graphical modeling environment for designing 
power electronic circuits. 


e Support for modeling various power converters, inverters, 
motor drives, and other power electronic systems. 


e Real-time simulation capabilities for rapid prototyping and 
testing of control algorithms. 
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e Integration with MATLAB/Simulink for system-level simu- 
lation and control design. 


* Code generation capabilities for automatically generating 
embedded code for microcontrollers. 


e Extensive library of pre-built power electronic components 
and control blocks. 


e Validation and verification tools for analyzing system behav- 
ior and performance. 


Overall, PLECS is a powerful tool for engineers and research- 
ers working in the field of power electronics, enabling them 
to model, simulate, and optimize the performance of power 
electronic systems efficiently. You can generate code for power 
converter control in PLECS for STM32 microcontrollers and 
implement your control algorithm on the microcontroller to 
regulate the power converter operation. In this chapter, we will 
investigate and study a practical example of implementing a 
power converter, designing a controller, and generating code 
for the STM32 Nucleo-G474RE microcontroller in detail using 
PLECS and MATLAB software. 


4.2 DESIGNING A CLOSED-LOOP 

CONTROLLER IN PLECS AND MATLAB 
In this section, we intend to design a buck converter in the PLECS 
software environment and simulate it. The designed buck con- 
verter is shown in Figure 4.1. 

The circuit parameter values are clearly visible in Figure 4.1. 
The switching frequency is 20 kHz and the triangular wave gen- 
erator amplitude changes between 0 and 1, and it is compared 
with constant duty cycle value with a logical comparator or rela- 
tional operator to generate the PWM signal and, finally, it is fed to 


Vin 
Vv: 12 


[5] 


Constant 


Relational 
Operator 


Triangular Wave 
Min: 0 
Max: 1 
f: 20e3 
DutyCyde: 0.5 
Delay: 0 


FIGURE 4.1 The designed buck converter in PLECS. 
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File Edit View Window Help 
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FIGURE 4.2 The output voltage waveform. 


the gate driver of the MOSFET. Also, there is a load change at 100 
ms and the load resistance is changing from 8 © to 4 Q at the time 
of 100 ms. The output voltage waveform is also depicted in Figure 
4.2 which has an undesirable overshoot at the startup and there is 
an offset or a deviation from the reference voltage value, 5 V, as 
shown in Figure 4.2. The simulation parameters window for this 
system is illustrated in Figure 4.3. 

For designing a PID controller for this power converter, we 
need to know the control to output frequency response of the 
system. Therefore, in library browser, control, and small signal 
analysis section, we do drag and drop the small signal pertur- 
bation and add it with the reference voltage value and also do 
drag and drop the small signal response as shown in Figure 4.4 
and connect it to the measured output voltage as shown in 
Figure 4.1. 
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Solver | Options | Diagnostics ] Initialization | 
Simulation time 


Start time: [0.0 | Time span: [0.2 


Solver 


Type: | Variable-step v | Solver: 


Solver options 


Max step size: [ie-6 B ] Relative tolerance: 


Initial step size: (auto | Absolute tolerance: 


Refine factor: |1 | 


Circuit model options 


Diode turn-on threshold: [o 


FIGURE 4.3 The simulation parameters window. 


To see the control to output frequency response of the system 
from Simulation tab in menu, we select the Analysis too as shown 
in Figure 4.5. 

In the opened Analysis Tools window, we do click on the plus 
(+) button and select the AC Sweep Analyses and do the settings 
as shown in Figure 4.6. 

To see the control to output frequency response, that is, Bode 
diagram of magnitude and phase of the system we do click on the 
Start analysis button in Analysis Tools window which is depicted 
in Figure 4.7. 

We will use MATLAB software to design a PID controller for 
closed-loop system so we need to export the frequency response 
data as a CSV file as illustrated in Figure 4.8. 
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Pl X 
File Window Help 


> System 
> Assertions 
4 Control 
> Sources 
> Math 
> Continuous 
> Delays 
> Discontinuous 
> Discrete 
> Filters 
> Functions & Tables 
> Logical 
> Transformations 
> Modulators 
> State Machine 
4 Small Signal Analysis 


© Small Signal Perturbation 


[~] Small Signal Response 


© Small Signal Gain 


> Electrical 

> Thermal 

> Magnetic 
> Mechanical 


FIGURE 4.4 Choosing the small signal analysis elements. 
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a! 


File Edit View | Simulation | Format Coder Window Help 


Start Ctrl+T 
Pause Space 
Simulation parameters... Ctrl+E 
Analysis tools... Ctrl- Shift- Y 
Simulation scripts... Ctrl+Shift+T 


FIGURE 4.5 Selecting the analysis tools... from simulation tab. 


Analysis type: AC Sweep 
Description: [AC Sweep 


Setup | Options | Steady-State Options | 


Operating point: | periodic 
System period: 5e-5 


Frequency range: |[10 10000] 


Amplitude: 1e-3 


Perturbation: 


Response: 


Accept Revert 


FIGURE 4.6 Settings of the analysis tools window. 


In the opened Export as window, we choose a name (data.csv) 
for it and click on the Save button to save it in the desired direc- 
tory as shown in Figure 4.9. 

Then we open the saved CSV file and delete the header or first 
row of it and save it again as depicted in Figure 4.10, since if we do 
not do that when we read the CSV file in MATLAB command 
window, we will get an error message. 
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| File Edit View Window Help 
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FIGURE 4.7 The control to output frequency response. 


Therefore, we use the following codes in the MATLAB com- 
mand window to read the saved data.csv file and plot the Bode 
diagram of the system in MATLAB. 


>> data = csvread('C:\plecs\stm32\data.csv’); 
>> w = 2*pi*data(:,1); 

>> val-10.^(data(:,2)/20).*exp(j*data(:,3) 
*pi/180) ; 

>> sys = frd(val,w); 

>> bode(sys),grid minor 


By running those codes in the MATLAB command window, 
the Bode diagram of the system is plotted in MATLAB as depicted 
in Figure 4.11. 
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FIGURE 4.9 Choosing a name and saving the CSV file. 


Digital Control Implementation with PLECS m 159 


File Edit Format View Help 


10, -0.73128178326590443, -0.009437196356044621 
-481056257175661, -0.73126256690536962, -0.010834883024415261 
.183915622940013, -0.73123720987151342, -0.01244204531966551 
-140045420136261 , -0.731203770672747, -0.014288209678226338 
. 376194613379671, -0.73115986778838571, -0.016398771957387523 


. 960079840319363 , -0.731101625503357, -0.018837703144581257 

. 90950744558992 , -0.73102528823058632, -0.021621836725465221 
- 315789473684209 , -0.73092405883380196, -0.024837444445185257 
. 211480362537763 , -0.73079110746094345, -0.028515555692431464 
- 662045060658578 , -0.73061678833250698, -0.032718152981864042 
. 840637450199203 , -0.73038384143149671, -0.037609191886669564 


Ln 1, Col 1 a 


FIGURE 4.10 Removing the first row of data.csv. 
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FIGURE 4.11 The Bode diagram of the system in MATLAB. 
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To get the transfer function of the system, we use the following 
code in the MATLAB command window. 


>> H = tfest(sys,2) 
H = 

2497 s + 5.186608 

s*2 + 4201 s + 5.639608 
Continuous-time identified transfer function. 
Parameterization: 
Number of poles: 2 Number of zeros: 1 
Number of free coefficients: 4 
Use "tfdata", "getpvec", "getcov" for 
parameters and their uncertainties. 
Status: 
Estimated using TFEST on frequency response data 
"Sys". 
Fit to estimation data: 99.2% (simulation focus) 
FPE: 7.625e-05, MSE: 7.07e-05 


As you can see above, estimated transfer function of the system 
with an accuracy of 99.2% is as below: 


8 
H- k - (4.1) 
s^ +4201s +5.639x10 


Then we use the following code in MATLAB command win- 
dow to design a closed-loop PID controller for the system. 


>> pidTuner (H) 


By running the code above, the PID Tuner window is opened 
and we select the Type as PID and Form as Parallel and also set the 
Response Time (0.1173) and Transient behavior (0.66), so that an 
appropriate step plot is achieved as shown in Figure 4.12. 

Then we click on the Export button in the PID Tuner window 
and in the opened Export Linear System we do enter a name (PID) 
in Export PID controller and click on the OK button as illustrated 
in Figure 4.13. 
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FIGURE 4.12 Setting the PID tuner window. 


I Export PID controller |PID 


Export plant model 


YOK $$ Cancel @) Help 


FIGURE 4.13 Entering a name (PID) in export PID controller window. 


Finally, in the MATLAB command window, we run the code 
below to get the closed-loop PID controller transfer function. 


>> PID 
PID = 
1 
Ki , e 
S 
with Ki = 18.5 
Continuous-time I-only controller. 
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4 Control 
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Discrete PID Controller 


FIGURE 4.14 Finding the continuous PID controller block. 


So, the designed closed-loop PID controller transfer function is 
as follows: 


PID(s) Qi (4.2) 
S 


In the Library Browser, we can find the Continuous PID 
Controller block in Control and Continuous section and drag and 
drop it to the uncompensated system as shown in Figure 4.14. 

Then, we do double click on the Continuous PID Controller 
block and in the Basic tab we enter the PID(s) transfer function 
values as depicted in Figure 4.15. 

Also, in the Anti- Windup tab of the Continuous PID Controller 
window, we set the upper saturation limit to 0.9 and lower satura- 
tion limit to 0.1, as shown in Figure 4.16, to limit the duty cycle 
between 0.1 and 0.9 in order to safe operation of the MOSFET. 
Then we can add the Continuous PID Controller to the closed- 
loop control system as illustrated in Figure 4.17. 

Finally, we simulate the compensated closed-loop control sys- 
tem and we can see the output voltage waveform on the Scope as 
shown in Figure 4.18. 

As you can see in this case the output voltage offset is omitted 
and the output voltage is equal to the reference voltage value 
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FIGURE 4.15 Setting the PID(s) transfer function values. 


Continuous PID Controller (mask) (link) 


Configurable PID-Controller (parallel form) in continuous time-domain with filtered 
derivative term and two different anti reset-windup methods. 


Assertions 


FIGURE 4.16 Setting the upper and lower saturation limits in the anti- 
windup tab. 
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FIGURE 4.17 Adding the PID(s) block to the closed-loop control system. 
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FIGURE 4.18 The output voltage waveform on the scope. 


which is 5 V. Also, the startup overshoot in output voltage is 
removed and following a step load resistance change from 8 Q to 
4 Q at 100 ms, the output voltage becomes stable in less than 5 ms. 


4.3 CODE GENERATION FOR THE STM32 
NUCLEO-G474RE IN PLECS 

Code generation for the STM32 Nucleo-G474RE in PLECS 

involves converting the simulation model created in PLECS into 

executable code that can be run on the specific microcontroller. 

Here is a general overview of the steps involved in code generation 

for the STM32 Nucleo-G474RE in PLECS: 


1. Design the control system in PLECS: Create and simulate 
the control system model in PLECS to ensure that it meets 
the desired performance specifications. 


2. Configure the target settings: Select the STM32 Nucleo- 
G474RE microcontroller as the target for code generation in 
PLECS by specifying the appropriate settings in the Target 
Setup window. 
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3. Generate code: Click on the Generate Code button in PLECS 
to initiate the code generation process. PLECS will convert 
the simulation model into C code that can be compiled and 
uploaded to the STM32 Nucleo-G474RE. 


4. Compile and upload code: Compile the generated C code 
using an Integrated Development Environment (IDE) such 
as the STM32CubeIDE. Upload the compiled code to the 
STM32 Nucleo-G474RE microcontroller using a program- 
ming tool such as ST-Link. 


5. Testthe control system: Run the control system on the STM32 
Nucleo-G474RE to verify its functionality and performance. 
Make any necessary adjustments to the code or simulation 
model in PLECS to optimize the system's performance. 


By following these steps, you can effectively generate code for the 
STM32 Nucleo-G474RE microcontroller in PLECS and imple- 
ment the control system on the microcontroller for real-world 
applications. Now, we are going to create the same buck converter 
using the STM32 PLECS library and then generate the code for 
the STM32 Nucleo-G474RE microcontroller in STM32CubeIDE. 
To do that, we create two subsystems one for the buck converter 
and the other one for the PID controller as shown in Figure 4.19. 

The circuit inside the Buck. Converter subsystem is depicted in 
Figure 4.20. 

As you can see in Figure 4.20, we have added a resistor divider 
with the values R5 = 100 kO and R6 = 10 kO, since we want to get 


FIGURE 4.19 Creating two subsystems for buck converter and PID 
controller. 
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FIGURE 4.20 The circuit inside the Buck Converter subsystem. 
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Duty 
CarrierType: Symmetrical (carrier counts up/down) 
CarrierFreq: 20e3 


FIGURE 4.21 The circuit inside the PID Controller subsystem. 


a voltage below 3.3 V to the input of ADCI pin of the STM32 
microcontroller. Therefore, we have the following equation: 


B6. inso VU (4.3) 


VO sense — mu 
R54 R6 11 


Also, the circuit inside the PID Controller subsystem is illus- 
trated in Figure 4.21. 

As you can see in Figure 4.21, we do multiply the output of 
ADCI by a Gain which is now equal to 11 to get the actual value 
of output voltage, Vo, and then it is compared with the reference 
voltage, 5 V. The ADCI block parameters window is shown in 
Figure 4.22. 

As you can see in Figure 4.22, we could set the Cont. conversion 
scale(s) to 11 instead of adding a Gain with the value of 11. The 
PWM block parameters window is depicted in Figure 4.23. 

Finally, from the Simulation tab in menu, we select the Start 
option as depicted in Figure 4.24, to see the results of simulation 
using the STM32 PLECS library elements and designed controller 
in the previous section as illustrated in Figure 4.25. 

Then we should enable the code generation option for the 
PID_Controller subsystem. Therefore, we do right click on the 
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FIGURE 4.22 The ADCI block parameters window. 


PWM (mask) (ink) 
Generate a complementary PWM signal pair. 


The carrier starts at 0 and varies between 0 and 1. The PWM output is active when the input is greater than the 
carrier. 


Channel i | Channel2 | Channel3 | Channel4 | Sync | Trigger | Protection | 
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Carrier type: Enable port: 
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FIGURE 4.23 The PWM block parameters window. 
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FIGURE 4.25 The results of simulation. 


PID_Controller subsystem, select Subsystem, and then choose the 
Execution settings... option as shown in Figure 4.26. 

In the opened Execution Settings window, we tick on the 
Enable code generation and enter the value of le-4 in Discretization 
step size field and then click on Apply and OK buttons, respec- 
tively, as shown in Figure 4.27. 
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FIGURE 4.26 Choosing the execution settings... option. 
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FIGURE 4.27  Ticking on the enable code generation. 
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File. Edit Source Refactor Navigate Search Project 
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FIGURE 4.28 Selecting open projects from file system... option. 


Then we run the STM32CubeIDE to enter it and from File icon 
we select Open projects from File System... option as depicted in 
Figure 4.28. 

In the Open projects from File System window, we click on the 
Archive... button as illustrated in Figure 4.29, to open the g474. 
zip file in the following directory as shown in Figure 4.30. 


PLECS4.7(64 bit) > tsp stm32 > projects > g474.zip 


Then we do click on the Finish button in the Import Projects 
from File System or Archive window as shown in Figure 4.31. 


Import Projects from File System or Archive 


This wizard analyzes the content of your folder or archive file to find projects and import them in the IDE. 


Import source: || 


[ Directory... || Archive. - 
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C Close newly imported projects upon completion 


Use installed project configurators to: 
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Detect and configure project natures 


Working sets 


C Add project to working sets 


Finish 


FIGURE 4.29 Clicking on the archive... button. 
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FIGURE 4.30 Opening the g474zip file. 
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FIGURE 4.31 Clicking on the finish button. 
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FIGURE 4.32 Selecting the properties option. 


In the Project Explorer, we find the cg folder and then right 
click on it and select the Properties option as depicted in Figure 
4.32. Therefore, the Properties for cg window is opened and we do 
copy the Location directory as illustrated in Figure 4.33. 

Then, we close the properties for cg window and come back to 
the PLECS environment and in the Coder icon of menu we select 
the Coder options... as shown in Figure 4.34, to open the Coder 
Options window and in the Target section we select the Target as 
STM32G4x and chip as G474xx, the Build type as generate code 
into STM32CubeIDE project and also we do paste the copied 
Location directory in Figure 4.33 to the STM32CubeIDE project 
directory and also tick use internal oscillator and finally click on 
the Accept and Build buttons, respectively, as shown in Figure 4.35. 
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FIGURE 4.33 Doing copy the location directory. 
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FIGURE 4.34 Selecting the coder options.... 


As shown in Figure 4.36, we get an error message since we have 
not defined the I/O components for the output ports, m and e v, 
of the PID Controller subsystem. We should note that we used 
both of them for simulation purpose and in code generation we 
can remove them from the PID Controller subsystem as depicted 
in Figures 4.37 and 4.38, respectively. 
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Build type: | Generate code into STM32CubelDE project 
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FIGURE 4.35 Settings of the coder options window. 
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FIGURE 4.36 Getting an error message. 


Again, in the Coder Options window, we do click on the build 
button as illustrated in Figure 4.35. The build is done successfully 
and we come back to the STM32CubeIDE and in the Project 
Explorer section we do right click and select the Refresh option as 
illustrated in Figure 4.39 to see the codes added to the cg folder as 
shown in Figure 4.40. 
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FIGURE 4.37 Removing e v and m ports from the PID Controller 
subsystem. 


FIGURE 4.38 The PID Controller subsystem without e v and m ports. 
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FIGURE 4.39 Selecting the refresh option. 
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FIGURE 4.40 The codes added to the ‘cg’ folder. 
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FIGURE 4.41 Clicking on the build all (Ctrl+B) button. 


[92 Problems J£ Tasks Œ Console X [Z] Properties x|2 oS] Ba - &[gl m a-nes B 
CDT Build Console [cube g474 v1 4] 
Finished building: cube g474 vi 4.1ist ^ 


Copy binary to codegen folder: 
cp "cube g474 vi 4.elf" "C:/plecs/stm32Nbuck loop stm32 codegen/PID Controller.elf" 


09:23:23 Build Finished. 0 errors, 1 warnings. (took 3s.204ms) 


« > 


FIGURE 4.42 The build finished message. 


Then we click on the Build All (Ctrl+B) button to build the 
project as shown in Figure 4.41. 

The build is finished with zero errors as depicted in Figure 4.42. 

Finally, we can run the program to upload the codes to the 
STM32 microcontroller as illustrated in Figure 4.43. 


Digital Control Implementation with PLECS m 179 


File Edit Source Refactor Navigate Search Project Run Window Help 


iy |$-&-Giwidbid 6 -($-G-:1 -[O|- 
Run 
[È Project Explorer X > fg 
ASY 3 


4 > g474.zip_expanded 
s lipe] cube_g474_v1_4 (in g474) 


FIGURE 4.43 Running the program. 


4.4 REVIEW OF SMALL SIGNAL ANALYSIS 
TOOLS FOR CONTROLLER DESIGN IN 
PLECS 


We're going to review how you can use the small signal analysis 
tools found in PLECS for controller design. Since the focus is 
going to be on using the analysis tools, we may not comment on 
every single detail within our PLECS circuits that we're going to 
be using. The electronic engineering software, PLECS, is quite 
widely used and adopted in various industries and in academic 
institutions worldwide and it is a complete power conversion sys- 
tem simulation package. It yields very robust and fast simula- 
tion results and it's a fully integrated power electronic system 
simulation package. PLECS standalone offers an independent 
solution. PLECS features a comprehensive component library, 
which covers not only the electrical aspects but also the mag- 
netic mechanical and thermal aspects of power conversion sys- 
tems and their associated controls. So, let's quickly go over the 
development tools that we'll actually be showing in this section. 
So, we'll be using PLECS standalone to demonstrate the capa- 
bilities of the built-in analysis tools models will contain com- 
ponents from the electrical component library, which includes 
both passive and active electrical components like resistors and 
semiconductors. The magnetic, thermal, and mechanical mod- 
eling domains are also available in PLECS. With the help of the 
built-in analysis tools and the controls library, we will design 
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and model our controller. PLECS does offer other tools that fur- 
ther support the simulation development and testing of power 
electronics controllers. In this section, the application that we're 
going to be using to demonstrate the built-in analysis tools is 
an auxiliary-based power supply based on the flyback converter 
topology as shown in Figure 4.44. 

So, because we will be saving some time by using a few prebuilt 
models, well take a moment here to review what's already been 
developed. As we open new models, well explain the changes 
that we've made. So hopefully it's easy for you to follow the flow. 
As an input to the converter, here in Figure 4.44, we have a single- 
phase rectifier. Were using a simplified high-frequency trans- 
former with three output windings, 5 V and plus and minus 12 V. 
The total output power is around 30 W. This is a very low-cost 
simple implementation. It's a popular choice for power supply 
designs. The topology provides isolation between the primary 
and the secondary sides, and it also then gives the designer the 
ability to provide multiple outputs for different voltage levels and 
a choice of positive or negative voltage for the output. So, during 
this section, we will add an inner peak current controller, which 
will regulate the peak MOSFET current and then also a voltage 
controller to regulate the 5-V winding. As we mentioned, were 
going to present these control analysis tools within the context of 
a design example. So, for the current and voltage control loop 
design of the flyback converter, we can use the procedures 
depicted in Figure 4.45. 

This section will use a design flow that works well within the 
PLECS modeling environment using these tools. There are of 
course other approaches and methodologies that are possible. 
So as a starting point, we will determine the small signal trans- 
fer function of the current-controlled converter in the form of 
generating a Bode plot, and this will then allow us to design a 
voltage controller and then subsequently calculate the loop gain 


C: 47e-6 
v_init: 300 


R:0.5 | | —-C5 


Isense 


FIGURE 4.44 The topology of auxiliary-based power supply. 
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Current control loop design 
Converter open loop 
transfer function 


Gain/phase margins met? 


Verify with simulation 


FIGURE 4.45 The current and voltage control loop design of the flyback 
converter. 


response to ascertain if that system is now stable. So as a quick 
review here, let's go over the peak current controller block, 
which is actually found in the PLECS library as illustrated in 
Figure 4.46. 

So, this is a modifiable controller. It's part of the controls cate- 
gory. Well be opening it in PLECS, allowing you to see how it is 
based on an SR latch or flip-flop storage circuit. It has also a slope 
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Peak current controller 


FIGURE 4.46 ‘The peak current controller. 


control to avoid oscillations when duty cycle, d, is greater than 
5096 (d » 5096). At the beginning of each switching cycle, a pulse 
is applied at the input at the set input S of the flip-flop, which also 
turns the switch on. The switch signal remains on until the maxi- 
mum duty cycle as specified by the user in the controller's param- 
eters is reached or the current sent signal exceeds the current 
reference signal. The slope control is also included in the peak 
current mode controller to avoid oscillations if the duty cycle is 
greater than 5096 and it does this by subtracting a ramp signal 
from the reference current input. Lets get to work now. We'll 
open up PLECS as we mentioned standalone and you should now 
see a model in front of you in the PLECS schematic as shown in 
Figure 4.47. 

So, here in Figure 4.47, we do have a model of that flyback con- 
verter with no closed-loop control at all. The gate of our MOSFET 
is connected to a pulse generator block and we'll first run an open- 
loop simulation here so you can see the output voltage waveform. 
We'll do this by going to simulation and then select the start option, 


D6 
vf: 0.3 
Ron: 0.01 


Cc 
C: 47e-6 
v_init: 300 


D4 


DutyCyde: 0.01 


FO aww 


Vin 


Input Voltage 
m 
MOSFET MOSFET Signals 


FIGURE 4.47 The flyback converter model in PLECS schematic. 
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FIGURE 4.48 The voltage waveform of +5 V scope. 


and as shown in Figure 4.48 that the startup transients are really 
not very smooth. Hence, we need to add a voltage controller. 

This is what we're going to be working on. In the PLECS library 
browser, were going to add components in our library to our 
schematic and to implement our control architecture. So first, we 
should replace the pulse generator block with the peak current 
controller that we've been talking about. So, we can use the search 
function and search for the component and then we drag it into 
our schematic. First, we'll connect the switch signal to the 
MOSFET, and now we can right click this component as with all 
components in the library that are subsystems. You can select the 
look under mask option as shown in Figure 4.49, and you can see 
the implementation as in Figure 4.46, which we've modeled it and 
you can then modify it as you would for your own application. 

So, we can see here the behavior described in Figure 4.46. If we 
double click the peak current controller, we can see the masked 
subsystem parameters, so the only default parameters here that 
we're going to change are just the duty cycle, which we will set to 


Cut 


Con copy 
Delete 


Comment out 


Comment through 
Format 


Create subsystem 


FIGURE 4.49 Selecting the look under mask option. 
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Ctrl+/ 
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Peak Current Controller (mask) (link) 


Used for implementing current mode control in a switching 
converter. At the beginning of each switching cyde, the 
output is set. When the Isense input exceeds the Iref input, 
the output is reset. 


Parameters 
Switching frequency [Hz]: 


Slope compensation [A]: 
lo 


Pe | ea lee 


FIGURE 4.50 Setting the peak current controller parameters. 


5096. We'll check the boxes, which will display the parameters on 
the schematic for us to see as depicted in Figure 4.50. 

There we go, so next, what we want to do is we want to create a 
current sense signal to feed into our I ense of the peak current con- 
troller. So, we're going to do that by simply placing a resistor com- 
ponent, which we'll grab from our passive electrical library and 
place that into our schematic, and we'll wire it up. We will change 
the resistor value to 0.5 ohms and we can display that. Now, we 
also need a voltmeter component. Again, utilizing the search 
function and we wire it up as illustrated in Figure 4.51. 

So, here we go, and we've routed now our Lense signal, and now 
for the reference current signal, I,. we're actually going to first use 
a constant block, so we're going to specify a set point value here, 
0.1 amps as you can see in Figure 4.51. This is going to limit our 
actual peak inductor current to 10 amps and it's going to create an 


FIGURE 4.51 


Wiring up the peak current controller signals. 
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File Edit View Window Help 


OOQ- m g&o 


FIGURE 4.52 Regulation on 5 V with peak current controller. 


output voltage of 5 V in other primary winding of flyback trans- 
former. So, we can run this simulation again, and now you can see 
the regulation on the 5-V output as shown in Figure 4.52. 

So, what’s impractical about the controller that we just imple- 
mented? So, what's impractical about the converter system is that 
an outer, not typically in a system, we're going to have an outer 
voltage loop that provides this reference current, 0.1 A, for the 
current controller in order to regulate the output. To design the 
voltage controller, now we'e going to use the small signal transfer 
function of the current-controlled converter to do that. We can 
also now get into these impulse response and AC sweep tools to 
get this function in the form of a Bode plot. So, to do that, we're 
going to open up another model as depicted in Figure 4.53. 

So, we've pre-configured this model to perform this impulse 
response analysis. We will explain the changes that we've made. 
So, first things first, we have changed the input rectifier and 
capacitor, and we've changed it to a 300-V DC source component, 
the dynamics of the input rectifier are much slower than those of 
the flyback converter and can therefore be ignored. Simple simpli- 
fications like this don't affect our controller design, and they make 


Constant Peak Current 
Controller 
fs: S0e3 Probe 
dmin: 1 
E dmax: 50 MOSFET 
iref Islope: 0 


Gain 


FIGURE 4.53 The impulse response analysis schematic. 
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this small signal analysis easier to execute. To obtain a Bode plot 
of the open voltage loop, we've added two special components to 
the model. So, you can see here, we've added what's called a small 
signal perturbation block (i,,'). So, this component is added to 
our current set point value. We also need to add a small signal 
response to the perturbation to measure the perturbation at the 
5-V output. So here's our small signal response block (vo') that 
we're using to measure the output. These components are found in 
the controls category, and then found down there in the small sig- 
nal analysis sub-category of the controls category. We've labeled 
these as i,,,/ and vo’, respectively, and now we can open the con- 
figuration settings ofthe analysis by going to the simulation menu 
and going to analysis tools, and then click on the plus (+) button 
to add an impulse response analysis as shown in Figure 4.54. 

So, there are a few settings to highlight as shown in Figure 4.55. 

We've got the impulse response analysis selected as our choice of 
analysis. We've labeled it simply as impulse response analysis. This is 
a configuration you add and remove analyses as you see fit for your 
model. The system here, we've specified a period length, and it's the 
least common multiple of all the periodic sources in the model, so 
which in our case is the switching signal or 1 over 50 kHz. The 


1] Analysis Tools: plecs2 

Analyses Analysis type: Impulse Response Analysis 

Impulse Response Analysis| Description: [Impulse Response Analysis 

Impulse R Analysi Ec EE 

pue Response Anas Setup | Options | Steady-State Options | 

Impulse Response Analysis 
mm P 
| + | ex _Show log . [Start analysis | Accept Revert g Help " " 


FIGURE 4.54 Adding an impulse response analysis. 


192 m Digital Control of Power Converters 


Analysis type: Impulse Response Analysis 
Description: — Impulse Response Analysis 


Setup | Options | Steady-State Options | 


System period: auto 


Frequency range: |[1 1e4] 
Amplitude: le-3 
Perturbation: ire? 


Response: vo" 


Start analysis Accept Revert 


FIGURE 4.55 The settings for impulse response analysis. 


frequency sweep range should be a vector that includes the lowest 
and the highest perturbation frequencies, and then on the perturba- 
tion and response section, we want to indicate the perturbation and 
response blocks that we have that we want to use for this analysis in 
our model. We only have one of each, so it's very easy for us as illus- 
trated in Figure 4.55. In the options tab as shown in Figure 4.56, we 
can specify the number of automatically distributed frequencies. 


Analysis type: Impulse Response Analysis 


[RENE re rt 


Setup | Options | Steady-State Options 


Simulation start time: {0 
Frequency scale: logarithmic 
Number of points: 21 


Start analysis Accept Revert 


FIGURE 4.56 The options tab of analysis tools window. 
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[7-1 Impulse Response Analysis > E 
File Edit View Window Help 
coana regeo 


10-— 
E 


Magnitude / dB. 


FIGURE 4.57 The Bode plot of impulse response analysis. 


So, were going to start this analysis, and then we're going to 
explain a few more things. So, that was very quick, we're done, 
that analysis was super-fast, and that's because the impulse 
response analysis is actually the fastest method that we have of 
generating a Bode plot as shown in Figure 4.57. 

When we run the impulse response analysis, it actually first 
runs a steady-state analysis in order to find the stable operating 
point ofthe system. The steady-state analysis is based on a Newton 
Raphson iterative technique, and it may not converge when the 
system is run from startup because the system is under damped, 
and the state variables in PLECS are going to be too far from the 
final operating point. So, the simplest method of addressing this 
problem is to simulate the system for a specified number of initial- 
ization cycles before the steady-state analysis is run. If we go back 
to our configuration settings and we go to the steady-state options, 
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Analyses Analysis type: Impulse Response Analysis 
Impulse Response Analysis | Description: [Impulse Response Analysis 


pupils Response Vine “Setup | Options | Steady-State Options 


Number of init. cydes: [10 


Termination tolerance: |1e-1 


Max. number of iterations: |20 


Rel. perturbation for Jacobian: [1e-1 


Jacobian calculation: [fast à 


Start analysis Accept Revert 


FIGURE 4.58 The steady-state options tab of analysis tools window. 


we can see here in Figure 4.58 that we have the ability to specify 
this number of initialization cycles. 

This then causes the Jacobian matrix to be calculated from 
an operating point that is closer to the steady-state operating 
point, which increases the likelihood of convergence in the 
under-damp system and thus allowing the analysis to run to 
completion. So, that's our analysis, so we can go back to this 
frequency response. We can see the frequency response here of 
our current controlled flyback converter. We're now going to 
save this data. So, we'll go file export as a CSV file as depicted 
in Figure 4.59. 

So, we'll save it as flyback.csv, so now at this stage, we are ready 
to start talking about that voltage controller. So, what we're going 
to use MATLAB for a PID controller design. So, in the MATLAB 
command window, we write the following codes. 


data = csvread('C:NCRC Press\flyback.csv’); 
w = 2*pi*data(:,1); 
val-10.^(data(:,2)/20).*exp(j*data(:,3) *pi/180); 
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| File | Edit View Window Help 


Load trace data... 


Save trace data... 


Page setup... Ctrl+Shift+P 
Print... Ctrl+P 


Export as PDF... 


as Bitmap... 


PLECS Preferences... 
All... 


PLECS Extensions... Cursor range... 


FIGURE 4.59 Exporting the Bode plot as a CSV file. 


sys = frd(val,w); 
bode (sys) ,grid minor 
G = tfest(sys,3) 


The estimated transfer function using MATAB with an accu- 
racy of 99.71% is as below. 


G= 
-4.024e04 s*2 + 5.83e09 s + 1.338e12 
s*3 + 2.596e05 s^2 - 4.657e08 s - 6.749e12 
Continuous-time identified transfer function. 
Parameterization: 
Number of poles: 3 Number of zeros: 2 
Number of free coefficients: 6 
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Use "tfdata", "getpvec", "getcov" for 
parameters and their uncertainties. 
Status: 
Estimated using TFEST on frequency response data 
"Sys". 


Fit to estimation data: 99.71% 
FPE: 1.067e-05, MSE: 8e-06 


Then, we can write the following code. 
pidTuner (G) 


We can adjust the response time and transient behavior slides 
to get the step response as illustrated in Figure 4.60. 

Then, we click on the Show Parameters button to see the PID 
controller parameters alternatively, we can click on the Export 
button and give it a name (Cg) and click on the OK button as 
shown in Figure 4.61. 


[ Pommer 
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Controller Parameters: Kp = 2.217, Ki = 4.784e«04, Kd = 1.226e-05 . 


FIGURE 4.60 The PID tuner window. 
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FIGURE 4.61 The show parameters and export windows. 


Next, in the command window of MATLAB, we can type the 
‘Cg’ and press the enter key to get the PID controller parameters as 
below. 


>> Cg 
Cg = 
1 
Kp + Ki * --- + Kd * 5s 
S 
with Kp = 2.22, Ki 4.78e+04, Kd = 1.23e-05 
Continuous-time PID controller in parallel 
form. 


Now, we can add the Continuous PID Controller block in 
Figure 4.53 to create a closed-loop voltage and peak current con- 
troller as depicted in Figure 4.62. 

The waveform of 5 V scope is illustrated in Figure 4.63. 


Constant Gain Peak Current 
Continuous PID Controller "E 
d E 
ki: 4.78e4 MOSFET MOSFET Signals 


kd: 1.23e-5 


FIGURE 4.62 The closed-loop voltage and peak current controller. 


SJ919Au0O?) JAMO JO [ONUC [eSI] m 961 


Digital Control Implementation with PLECS m 199 


FIGURE 4.63 The waveform of 5 V scope. 


Then, we add a load step to the model, to evaluate the tran- 
sient behavior of closed-loop controller. So, at the time 25 ms, 
we add a parallel resistor (2 Q) in the 5 V circuit as shown in 
Figure 4.64. 

The output voltage waveforms during the load step at the time 
25 ms are shown in Figure 4.65. 

Therefore, our closed-loop controller is stable. In order to code 
generation for the closed-loop voltage and peak current controller, 
we use a subsystem as shown in Figure 4.66 and take our control- 
ler inside of it. Then we use STM32 target elements such as ADC 
and digital out as depicted in Figure 4.67. 

If we run the schematic of Figure 4.66, the output voltage wave- 
forms remain the same as Figure 4.65 as shown in Figure 4.68. 

Then, we do right click on the subsystem and select Subsystem 
and Execution settings... and in the popped-up window we tick 
the Enable code generation and choose the Discretization step 
size as 2e-7 as illustrated in Figures 4.69 and 4.70, respectively. 


THEN 5V/+12V/-12V 


MOSFET MOSFET Signals 


Continuous PID 
Controller 
kp: 2.22 
ki: 4.78e4 
kd: 1.23e-5 


FIGURE 4.64 Adding a load step to the model. 
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File Edit View Window Help 
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FIGURE 4.65 The output voltage waveforms during the load step. 


In the Open projects from File System window, we click on the 
Archive... button as illustrated in Figure 4.29, to open the g474. 
zip file in the following directory as shown in Figure 4.30. 


PLECS4.7(64 bit) > tsp stm32 > projects > g474. 
zip 


Then we run the STM32CubeIDE to enter it and from File 
icon we select Open projects from File System... option as 
depicted in Figure 4.28. Then we do click on the Finish button in 
the Import Projects from File System or Archive window as 
shown in Figure 4.31. In the Project Explorer, we find the cg 
folder and then right click on it and select the Properties option 
as depicted in Figure 4.32. Therefore, the Properties for cg win- 
dow is opened and we do copy the Location directory as illus- 
trated in Figure 4.71. 


FIGURE 4.66 Taking our controller inside the subsystem. 
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Input: [1] 


Peak Current Digital Out 
Controller Port: B 
Pin: [0] 

OutputType: push-pull 


e PID() v 
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Controller 
kp: 2.22 
ki: 4.78e4 
kd: 1.23e-5 


Input: [3] 


V 
Input: [2] 


FIGURE 4.67 The subsystem controller and STM32 target elements. 
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FIGURE 4.68 The output voltage waveforms with STM32 target elements. 
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FIGURE 4.69 Selecting the execution settings... option. 
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Execution Settings 


[V] Treat as atomic unit 


[_] Minimize occurrence of algebraic loops 
Sample time: 


Code Generation 


[V] Enable code generation 


Discretization step size: 
Simulation mode: 


Open coder options 


2e-7 
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Apply 


FIGURE 4.70 Enabling code generation and setting discretization step 
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FIGURE 4.71 Doing copy the location directory. 


Apply and Close Cancel | 
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System General | Parameterinüning | Target | Scheduling | External Mode 
plecs pcc. pid3 


Subsystem Target: STM32G4x 


General | Task Scheduler | Protection | External Mode 


Chip: G474xx 


Chip variant: [Garin 
System dock frequency (SYSCLK) [MHz]: |170 

[V] Use internal oscilator 

External ADC reference (Vref+) M]: — |3.3 


Step size tolerance: [Enforce exact value 
Build type: [Generate code into STM32CubeIDE project 
STM32CUbelIDE project directory: |C:\plecs\stm32\erc\g474.zip_expanded\g474\cg 


Generate pinmap file: Disabled 


FIGURE 4.72 Settings of the coder options window. 


Then, we close the Properties for cg window and come back to 
the PLECS environment and in the Coder icon of menu we select 
the Coder options..., to open the Coder Options window and 
in the Target section we select the Target as STM32G4x and chip 
as G474xx, the Build type as generate code into STM32CubeIDE 
project and also we do paste the copied Location directory as in 
Figure 4.71 to the STM32CubeIDE project directory and also tick 
use internal oscillator and finally click on the Accept and Build 
buttons respectively as shown in Figure 4.72. 

The build is done successfully and we come back to the 
STM32CubeIDE and in the Project Explorer section we do 
right click and select the Refresh option to see the codes added 
to the cg folder, then we click on the Build All (Ctrl+B) button 
to build the project as shown in Figure 4.73. The build is fin- 
ished with zero errors and finally, we can run the program to 
upload the codes to the STM32 microcontroller as illustrated in 
Figure 4.74. 
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FIGURE 4.73 Clicking on the build all (Ctrl+B) button. 
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FIGURE 4.74 Successfully building and running the program. 
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crossover frequency, 108 
current and voltage controller 
design, 110 
current and voltage controller 
transfer functions, 
110-111 
gain and phase margins, 
111-117 
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parameters, 107 
small signal output current, 
108-109 
Z-domain transfer functions, 
113 
crossover frequency, 91-93, 95 
digital control implementation 
average current mode control 
of, 139-149 
firmware implementation of 
current control, 134-139 
firmware implementation of 
voltage control, 
129-134, 131 
firmware peripherals 
initialization with CubeMX 
software, 121-124 
hardware, 118-121, 119, 120 
manual firmware initialization, 
124-129, 125, 129 
duty cycle and inductance values, 82 
inductor current ripple, 82-83 
MATLAB, current/voltage control 
loop 
Bode plot, 96-103, 96-105 
c2d function in, 86, 87, 89, 89 
Google Colab code, 95 
PID tuner window, 84-86, 85, 
87, 88 
step response of, 105, 106 
transient response of, 105, 105 
parameters of, 81, 82, 84 
schematic of, 81, 82 
small signal transfer function, 83 
switching frequency, 90 
voltage and current sensors, 81 
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CCM, see continuous conduction 
mode (CCM) 

closed-loop controller in PLECS and 
MATLAB 


analysis tools selection, 154, 156 

analysis tools window settings, 
154, 156 

Bode diagram, MATLAB, 157, 159 

choose name and saving the CSV 
file, 156, 158 

continuous PID controller block, 
162, 162 

control to output frequency 
response, 154, 157 

designed buck converter, 151, 152 

enter a name PID controller 
window, 160, 161 

first row of data.csv. remove, 156, 
159 

frequency response data, CSV file, 
154, 157 

MATLAB command window, 160 

output voltage waveform, 153, 153 

output voltage waveform, scope, 
162, 165 

PID(s) block to the closed-loop 
control system, 162, 164 

PID(s) transfer function values 
setting, 162, 163 

PID tuner window setting, 160, 
161 

simulation parameters window, 
153, 154 

small signal analysis elements, 
153, 155 

upper and lower saturation limits, 
anti-windup tab setting, 
162, 163 

code generation for the STM32 

Nucleo-G474RE 

ADCI block parameters window, 
168, 169 

buck converter and PID controller 
creation, 166, 166 

build finished message, 178, 178 

circuit inside the Buck_Converter 
subsystem, 166, 167 


circuit inside the PID_Controller 
subsystem, 166, 168, 168 

clicking on the archive, 172 

click the build all (Ctrl+B) button, 
178, 178 

click the finish button, 172, 173 

coder options selection, 174, 175 

coder options window setting, 
174, 176 

codes added to the cg folder, 176, 
178 

copy the location directory, 174, 
175 

enable code generation, 170, 171 

error message, 174, 176 

open projects from file system 
selection, 172 

open the g474.zip file, 172, 173 

PID Controller subsystem 

e vandm ports remove, 175, 
177 
without e v and m ports, 175, 

177 

program run, 178, 179 

properties option selection, 174, 
174 

PWM block parameters window, 
168, 169 

refresh option selection, 176, 177 

results of simulation, 168, 170 

start option selection, 168, 170 

steps, 165-166 


continuous conduction mode (CCM), 


90, 139, 145, 147 
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DC motor speed control 


Arduino demonstration 
Arduino script, 70-75 
DC gearmotor’s, Serial Plotter, 
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Measured Speed, 77, 77 

set point speed tracking, 70, 70 

simulated data, 78, 79 

Simulink block diagram, 
closed-loop feedback 
system, 69, 70 

theoretical result and measured 
result, 78, 80, 80 

controller design 

approximate transfer function 
for DC gearmotor, 60, 61 

classical control techniques, 66 

closed-loop transfer function, 
63 

continuous-time block 
diagram, 63, 63 

continuous-time controller, 
68, 68 

difference equation, 69 

explicit transfer functions, 64, 
64 

rational transfer functions, 
65, 65 

sampling time, T, 61, 62, 63 

Show Parameters tab, 67, 67 

transfer functions P(s) and 
Gzouy(s), 66 

Tustin’s method, 68 

system modeling 

Arduino script, 55-58 

block diagram, control loop, 
49, 49 

data variables, 58, 58 

Kirchhoff’s voltage law, 52 

Newton’s second law, 52 

open-loop transfer function, 52 

parameter values, 58-59, 59 

plant’s transfer function, 60 

procedures, 50-51 

setup, schematic of, 49, 49 

Simulink model, 54, 54-55 

transfer function, 51, 51-52, 53 
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average current mode control 
of, 139-149 

firmware implementation of 
current control, 134-139 

firmware implementation of 
voltage control, 129-134, 
131 

firmware peripherals 


initialization with CubeMX 


software, 121-124 
hardware, 118-121, 119, 120 


manual firmware initialization, 


124-129, 125, 129 
with PLECS, see PLECS software 
digital controller design 
approximate sampling time, T, 
41, 48 
bilinear transform/Tustin's 
method, 48 
implementing, 49 
inverse Z-transform, 48 
transfer function as P(s), 40-42 
ZOH/sampling, 42, 43, 44-45, 
45-46 
digital control with Arduino 
advantages of 
analog/continuous-time 
control system, 3, 3 
buck converter, 3, 3-4 
digital control scheme, 6, 6 
type II compensator, 4, 4 
DC motor speed control, see DC 
motor speed control 
digital controller design 
approximate sampling time, T, 
41, 48 
bilinear transform/Tustin's 
method, 48 
implementing, 49 
inverse Z-transform, 48 
transfer function as P(s), 
40-42 


ZOH/sampling, 42, 43, 44-45, 
45-46 
digital feedback control algorithm 
control structure, 38, 38 
“if? statement, 37 
implementation, 35, 35-36 
impulse function, 38, 38 
ZOH block operation, 37, 
37-38 
disadvantages of analog/ 
continuous-time 
controllers, 5 
discrete-time system 
Arduino IDE’s Serial Plotter, 
29, 29 
data variable creation, 32, 32 
frequency response of, 20, 21 
frequency responses, 22, 23, 24 
From Workspace block 
properties, 32, 33 
if statement, 27-28 
input step signal, 26, 26 
inverse Z-transform, 25 
MATLAB, 19-20 
rename the data variable as 
Output Data, 32, 33 
scope and running the 
simulation, 30, 31 
second-order continuous-time 
system, 19, 20 
second-order system, 24 
setting, 30, 30 
step input block, 30, 30 
step response in scope, 31, 31 
theoretical simulated output 
against the implementation 
output, 33, 34 
transient response of, 21-22, 22 


discretization of continuous-time 


systems 
Bode diagram of, 18, 19 
frequency response analysis 
technique, 15 


MATLAB, 18 
natural logarithm, 15-16 
root locus, 14 
Tustins method, 14, 16-17, 17 
Z-transform, 10 
closed-loop feedback diagram, 
13 
difference equation, 14 
Kirchhoff's Current Law 
(KCL), 10 
Laplace transform, 8, 8-11, 11 
s-plane and z-plane, 11, 12 
type II compensator, 4, 9 
digital feedback control algorithm 
control structure, 38, 38 
“if? statement, 37 
implementation, 35, 35-36 
impulse function, 38, 38 
ZOH block operation, 37, 37-38 
double loop control, 139, 141 
duty cycle method, 140 
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explicit transfer functions, 64, 64 
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firmware 
implementation 
of current control, 134-139 
of voltage control, 129-134, 131 
manual initialization, 124-129, 125 
peripheral initialization with 
CubeMX, 121-124 
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Google Colab code, 95 
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Integrated Development Environment 
(IDE), 18 

inverse Z-transform, 16, 25, 48, 68-69 
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Kirchhoff's Current Law (KCL), 10 
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MATLAB, 18-19, 22, 58, 61, 157 
Bode diagram of, 159 
‘bode function, 20 
c2d function, 68, 86, 89 
closed-loop controller in, see 
closed-loop controller in 
PLECS and MATLAB 
closed-loop PID controller, 
160-161 
control design tools, 64, 66 
current control loop, 84 
frequency responses, 24 
measured output speed, 77, 77 
Octave, 108 
phase margin-based frequency 
response, 95 
PI controllers, pidTuner tool, 95, 
97, 101, 105, 105, 106, 113 
PID controller design, 194, 197 
rational transfer function, 65 
two-dimensional variable, 31 
MATLAB Simulink, 30, 54-55, 77, 105 
Metal-Oxide-Semiconductor Field- 
Effect Transistor (MOSFET), 
108, 118-120, 139-140, 153, 
162, 180, 183, 185 
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Newton’s second law, 52 
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Octave software, PI controller design 
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crossover frequency, 108 
current and voltage controller 


design, 110 


current and voltage controller 


transfer functions, 110-111 


gain and phase margins, 


111-117 


parameters, 107 
small signal output current, 


108-109 


z-domain transfer functions, 113 


Pade approximation, 64 
PI controller design, see proportional 


integral (PI) controller 
design 


pidTuner(sys) function, 66 
PLECS software 


closed-loop controller 


analysis tools selection, 154, 
156 

analysis tools window settings, 
154, 156 

Bode diagram, MATLAB, 157, 
159 

choose name and saving the 
CSV file, 156, 158 

continuous PID controller 
block, 162, 162 

control to output frequency 
response, 154, 157 

designed buck converter, 151, 
152 

enter a name PID controller 
window, 160, 161 

first row of data.csv. remove, 
156, 159 

frequency response data, CSV 
file, 154, 157 


MATLAB command window, 
160 

output voltage waveform, 153, 
153 

output voltage waveform, 
scope, 162, 165 

PID(s) block to the closed-loop 
control system, 162, 164 

PID(s) transfer function values 
setting, 162, 163 

PID tuner window setting, 160, 
161 

simulation parameters window, 
153, 154 

small signal analysis elements, 
153, 155 

upper and lower saturation 
limits, anti-windup tab 
setting, 162, 163 


code generation for the STM32 


Nucleo-G474RE, 165-179 

ADCI block parameters 
window, 168, 169 

buck converter and PID 
controller creation, 166, 166 

build finished message, 178, 
178 

circuit inside the Buck . 
Converter subsystem, 166, 
167 

circuit inside the PID_ 
Controller subsystem, 166, 
168, 168 

clicking on the archive, 172 

click the build all (Ctrl+B) 
button, 178, 178 

click the finish button, 172, 173 

coder options selection, 174, 
175 

coder options window setting, 
174, 176 

codes added to the cg folder, 
176, 178 


copy the location directory, 
174, 175 

enable code generation, 170, 171 

error message, 174, 176 

open projects from file system 
selection, 172 

open the g474.zip file, 172, 173 

PID Controller subsystem, 
175, 177 

program run, 178, 179 

properties option selection, 
174, 174 

PWM block parameters 
window, 168, 169 

refresh option selection, 176, 
177 

results of simulation, 168, 170 

start option selection, 168, 170 

steps, 165-166 


features of, 150-151 
small signal analysis tools 


add an impulse response 
analysis, 191, 191 

Bode plot of impulse response 
analysis, 193, 193 

dick the build all (Ctrl+B) 
button, 206, 207 

closed-loop voltage and peak 
current controller, 197, 198 

coder options window setting, 
206, 206 

controller inside the subsystem, 
199, 202 

copy the location directory, 
205, 206 

current and voltage control 
loop design, flyback 
converter, 180, 182 

enable code generation and 
setting discretization step 
size, 199, 204 

execution settings selection, 
199, 204 
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file export as a CSV file, 194, 
195 

flyback converter model, 183, 
184 

impulse response analysis 
schematic, 189, 190 

impulse response analysis 
setting, 191, 192 

look under mask option 
selection, 185, 186 

options tab of analysis tools 
window, 192, 192 

output voltage waveforms 
during the load step, 199, 
201 

output voltage waveforms with 
STM32 target elements, 
199, 203 

peak current controller, 182, 
183 

peak current controller 
parameters setting, 187, 187 

PID tuner window, 196, 196 

regulation on 5 V with peak 
current controller, 189, 189 

run the program, 206, 207 

show parameters and export 
windows, 196, 197 

steady-state options tab of 
analysis tools window, 194, 
194 

topology of auxiliary-based 
power supply, 180, 181 

voltage waveform of +5 V 
scope, 185, 185 

waveform of 5 V scope, 197, 
199 

wiring the peak current 
controller signals, 187, 188 


proportional integral (PI) controller 
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Bode plot, 96-103, 96-105 

c2d function in, 86, 87, 89, 89 

vs. current and voltage control 
loops, 105, 105 

Google Colab code, 95 

phase margin calculations, 105 

PID tuner window, 84-86, 85, 
87, 88 

step response of, 105, 106 

Octave software 

crossover frequency, 108 

current and voltage controller 
design, 110 

current and voltage controller 
transfer functions, 110-111 

gain and phase margins, 
111-117 

parameters, 107 

small signal output current, 
108-109 

Z-domain transfer functions, 
113 


Serial Plotter, 29, 29, 76, 76 
small signal analysis tools, PLECS for 


controller design 

add an impulse response analysis, 
191, 191 

Bode plot of impulse response 
analysis, 193, 193 

click the build all (Ctrl+B) button, 
206, 207 

closed-loop voltage and peak 
current controller, 197, 198 

coder options window setting, 
206, 206 

controller inside the subsystem, 
199, 202 

copy the location directory, 205, 
206 


current and voltage control loop 
design, flyback converter, 
180, 182 

enable code generation and setting 
discretization step size, 199, 
204 

execution settings selection, 199, 
204 

file export as a CSV file, 194, 195 

flyback converter model, 183, 184 

impulse response analysis 
schematic, 189, 190 

impulse response analysis setting, 
191, 192 

look under mask option selection, 
185, 186 

options tab of analysis tools 
window, 192, 192 

output voltage waveforms during 
the load step, 199, 201 

output voltage waveforms with 
STM32 target elements, 
199, 203 

peak current controller, 182, 183 

peak current controller parameters 
setting, 187, 187 

PID tuner window, 196, 196 

regulation on 5 V with peak 
current controller, 189, 189 

run the program, 206, 207 

show parameters and export 
windows, 196, 197 

steady-state options tab of analysis 
tools window, 194, 194 

topology of auxiliary-based power 
supply, 180, 181 

voltage waveform of +5 V scope, 
185, 185 

waveform of 5 V scope, 197, 199 

wiring the peak current controller 
signals, 187, 188 
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STM32CubeIDE, 121, 123-124, 166, 
172, 174, 176, 201, 206 
STM32CubeMX, 121 
STM32 microcontroller ADC pin, 127 
STM32 Nucleo-G474RE, 118-121; see 
also code generation for the 
STM32 Nucleo-G474RE 
system modeling, DC motor speed 
control 
Arduino script, 55-58 
block diagram, control loop, 49, 49 
data variables, 58, 58 
Kirchhoff’s voltage law, 52 
Newton’s second law, 52 
open-loop transfer function, 52 
parameter values, 58-59, 59 
plant’s transfer function, 60 
procedures, 50-51 
setup, schematic of, 49, 49 
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Z-transform, 10 
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difference equation, 14 
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s-plane and z-plane, 11, 12 
type II compensator, 4, 9 


